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ABSTRACT 

A  personalized  system  of  instruction  combining  the 
principle  of  auto-tutorial  instruction  with  modified  self- 
pacing  was  applied  to  a  course  in  the  fundamentals  of  gas 
dynamics  taught  to  thirteen  students  in  the  Department  of 
Aeronautics  at  the  Naval  Postgraduate  School,  Monterey, 
during  the  Spring  quarter  of  1972.   The  course  results  are 
summarized  and  the  instructional  materials  developed  for  the 
course  are  included  as  appendices.   While  some  degradation 
was  apparent  due  to  modifications  made  to  the  method,  the 
results  tended  to  affirm  the  substantial  value  of  this 
instructional  method  and  its  superiority  to  the  conventional 
lecture  method. 


TABLE  OF  CONTENTS 

I.  THE  EDUCATION  PROBLEM  5 

II.  THE  CONVENTIONAL  LECTURE  METHOD  9 

III.  ALTERNATIVES  TO  THE  CONVENTIONAL 

LECTURE  METHOD  11 

IV.  PERSONALIZED  INSTRUCTION  13 

V.  APPLICATION  OP  PSI  TO  GAS  DYNAMICS 16 

VI.  EVALUATION  OF  RESULTS 17 

VII.  CONCLUSIONS  AND  RECOMMENDATIONS  18 

APPENDIX  A   STUDY  GUIDES  21 

APPENDIX  B   STUDENT  EVALUATION  QUESTIONNAIRE  176 

APPENDIX  C   SELECTED  COMMENTS  AND  CRITIQUE  l8l 

LIST  OF  REFERENCES 182 

INITIAL  DISTRIBUTION  LIST lQl\ 

FORM  DD  1473 185 


ACKNOWLEDGEMENTS 

The  author  wishes  to  express  his  appreciation  to 
Associate  Professor  Robert  D.  Zucker  for  his  patience  and 
guidance  as  advisor  to  the  author  in  the  development  of 
this  thesis. 


I.   THE  EDUCATION  PROBLEM 

Faced  with  the  spectre  of  the  "post-industrial  age" 
[Ref.  1]  with  its  exponential  rise  in  population  and  complexity 
of  life,  the  engineering  profession  is  being  tasked  to  solve 
problems  of  vastly  increased  scope  with  reduced  assets. 
Rising  costs,  increasingly  critical  public  scrutiny,  and  a 
disturbing  downward  trend  in  the  number  of  candidates  for 
engineering  training  are  indicative  of  an  impending  crisis  in 
technical  education.   Some  of  the  more  astute  observers, 
including  Vice  Admiral  Hyman  Rickover,  have  been  forecasting 
such  problems  for  more  than  a  decade  [Ref.  2].   Recognizing 
these  trends,  engineering  educators  are  urgently  exploring 
new  methods  of  instruction  which  will  not  only  provide  more 
depth  to  the  engineering  student's  education,  but  will  also 
provide  for  more  efficient  time  utilization  by  both  students 
and  teachers. 

R.  Buckminster  Fuller  maintains  that  it  will  not  suffice 
to  merely  expand  our  present  educational  systems,  but  rather 
that  a  "transformation"  of  great  magnitude  is  essential  [Ref. 
3].   The  early  stages  of  this  transformation  are  already  in 
progress  and  have  provided  some  of  the  impetus  and  direction 
to  this  project. 

Innovations  in  engineering  education  almost  invariably 
represent  a  departure  from  the  venerable  and  currently  preva- 
lent mode  of  teaching  known  as  "lecturing,"  wherein  the 


instructor  devotes  a  major  portion  of  his  contact  hours  with 
the  students  to  the  role  described  by  Professor  Flammer  of 
Utah  State  as  that  of  a  "dispenser  of  information"  [Ref.  4], 
He  sees  a  dire  need  for  the  role  of  the  professor  to  change  — 
not  diminishing  as  some  observers  fear  with  the  advent  of 
various  other  media  for  the  dispensing  of  information,  but 
rather  assuming  vastly  greater  importance  as  the  teacher 
becomes  more  of  a  diagnostician  and  prescriber  of  learning 
experiences  for  each  individual  student  [Ref.  4],   As  Dr. 
Billy  Koen  has  suggested,  "...  it  is  not  unreasonable  to 
expect  that  engineering  educators  ought  to  be  able  to 
'engineer'  their  courses."  [Ref.  5] 

Learning  may  be  thought  of  as  resulting  from  the  auspi- 
cious combination  of  several  interactive  factors  which  can 
either  promote  or  hamper  successful  learning.   The  three 
principal  influences  are  motivation,  operant  span  and  stimulus 
control .   Each  of  these  is  in  turn  influenced  to  a  greater  or 
lesser  degree  by  several  other  educational  factors.   To  draw 
an  analogy  from  the  realm  of  control  theory,  the  problem  in 
education  is  to  introduce  methods  which  generate  feedback  in 
order  to  modify  the  influence  of  each  of  these  independent 
factors,  so  as  to  optimize  learning.   One  basic  premise  of 
the  concept  of  'personalized  instruction'  is  that  this  more 
effectively  provides  the  required  feedback. 

The  influencing  factors  are  defined  as  follows: 

Motivation  -  The  force  which  drives  a  student  to  seek  to 
learn. 


Stimulus  Control  -  The  effect  exerted  by  the  learning 

surrounds  on  the  individual. 

Operant  Span  -  The  maximum  scope  (in  size  or  complexity) 

of  a  single  block  of  material  which  may  be 
assimilated  by  a  given  student  at  a  time. 

Structure  -  The  concept  and  method  of  presentation  of  the 
course  material. 

Facilities  -  The  physical  surrounds  (including  graphic  aids 
and  other  media). 

Pacing  -  Rate  of  delivery. 

Intellectual  Ability  -  The  student's  inherent  talent. 

Level  of  Development  -  The  student's  specific  level  of 

achievement  or  skill  in  a  given 
subject  area. 

Quotidian  Variability  -  A  psychometric  term  used  to  identify 

daily  variations  in  the  student,  not 
intrinsic  to  the  educational  method, 
such  as  illness,  fatigue,  or  behavior 
arising  from  idiosyncratic  response 
to  some  unique  aspect  or  feature  of 
the  educational  environs.  [Ref.  6] 

The  interaction  of  these  factors  forms  the  basis  of  a 
substantial  portion  of  current  educational  research.   Detailed 
investigation  of  these  interactions  represents  subject  matter 
well  beyond  the  scope  of  this  paper.   Structure  and  pacing 
are  the  most  obvious  areas  of  concern  in  a  system  of 
personalized  instruction,  though  each  factor  must  be  considered 

In  the  development  of  better  methods,  two  unique  aspects 
of  engineering  education  must  be  considered: 

a.  The  professional  skills  and  concepts  involve  principally 
the  cognitive  domain. 

b.  A  hierarchy  of  skills  exists,  which  causes  mastery  of 
most  areas  to  be  dependent  on  a  thorough  grounding  in 
preceding  material. 


Failure  to  thoroughly  master  the  required  hierarchy  of 
skills  (a  common  tendency  among  many  students)  can  lead  to 
the  development  of  a  condition  referred  to  as  "cumulative 
ignorance,"  which  probably  accounts  for  the  majority  of  the 
failures  and  "drop  outs"  in  the  engineering  educational 
process.   If  no  other  gain  be  realized,  the  possibility  of 
precluding  some  cumulative  ignorance  justifies  the  efforts 
expended  in  searching  for  better  methods.   However  there  is 
mounting  evidence  that  methods  have  been  devised  and  are 
proven  which  not  only  help  solve  the  problem  of  cumulative 
ignorance,  but  also  provide  far  more  optimal  interaction  of 
the  various  parameters  affecting  the  learning  process. 


II.   THE  CONVENTIONAL  LECTURE  METHOD 

The  standard  with  which  all  other  educational  methods  are 
usually  compared  (or  contrasted)  is  the  "conventional  'lecture 
method."   The  term  as  used  herein  connotes  a  formal  discourse 
wherein  the  instructor  ("lecturer")  presents  the  course 
material  through  the  media  of  the  spoken  word,  using  no 
materials  or  training  aids  other  than  a  chalk  board.   The 
lecture  method  normally  exhibits  the  following  characteristics: 

a.  The  lecture  is  time-constrained  by  the  standard 
scheduled  period. 

b.  The  student  is  in  a  passive  role  and  is  expected  to 
transcribe  by  notes,  or  otherwise  record,  whatever 
information  he  deems  essential. 

c.  The  pace  is  essentially  fixed  by  the  lecturer  based  on 
the  scheduled  contact  hours/course  coverage. 

d.  There  is  virtually  no  "real-time"  feedback.   Due  to  the 
time  constraints,  questions  or  individual  difficulties 
are  usually  deferred.   This  leads  to  the  not  unusual 
situation  of  the  student  missing  vital  succeeding 
concepts  while  momentarily  pondering  a  puzzling  point. 

e.  The  clarity  of  presentation  is  strongly  dependent  on 
the  expository  ability  of  the  lecturer. 

The  conventional  lecture  method  is  seldom  utilized  exactly 

as  described  above.   It  is  normally  modified  by  one  (or  more) 

of  many  variations  in  structure  or  facilities  ranging  from 

simple  devices,  such  as  the  use  of  colored  chalk,  to  the  complex, 

such  as  closed  circuit  television  with  printed,  pre-distributed 

notes.   Where  class  size  and  time  permit,  most  lecturers 

seek  feedback  by  encouraging  questions  from  the  class  at 

any  time.   While  these  variations  mitigate  somewhat  the 

shortcomings  of  the  lecture,  the  basic  weaknesses  remain. 


The  more  capable  students  are  held  back  and  the  students 
suffering  from  any  substantial  amount  of  "cumulative  ig- 
norance" find  their  burden  of  cumulative  ignorance  is  growing. 

The  lecture  method  has  three  advantages  which  have  acted 
to  foster  its  acceptance  for  hundreds,  if  not  thousands  of 
years.   On  closer  examination  these  'so-called'  advantages 
prove  to  contribute  to  administrative  expediency  -  but  often 
at  the  expense  of  learning.   There  is  truth  in  the  statement 
that  students  often  learn  the  material  of  a  course  in  spite 
of  the  teaching  methods  rather  than  because  of  them.   Firstly, 
the  conventional  lecture  provides  for  the  most  rapid  and 
inexpensive  output  of  information  by  the  lecturer  with  the 
minimum  amount  of  equipment  or  materials.   This  says  nothing 
about  the  amount  of  information  actually  received  by  the  stu- 
dent.  Secondly,  the  elegance  of  presentation  and  the  influence 
of  the  lecturer  on  his  class  can  serve  to  "sell"  the  concepts 
involved  to  the  class.   Thirdly,  this  method  is  accepted. 
This  is  a  very  real  consideration  when  one  must  deal  with 
the  problem  of  institutional  inertia. 
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III.   ALTERNATIVES  TO  THE  CONVENTIONAL  LECTURE  METHOD 

Many  of  the  variations  or  innovations  associated  with 
non-lecture  methods,  such  as  written  course  objectives, 
instructor  tutoring,  peer  tutoring,  visual  aids  and  printed 
notes  or  outlines,  can  be  effectively  used  to  enhance  the 
lecture.   However,  while  improving  the  lecture  method,  these 
innovations  do  not  directly  contribute  to  the  basic  need  of 
each  student  to  explore  the  topic  a_t  his  own  best  speed . 
This  can  only  be  achieved  by  abandoning  the  lecture  as  the 
primary  information  vehicle. 

Currently  in  vogue  among  engineering  educators  are 
several  varieties  of  non-lecture  methods,  nearly  all  based 
on  the  reinforcement  theories  first  expounded  by  Sidney 
Pressey  and  B.  F.  Skinner  [Ref.  7].   Their  research  resulted 
first  in  the  development  of  "programmed  instruction".   The 
Programmed  Text  and  the  Teaching  Machine  are  direct  applica- 
tions of  this  reinforcement  theory  and  are  now  widely  used 
(with  proven  effectiveness)  for  the  instruction  of  material 
involving  verbal  or  skill  knowledge.   These  methods  suffer 
from  relatively  inflexible  structure  and  the  tendency  to 
design  the  instructor  out  of  the  picture.   Also,  neither 
method  has  proven  quite  so  adaptable  to  complex  concepts 
involving  cognitive  strategy,  such  as  engineering  problem 
solving  or  design.   For  this  reason,  engineering  educators 
have  been  among  the  last  to  experiment  with  new  teaching 
techniques.   However,  these  shortcomings  have  recently 


1 1 


prompted  the  development  of  several  closely  related  methods 
known  variously  as  "Auto-Tutorial  Instruction"  [Ref.  8], 
"Self-Paced  Individually  Prescribed  Instruction"  [Ref.  9], 
"Personalized  (or  Proctorial)  System  of  Instruction" [Ref .  10, 
11]  or  some  similar  designation,  all  of  which  are  more  or 
less  self-paced  and  vary  primarily  in  the  method  by  which 
reinforcement  is  achieved  and  progress  is  evaluated.   However, 
in  all  cases  the  teacher,  in  the  role  of  diagnostician,  tutor, 
counselor  and  designer  of  learning  experiences  [Ref.  12]  is 
the  key  to  success  of  these  systems. 
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IV.   PERSONALIZED  INSTRUCTION 

Professor  F.  H.  Keller  is  generally  regarded  as  the 
principal  architect  of  this  method  which  he  calls  the 
"Proctorial  (or  personalized)  System  of  Instruction"  (PSI), 
though  many  of  the  features  were  concurrently  and  indepen- 
dently developed  by  Professor  S.  H.  Postlewaite  at  Purdue 
University  in  Indiana.   The  basic  features  of  this  method 
are  succinctly  described  by  Professor  B.  V.  Koen  of  the 
University  of  Texas  as  follows:  [Ref.  13] 

1.  The  ' go-at-your-own-pace '  feature,  which  permits  a 
student  to  move  through  the  course  at  a  speed  commensu- 
rate with  his  ability  and  other  demands  upon  his  time; 

2.  The  unit-perfection  requirement  for  advance,  which  lets 
the  student  go  ahead  to  new  material  only  after 
demonstrating  mastery  of  that  which  preceded; 

3.  The  use  of  lectures  as  vehicles  of  motivation,  rather 
than  sources  of  critical  information; 

4.  The  related  stress  upon  the  written  work  in  teacher- 
student  communication;  and,  finally: 

5.  The  use  of  proctors,  which  permits  repeated  testing, 
immediate  scoring,  almost  unavoidable  tutoring,  and 

a  marked  enhancement  of  the  personal-social  aspect  of 
the  educational  process. 

This  summary  by  Dr.  Koen  clearly  illustrates  the  particular 

desirability  and  applicability  of  this  method  to  engineering 

education  and  accounts  for  the  growing  success  of  this  method 

The  particularly  agonizing  problem  of  "cumulative  ignorance" 

is  effectively  overcome  by  the  unit-perfection  requirement  of 

this  method. 

The  concept  of  PSI  is  predicated  on  the  validity  of 

'reinforcement  theory . '   The  course  is  divided  into  an 
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appropriate  number  of  units,  depending  on  major  subject  areas 
and  complexity  of  the  material.   Each  student  is  provided 
with  a  Study  Guide  for  the  unit  of  instruction  on  which  he  is 
working  which  will  contain: 

a.  Course  objectives  -  defined  in  terms  of  measurable 
behavior. 

b.  Reading  assignments  (if  a  text  is  used). 

c.  Discussion  of  theoretical  concepts  and  application 
to  specific  situations. 

d.  Example  problems  with  a  discussion  of  the  general 
philosophy  and  specific  techniques  of  solution. 

e.  Practice  problems,  with  answers  and  guidance  where 
necessary. 

The  student's  initial  exposure  to  new  concepts  can  thus  proceed 

at  his  own  pace.   A  fundamental  concept  of  PSI  is  that  the 

instructor  holds  the  key  to  success  of  the  system  in  his  role 

as  "trouble-shooter."   Ideally  the  Study  Guides  will  provide 

a  sufficiently  clear  exposition  of  the  concepts  so  that  the 

teacher  will  be  freed  to  spend  most  of  his  time  dealing  with 

the  individual  problems  which  require  personal  attention. 

The  principal  problems  encountered  by  practitioners  of 

this  method  were  summarized  recently  by  Dr.  Gilmour  Sherman 

of  Georgetown  University  —  who,  as  an  associate  of  Professor 

Keller  in  the  development  of  the  PSI  method,  has  received 

considerable  feedback  from  various  educators  implementing  PSI 

in  one  or  another  form.   Dr.  Sherman's  list  can  be  summarized 

as  follows:   [Ref.  14] 

1.  Professor  works  harder  than  ever. 

2.  Students  work  harder. 

3.  Substantial  increase  in  logistic  and  administrative 
workload . 
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4.  Staff  (if  paid)  becomes  expensive. 

5.  Class  size  for  effective  implementation  is  limited. 

6.  Production  lead  time  for  materials  is  burdensome. 

7.  Procrastination  by  students  creates  problems  with 
the  administration. 

8.  How  to  assign  grades? 

Of  these,  the  last  two  problem  areas  seem  the  most  pernicious. 
The  advocates  of  PSI  with  experience  have  generally  shown 
that  the  system  can  be  engineered  to  minimize  all  of  these 
problems  with  a  few  iterations,  resulting  in  proven  increases 
in  effectiveness  of  instruction.   For  an  in-depth  current 
assessment  of  the  Keller  Method  of  PSI,  the  reader  is  referred 
to  Ref.  15,  in  its  entirety. 
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V.   APPLICATION  OF  PSI  TO  A  GAS  DYNAMICS  COURSE  AT  NPS 

An  initial  consideration  in  the  development  of  PSI 
materials  for  the  Gas  Dynamics  course  offered  at  NPS  was 
to  insure  compatibility  with  the  quarter  system.   It  was 
therefore  determined  that  while  retaining  the  self-paced, 
self-guided  features  of  the  course  for  each  instructional 
unit,  the  pace  would  be  controlled  in  a  longterm  sense  by 
scheduling  the  check  tests  (to  determine  mastery  of  each 
unit)  to  be  taken  at  regular  intervals  by  all  (13)  students 
involved  in  the  PSI  course.   This  also  permitted  adequate 
lead  time  for  the  development  of  the  Study  Guides  for  each 
unit  and  allowed  some  feedback  of  information  to  modify 
later  units  to  improve  any  notable  weaknesses  in  format 
discovered  in  early  units.   The  study  guides  are  included 
as  Appendix  A. 

Two  major  mid-course  modifications  were  instituted.   The 
first,  a  decision  to  discontinue  the  practice  of  repeating 
check  tests  (except  in  an  extreme  case)  was  made  by  the 
instructor  when  it  was  determined  that  time  constraints  may 
have  prevented  the  class  from  finishing  within  the  quarter. 
Secondly,  in  order  to  appease  a  number  of  students  who  found 
it  discomforting  to  have  no  group  sessions,  the  instructor 
instituted  the  practice  of  group  review  discussions  in  order 
to  re-inforce  the  material.   This  provided  tutoring  on  a 
group  basis,  for  those  who  desired  to  attend  the  problem 
sessions . 
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VI.   EVALUATION  OF  RESULTS 

The  success  of  this  project  was  intended  initially  to 
be  measured  by  three  principal  means: 

1.  Comparison  with  a  "control"  group  in  the  form  of  a 
class  of  twelve  students  taught  by  the  instructor 
using  the  conventional  lecture  method  during  the 
same  quarter. 

2.  A  comprehensive  student  questionaire  administered 
at  the  end  of  the  course  to  the  thirteen  students 
involved  in  the  "PSI"  class. 

3.  A  critique  session  concerning  the  course,  to  be 
moderated  by  the  writer  at  the  end  of  the  quarter. 

Due  to  the  lack  of  statistical  significance  of  any  data 
involving  groups  of  this  size,  grade  comparisons  are 
considered  to  have  little  meaning.   However  the  instructor 
felt  that  there  was  a  tendency  toward  better  overall 
comprehension  by  the  students  in  the  PSI  group  as  determined 
by  major  quiz  and  final  examination  averages. 

The  questionaire  (Appendix  B)  which  was  given  to  each 
student  in  the  PSI  group  revealed  several  key  facts.   Twelve 
of  the  thirteen  students  returned  the  questionaire  and 
where  appropriate  the  numerical  average  is  tabulated  in  the 
right  hand  column  beside  the  question.   Noteworthy  are  the 
student's  apparent  high  regard  for  the  concept  of  defined 
objectives,  study  guides  and  check  tests.   Significant 
comments  from  the  last  questions  are  included  in  Appendix  C. 

The  critique  session  tended  to  corroborate  the  comments 
in  the  questionaires .   A  noteworthy  comment  from  this  session 
is  included  in  Appendix  C. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

In-so-far  as  this  course  followed  the  model  of  PSI  as 
conceived  by  Professor  Keller,  the  results,  based  on  the 
subjective  evaluation  of  the  students  involved  and  the 
instructor,  tend  to  substantiate  the  value  of  this  approach. 
Given  the  bastardization  of  the  long-term  self-paced  concept 
which  resulted  from  controlling  the  administration  of  the  check 
tests  in  the  manner  described,  the  hoped-for  uniform  high 
level  of  mastery  was  not  achieved.   However,  several  comments 
by  students  in  the  course  indicated  the  'depth'  of  knowledge 
achieved  was  generally  greater  than  normal  and  most  students 
admitted  they  worked  harder  than  in  previous  courses.   The 
controls  applied  to  the  pacing  of  the  course  negated  a  true 
evaluation  of  the  problem  of  student  procrastination;  however, 
comments  from  two  students  who  fared  poorly  in  the  course 
attributed  this  to  a  personal  lack  of  interest  in  the  subject 
area.   It  is  also  noteworthy  that  only  four  of  the  twelve 
students  who  returned  questionaires  sought  tutorial  aid  from 
the  instructor.   This  resulted  in  part  from  the  students' 
dependence  on  the  review  discussion  sessions  to  iron  out 
problems  and  in  at  least  one  instance  was  the  admitted  result 
of  a  student's  personal  antipathy  toward  the  course  and 
instructor.   This  is  a  fairly  obvious  example  of  quotidian 
variability  adversely  affecting  motivation,  hence  learning. 
The  student  in  question  is  otherwise  average  but  fared  poorly 
in  this  course. 
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It  is  recommended  that  in  the  next  offering  of  Gas- 
Dynamics  the  course  be  made  truly  self-paced  with  the  use  of 
a  proctor  drawn  from  one  of  the  students  who  have  taken  the 
course  under  this  project.  Eleven  of  the  twelve  students 
who  responded  to  the  questionaire  indicated  their  willingness 
to  serve  as  a  tutor.   While  it  is  not  unreasonable  to  assume 
that  any  student  volunteering  to  serve  in  this  capacity  would 
be  substantially  compensated  by  the  challenge  involved  in 
this  task,  the  time  constraints  and  the  inherent  demands  of 
the  academic  load  certainly  justify  some  more  tangible  reward. 
Many  institutions  are  paying  the  graduate  students  utilized 
in  this  capacity.   This  being  inappropriate  at  this  institution, 
it  is  strongly  recommended  that  some  sort  of  academic  credit 
be  granted  to  the  student  undertaking  this  task  as  tutor  or 
proctor . 

V/ith  the  materials  and  experience  resulting  from  this 
first  cycle,  the  foundation  exists  to  truly  exploit  this 
proven  method  in  succeeding  offerings  of  the  Gas  Dynamics 
course  at  the  Naval  Post-graduate  school.   Moreover  the 
flexibility  and  effectiveness  of  this  method  indicate  the 
desirability  of  applying  this  method  to  other  courses. 
Perhaps  the  most  concise  statement  to  date  concerning  the 
underlying  strength  of  the  Keller  method  was  made  by 
Professor  A.  J.  Dressier  in  the  foreword  to  a  recent  work- 
shop on  the  Keller  Method.   [Ref.  15] 
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"The  rapid  spread  of  the  Keller  method  is  remarkable 
when  one  considers  the  handicaps  it  has  had  to  over- 
come : 

(1)  It  is  a  new  idea  that  challenges  a  teaching 
technique  (the  lecture  method)  that  is  hundreds,  if 
not  thousands,  of  years  old.   New  ideas  are  not  usually 
received  with  delight. 

(2)  It  is  based  on  a  controversial  theory;  the 
reinforcement  theory  associated  with  the  name  of 

B.  F.  Skinner. 

(3)  It  involves  more  work  for  both  the  teacher 
and  the  student.   More  work  is  something  that  most 
people  feel  they  do  not  need. 

In  spite  of  these  negative  features,  the  Keller  method 
has  spread  at  a  remarkable  rate  for  one  simple  reason  - 
it  works  I " 
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APPENDIX  A  -  STUDY  GUIDES 

The  Study  Guides  written  for  this  course  were  developed 
to  be  used  in  conjunction  with  the  text  Gas  Dynamics  by 
Professor  James  E.  A.  John  [Ref.  16],   Additional  reference 
material  and  assistance  was  gained  from  Refs.  17  -  20. 
The  Course  Objectives  were  developed  using  the  methods  in 
Refs.  21  and  22. 

Some  additional  handouts  were  prepared  by  the  instructor 
and  distributed  to  assist  the  students  with  certain  points. 
Some  of  this  material  will  eventually  be  incorporated  in 
the  future  revisions  of  the  Study  Guides. 
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INTRODUCTION 


It  is  assumed  that  before  entering  the  world  of  Gas 
Dynamics  you  have  had  a  reasonable  background  in  Mathematics 
(through  Calculus),  together  with  an  exposure  to  basic  Fluid 
Dynamics  and  a  course  in  elementary  Thermodynamics. 

Prior  to  undertaking  new  material  it  will  prove  fruitful 
to  review  some  fundamental  concepts  and  facts  from  which 
the  study  of  Gas  Dynamics  will  proceed. 

First  read  Sections  1.1  through  1.5  in  the  text.    You 
may  also  find  it  useful  to  review  your  notes  or  textbooks 
used  in  previous  Fluid  Mechanics  and  Thermodynamics  courses. 
Then  answer  the  following  review  questions  as  completely  as 
possible  using  any  source  you  wish.   Your  answers  are  to  be 
submitted  to  the  instructor  at  which  time  you  will  be  pro- 
vided correct  answers  to  these  review  questions.   No  attempt 
should  be  made  to  continue  with  this  course  until  you  feel 
that  you  understand  the  answers  to  all  of  these  review 
questions . 


All  references  to  sections,  equations,  diagrams,  etc.  in 
the  'Text'  refer  to  GAS  DYNAMICS  by  James  E.  A.  John.   Any 
other  referenced  sources  will  be  referred  to  by  name. 
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REVIEW  QUESTIONS 

1.  Can  you  recall  the  Taylor's  Series  Expansion,  it's 
applications  and  limitations? 

2.  How  is  a  derivative  such  as  ^  defined? 

3.  Cite  the  basic  units  of  measurement  in  the  MLT  system. 

(Metric  and  English) 

M.   What  is  the  definition  of  a  one  pound  force  (in  terms 
of  pounds-mass,  feet  and  seconds)? 

5.  Explain  the  significance  of  gQ  in  Newton's  Second  Law. 

6.  What  are  a  consistent  set  of  units  for  gQ . 

7.  What  is  the  distinguishing  characteristic  of  a  fluid 
(as  compared  to  a  solid)? 

8.  Describe  the  basic  difference  between  gases  and  liquids. 
VJhen  can  a  gas  be  treated  as  a  continuum? 

9.  Explain  the  difference  between  a  microscopic  and  a 
macroscopic  approach  to  the  analysis  of  fluid  behavior. 

10    Describe  a  FLUID  PROPERTY  and  give  three  or  more  examples 
of  a  fluid  property.   (Note:  We  will  use  the  term  'Fluid 
property'  to  mean  'Gross'  or  'Bulk'  fluid  property  as 
used  in  the  text . ) 

11.  Fluid  properties  may  be  categorized  as  either  INTENSIVE 
or  EXTENSIVE.  Define  what  is  meant  by  each  and  list  as 
many  examples  of  each  type  of  property  as  you  can. 

12.  What  does  the  term  "specific"  imply?  (example:  specific 
enthalpy) 

13.  What  is  the  relationship  between  density  and  specific 
volume? 

14.  Describe  an  EQUATION  OF  STATE  and  recall  one  with  which 
you  are  familiar. 

15.  State  the  PERFECT  GAS  LAW.   When  is  it  valid?   (i.e.  what 
are  it's  limitations?) 

16.  Give  a  consistent  set  of  units  for  each  term  in  the 
Perfect  Gas  Law. 

17.  Define  a  system.   Describe  the  difference  between  an 
open  system  and  a  closed  system. 
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18.  Describe  the  CONTROL  VOLUME  APPROACH  to  a  problem  and 
compare  It  with  the  CONTROL  MASS  APPROACH.   Can  you 
identify  the  INDEPENDENT  and  DEPENDENT  variables  in 
each  case? 

19.  Define  a  STATE  FUNCTION  and  a  PATH  FUNCTION.   Give 
examples  of  each. 

20.  What  is  a  PROCESS?  '  What  is  a  CYCLE? 

21.  How  is  the  Zero  th  LAW  of  THERMODYNAMICS  related  to 
temperature? 

22.  State  the  FIRST  LAW  of  THERMODYNAMICS  for  a  closed 
system. 

23.  Give  a  form  of  the  FIRST  LAW  applicable  to  an  open 
system. 

2k.       State  any  form  of  the  SECOND  LAW  of  THERMODYNAMICS. 

25..  Define  ENTROPY. 

26.  Give  equations  that  relate  entropy  to 

a. )   internal  energy, 
b . )   enthalpy. 

27.  Define  a  REVERSIBLE  PROCESS  for  a  THERMODYNAMIC  system. 
Is  any  real  process  ever  reversible?  Of  what  practical 
value  is  the  concept  of  reversible  processes? 

28.  What  are  some  IRREVERSIBLE  EFFECTS? 

29.  What  is  an  ADIABATIC  process? 

30.  If  a  process  is  BOTH  REVERSIBLE  and  ADIABATIC  it  is 
also  something  else.   What? 

31.  What  STATE  VARIABLE  is  constant  in  an  ISENTROPIC 
PROCESS?   Is  it  possible  to  have  an  isentropic  process 
that  is  not  a  reversible  adiabatic  process? 

32.  If  we  say  dS  =  0,  does  this  automatically  imply  an 
ISENTROPIC  situation?   (where  "S"  is  ENTROPY). 

33.  Define  (in  the  form  of  a  partial  derivative)  the  SPECIFIC 
HEAT'S  (Cv  and  Cp).   Are  these  expressions  valid  for 

a  material  in  any  state?   Are  they  valid  for  any  process? 

34.  For  a  PERFECT  GAS,  SPECIFIC  INTERNAL  ENERGY  is  a 
function  of  which  state  variables? 

35.  For  a  PERFECT  GAS,  SPECIFIC  ENTHALPY  is  a  function  of 
which  state  variables? 
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36.  For  the  equation:   Q  =  W  +  AE 

How  are  the  signs  conventions  defined  for  Heat  and  Work? 

37.  Process  Plots 


PI      increasing 

II 

(p  =  <£)     n=D 

^ia 

1 

II 

c 

^ 

Pvn  -a 


;  where  (£  is  a  "constant" 

n  is  any  positive  number. 


This  is  a  General  Polytropic  Process. 

a)  In  the  p-v  diagram  place  each  of  the  following  in 
the  correct  box:   0 ,  (p  ,  S. 

b)  In  the  T-s  diagram,  label  each  process  line  with 
the  correct  value  for  n  and  identify  which  fluid 
property  is  held  constant. 

c)  Are  these  plots  general,  or  are  they  restricted 
to  a  Perfect  Gas?   Why? 
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COMMENTS  OF  THE  REVIEW  QUESTIONS 


From  CRC  Tables  (17th  edition) 

TAYLOR'S  SERIES: 
f(x)-f(a)+  (x  -  a)f'(a)  +  (x  I  a)   f"(a) 


+ 


The  Taylor's  Series  expansion  enables  one  to  examine  the 
unknown  fluid  properties  at  an  arbitrary  point  (x)  in  a 
continuous  flow  field  by  an  expansion  of  the  derivatives 
of  the  properties  at  a  point  (a)  where  these  are  known. 
Note  that  all  of  the  derivatives  are  evaluated  at  point 
(a). 


dy    lim 
dx    Ax-*0 


f(x  +  Ax)  -  f(x) 


Ax 


The  important  point  here  is  that  the  derivative  is  a 
process  which  becomes  valid  in  the  limit  as  the  incre- 
mental change  in  the  state  variable  becomes  vanishingly 
small. 


3. 

SYSTEM 

(MKS) 

Absolute 
Metric (CGS' 

British 
Gravitat'nl 

Absolute 

English 

English 
Engineering 

Mass 

kilogram 

gram 

slug 

pound- 
mass 

pound- 
mass 

Length 

meter 

centimeter 

foot 

foot 

foot 

Time 

second 

second 

second 

second 

second 

Force 

newt on 

dyne 

pound- 
force 

poundal 

pound- 
force 

gC 

1  kg-m 

newton-sec^ 

1  gm-cm 
dyne-sec^ 

,  slug-ft 
lbf-sec^ 

x  lbm-ft 
poundal- 
sec2 

32.17^ 

ft-lbm 
lbf-sec^ 

*J.   A  one  pound  force  will  give  a  one  pound  mass  an  accelera- 
tion of  32. 17^  ft/sec2. 

5-   FORCE  «  MASS  times  ACCELERATION 

FORCE  =  (Constant) (Mass) (Acceleration) 
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6. 


The  CONSTANT  rectifies  the  dimensions  involved.   If  we 
assign  the  symbol  "1/g  "  to  this  constant,  Newton's  Lav; 
can  be  stated 

p  =  Ea 
Sc 

and  g  will  have  the  values  indicated  in  the  table  above 


-   1  (lbm)   32. 171*  (ft/sec2)      00    ,„,.   ft-lbm 
gc   "   1  (lbf) "   32>m   lbf-sec2 


7.  It  will  continuously  deform  under  the  influence  of  an 
applied  SHEARING  FORCE. 

8.  There  is  no  rigorous  distinction  between  gases  and  liquids 
but  generally  we  think  of  liquids  as  being  nearly  incom- 
pressible (actually  they  are  slightly  compressible)  whereas 
gases  are  very  compressible  (thus  gases  assume  the  volume 
of  a  closed  container).   The  problem  in  making  a  clear 
distinction  becomes  impossible  when  one  considers  a 
substance  which  is  above  its  critical  point. 

9.  Refer  to  Section  1.2  in  the  text  (Gas  Dynamics ,  John) 

10.   Considering  the  macroscopic  approach  of  a  small,  but  finite 
volume  of  fluid  consisting  of  a  sufficiently  large  number 
of  molecules,  we  may  define  a  FLUID  PROPERTY  by  using  the 
limit  concept  of  two  relationships  such  as: 


AF 
lim   /    normal  «, 

AA+a  ^     AA     ' 


where  "a"  is  a  small 
but  finite  area. 


In  any  case  the  resulting  characteristic  exhibited  by  the 
fluid  will  be  as  though  we  were  considering  a  single  point 
in  the  fluid.   Such  a  characteristic  (a  property)  is 
independent  of  the  path  by  which  the  point  was  reached. 

Examples  are:   pressure,  density,  temperature,  entropy, 

enthalpy,  etc. 

11.  Intensive  properties  exhibit  the  same  characteristics 
independent  of  the  amount  of  mass  present.   Temperature, 
pressure  and  specific  properties  are  INTENSIVE. 
Extensive  properties  are  dependent  on  the  amount  of 
mass  present.   Mass,  energy,  entropy  and  enthalpy  are 
examples . 

12.  The  term  "specific"  refers  to  the  amount  of  an  extensive 
property  per  unit  mass. 
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13.   Density  is  the  reciprocal  of  specific  volume: 

1 
p  =  v 

lk.      An  equation  of  state  is  a  relationship  among  properties. 
The  most  familiar  is  the  Perfect  Gas  Law: 

P  =  pRT 


where 

P 

is 

pressure 

P 

is 

density 

R 

is 

a  gas  constant 

T 

is 

temperature 

15.  The  Perfect  Gas  Lav;  is  derived  from  kinetic  theory  and 
neglects  molecular  volume  and  inter-molecular  forces. 
Therefore,  the  Perfect  Gas  Law  can  be  used  to  deal  with 
gases  under  conditions  of  low  density  (consider  the 
variables)  wherein  these  two  phenomena  are  not  significant 
factors.   Gases  at  low  pressures  and  high  temperatures 
are  closely  described  by  the  Perfect  Gas  Law. 

16.  P  —  lbf/ft^  (absolute  pressure) 
p  -  lbm/ft3 

T  —  °R  (absolute  temperature) 
ft-lbf 


R 


lbm-°R 


17.  The  term  system  is  used  in  a  broad  sense  to  identify 
clearly  whatever  it  is  that  is  under  consideration.   A 
closed  system  is  one  in  which  mass  cannot  enter  or  leave; 
an  open  system  is  one  in  which  mass  may  enter  and/or 
leave . 

18.  The  CONTROL  MASS  approach  to  a  fluid  flow  problem  utilizes 
a  closed  system,  whereas  the  CONTROL  VOLUME  approach 
considers  a  fixed  volume  in  space.   Energy  can  cross  the 
boundaries  of  either  system  whereas  mass  can  only  cross 
the  boundaries  of  an  open  system. 

19.  A  state  function  (also  called  a  point  function  or  property) 
is  an  observable  characteristic  which  describes  the 
physical  or  thermal  state  of  a  system.   Two  or  more  state 
functions  fix  the  state  of  a  system. 

Path  functions  are  transient  phenomena.   Consider  Heat 
and  Work : 

—  Systems  never  possess  heat  or  work,  but  either  or  both 

cross  the  system  boundary  when  a  system  undergoes  a 

change  of  state. 
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—  Both  heat  and  work  are  boundary  phenomena.   Both 
are  observed  only  at  the  boundaries  of  the  system, 
and  both  represent  energy  crossing  the  boundary 

of  the  system.  %^ 

—  Both  are  PATH  functions  and  inexact  differentials. 

20.  A  process  is  a  change  of  state  in  a  system  which  has  a 
unique  sequence  or  relationship  between  the  state 
variables.   A  cycle  consists  of  several  processes  which 
vary  the  state  of  a  system  and  return  it  to  the  original 
state . 

21.  The  so  called  "zeroth"  law  of  thermodynamics  defines  the 
state  variable  temperature  in  terms  of  two  systems  which 
are  in  thermal  equilibrium  with  a  third  system. 

22.  For  cyclic  operation:   EQ  =  EW   or   idQ  =  fd\J 

For  non-cyclic  operation:    Q  =  W  +  AE 

6Q  =  6V7  +  dE 

where  heat  flow  into  the  system  is  positive  and 
work  transferred  out  of  the  system  is  positive. 


Thus 

a 

new  property  has 

been 

de 

fined 

E, 

the 

total 

ener 

gy 

of 

the 

system, 

i 

dQ 
dt 

dW 
dt 

+ 

9 
3t 

/  e  P 

c.  V 

dv 

+ 

c 

/ 

.  s . 

e  p 

(V- 

dA) 
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2k .      It  is  impossible  to  make  any  transformation  whose  only 
final  result  is  the  exchange  of  a  non-zero  amount  of 
heat  with  less  than  two  heat  reservoirs  and  the  appear- 
ance of  a  positive  amount  of  work  in  the  surroundings, 
(from  Lord  Kelvin's  statement  of  the  second  law  of 
thermodynamics,  1851) 

»•     ds  .  !Sev 

26.   For  any  process: 

Tds  =  du  +  pdv 

where  "u"  is  specific  internal  energy,  "s"  is  specific 
entropy  and  "v"  is  specific  volume.   T  &  p  are 
temperature  and  pressure. 

Tds  =  dh  -  vdp 
where  "h"  is  specific  enthalpy. 


INTRODUCTION  TO  THERMODYNAMICS,  Sonntag  &  Van  Wylen. 
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27.  Conceptually,  a  reversible  process  is  one  in  which  both 
the  system  and  its  surroundings  can  be  returned  to  their 
original  conditions.   Although  NO  real  process  is 
reversible  in  a  strict  sense,  processes  which  closely 
approximate  this  ideal  may  be  realized. 

28.  An  irreversible  effect  causes  a  disturbance  in  the 
equilibrium  state  of  a  system.   To  quote  from  Elements 
of  Gas  Dynamics,  by  Liepmann  and  Roshko: 

...  A  system  is  in  equilibrium  if  it  is  free  of 
"currents".   The  term  current  refers  to  the  flux  of 
a  quantity  like  heat,  mass,  momentum,  etc.   A  current 
of  heat  flows  if  there  is  a  finite  temperature 
difference;  a  current  of  mass  flows  if  there  is  a 
difference  in  concentration  of  one  component;  a 
current  of  momentum  flows  if  there  is  a  difference 
in  velocity. 

Thus  stirring  a  fluid,  sudden  heating  or  any  process 
involving  friction  will  introduce  irreversibilities. 

29.  A  process  involving  no  transfer  of  heat  to  or  from  the 
system  is  called  adiabatic. 

30.  Isentropic. 

31.  Entropy.   It  is  possible  to  balance  the  increase  of 

entropy  within  a  system  due  to  irreversibilities  with  a 

decrease  in  entropy  due  to  heat  flow  and  achieve  a 

resultant  change  in  entropy  of  zero  (or  an  ISENTROPIC 

process).   This  process  is  obviously  neither  adiabatic 

nor  reversible. 

Recall:   dS  =  dS   +  dS. 

e     1 

32.  Yes. 

33. 

c   =  (  ^  )  c   =  (  2±) 

cp  -  v  9T  ;  lv  '  l  3T;v 

These  definitions  are  valid  for  any  process  and  for  any 
material  in  any  state. 

34.  Temperature  ONLY. 

35.  Temperature  ONLY. 

36.  Heat  into  the  system  Is  positive.   Work  done  by  the 
system  is  positive. 
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37. 


D  increasing 


The  process  lines  with  n  =  1   referring  to   T  = 
constant  ((£)   and   n  =  y      referring  to  isentropic  are 
restricted  to  Perfect  Gases.   Why? 
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UNIT  1  -  BASIC  EQUATIONS  -  OBJECTIVES 


The  student  shall  be  able  to: 

1.  State  the  basic  concepts  from  which  the  study  of  Gas 
Dynamics  proceeds. 

2.  Explain  the  difference  between  a  MATERIAL  DERIVATIVE 
(or  TOTAL  DERIVATIVE)  and  a  PARTIAL  DERIVATIVE  with 
respect  to  time. 

3.  Obtain  the  properties  at  one  point  in  a  flow  field  in 
terms  of  properties  and  their  derivatives  at  another 
point  by  means  of  a  Taylor's  Series  Expansion. 

M.   State  the  equation  used  to  relate  the  material  derivative 
of  any  extensive  property  to  the  properties  inside  of  and 
crossing  the  boundaries  of  a  control  volume.   Interpret 
each  term  in  the  equation. 

5.  Starting  with  the  basic  concepts  or  equations  which  are 
valid  for  a  CONTROL  MASS,  obtain  the  integral  forms  of 
the  ENERGY,  CONTINUITY,  and  MOMENTUM  equations  for  a 
CONTROL  VOLUME. 

6.  Simplify  the  integral  forms  of  the  ENERGY,  CONTINUITY 
and  MOMENTUM  equations  for  a  control  volume  for  the 
conditions  of  steady,  one-dimensional  flow. 

7.  Apply  the  simplified  forms  of  the  ENERGY,  CONTINUITY, 
and  MOMENTUM  equations  to  differential  control  volumes. 

8.  Show  that  by  introducing  the  concept  of  ENTROPY  and  the 
definition  of  ENTHALPY,  the  path  function  HEAT  (dQ)  may 
be  removed  from  the  ENERGY  equation  to  yield  the  expres- 
sion sometimes  called  the  'PRESSURE  ENERGY  Equation.' 

IP   +   d(vii     gdjl  +   TdS.  +  6W   =0 
P       2gc       gc         i      s 

9.  Simplify  the  PRESSURE  ENERGY  Equation  to  obtain  BERNOULLI'S 
Equation . 

10.   Apply  Newton's  Second  Law  of  Motion  to  a  differential 
control  volume  and  develop  the  MOMENTUM  equation  for 
Steady,  One-dimensional  flow. 


32 


11.   Identify  each  term  and  explain  the  restrictions  and 
assumptions  for  each  of  the  following  equations: 


1)      dp 


dA 
A 


^      =      0 


2)       <5q      =      d(Pv)       +    6vi 


3) 


dP 


d(V2) 
2gc 


shaft 

gdz      + 
g~ 


du 


+     d(V    )      +     Mi 


Tds.      + 

i 


6w. 


=       0 


6, 


P 


5) 


p 


f(V2) 
DE2Sc 

2gc 


dx      +      Sdz 
gc 


gZl 


p 


VdV 
g~ 


=      0 


XL 

2g, 


g^2 
g~ 
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Study  Guide  Unit  1 
Basic  Equations 

1.1   Equations  of  Motion  - 

a)  Review  Section  l.k    of  the  Text.   List  the  unknowns  and 
the  Equations  (by  name)  which  will  be  used  to  solve 
for  these  unknowns  in  the  study  of  compressible  fluid 
flow.   How  does  the  problem  simplify  if  the  flow  is 
considered  incompressible? 

b)  State  in  your  own  words  what  is  meant  by  (1)  the 
conservation  of  Mass,  (2)  the  Conservation  of  Energy 
and  (3)  the  Conservation  of  Momentum.   Is  Momentum 
actually  conserved? 


The  Basic  Relationships  with  which  we  will  analyze  Gas 

Dynamics  will  be: 

(1)  CONTINUITY  (Derived  from  Conservation  of  Mass) 

(2)  ENERGY     (Derived  from  Conservation  of  Energy) 

(3)  MOMENTUM    (Derived  1 

?rom  NEWTON'S  SECOND  LAW) 

(4)  An  EQUATION  OF  STATE 

P  =  $ 

p  =  pRT 

others 

(5)  THERMODYNAMIC  CONSIDERATIONS 
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1.2   Transformation  of  a  Material  Derivative  of  any  Extensive 
Property  to  a  Control  Volume  Approach. 


Let  X  be  the  total  amount  of  any  extensive  property  in 
a  given  mass. 


Let  x  be  the  amount  of  X  per  unit  mass. 


Thus:.   X  =  f  xdm  =   //Tpxdv   =  f 


pxdv 


(*1.21) 


where  dv  =  incremental  volume  element. 


1 

At  time  t  +  At  the 

1     / 

same  mass  particles 

occupy  this  region 

m  space. 

At  time  t  the  given  mass  occupies  this  region. 


dX 


Consider  what  happens  to  a  material  derivative  such  as  -tt- 


dX  _  lim 
dt  "  At->0 


(Final  value  of  X),  . . .  -  (Initial  value  of  X), 

t  +  At t 

At 


(*1.22) 


Final  X  ->  X  of  regions  2  &  3  computed  at  time  t  +  At  ~ 
Initial  X  ->  X  of  regions  1  &  2  computed  at  time  t  ~ 

LI   xdm  +  /  xdm      -   /  xdm  +  /  xdm 
2      J3  -J  t  +  At   L<l       J2  J  t 


dX 

dt    At->0 


(*1.23) 


At 
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Consider  the  two  terms 


lim 
At+0 


U 


xdm 


Jt  +  At     L^2 


£■ 


dm 


At 


But  since 


/ 


Xp  =  /  xdr 


lim 
At->0 


X2(t  +  At)  -  X2(t) 
At 


3X, 


3t 


1. 


2. 


Partial  derivative  notation  is  used  since  the  region  is 
fixed. 


As  At  -^  0,  region  (2)-*-  the  original  confines  of  the  mass 
We  shall  call  this  region  the  control  volume. 


Thus:   lim 
At+0 


U 


xdm 


Jt  +  At 


/ 


xdm 


3X        3 
*-    2  -it  c  .  v. 


At 


3t 


3t 


/ 

C  .  V 


pxdx  (*1.2* 


Consider  the  term 


/   xdm 
lim  Ly3 


t  +  At 


At+0" 


At 


Note  -  by  definition 
region  (3)  is  formed 
by  the  fluid  moving 
out  of  the  control 
volume . 


This  integral  represents  the  amount  of  X  in  region  (3) 
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© 


/ 


Let: n  be  a  unit  outward  normal 


Note:   n  ds  is  equivalent  to  dA  in  the  Text, 

Section  1.6. 


ds   be  an  increment  of  surface  dividing 
regions  (2)  &  (3) 

V.n  = component  of  V   to  ds 

(V.n)  ds  = volumetric  flow  rate 

p(V.n)  ds=mass  flow  rate 

p(V.n)  ds  At  = amount  of  mass  that  crossed  ds  in  time  At 

px(V.n)  ds  At  = amount  of  X  that  crossed  ds  in  time  At 


Thus 


out 


px(V.n)dsAt  ~  total  amount  of  X  in  region   3    (*1.25) 
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and 


lim 
At->0 


/xdm 


Jt  +  At 


At 


out 


xp(V.n)  ds 


(*1.26) 


where  S   ,  is  the  surface  where  fluid  leaves  the  control  volume 
out 


WHY  THE  LIMIT?   All  properties  are  going  to  be  evaluated  at 

the  surface  S,  thus  equation  (*1.25)  is  only 
an  approximation.   It  becomes  exact  in  the 
limit  At-*0. 

This  integral  is  called  a  flux  or  rate  of  X  flow  out  of  the 
control  volume. 


Similarly : 


with  n'  a  unit  inward  normal 


n 


lim 
At+0 


./■ 


xdm 

i_r       j  t 


At 


/ 


in 


xp(V.n» )ds 


or  this  represents  an  "X  Flux" 
into  the  control  volume. 


Note  that  this  term  appears  in  equation  (*1.23)  with  a  minus 
sign.   If  we  use  n  (outward)  here,  since  n  =  -n',  the 
combination  of  these  two  terms  can  be  obtained  by  integrating 
the  following  over  the  entire  surface: 


J, 

out 


/ 


xp(V.n)  ds   -    /  xp(V.n')  ds   = 


m 


c .  s . 


xP 


(V.n)  ds    (*1.27) 


where  c.s.  represents  the  control  surface  surrounding  the 
control  volume. 
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Final   result    becomes : 


dX 
dt 


f    «    +     f 

J      xpdv     +  / 


material    =  it  J      xpdv  +   J      «(V.n)ds    (.1.28) 

j   •   x.  •  c.v.  c.s. 

derivative 


Triple  integral    ^—Double  integral 


AGAIN:   This  applies  to  any  extensive  (Mass  dependent)  property 
whether  it  is  a  scalar  or  a  vector. 


In  Words:   The  rate  of  change  of  X  for  a  given  mass  as  it  is 
moving  around  is  equal  to  the  rate  of  change  of  X 
inside  the  control  volume  plus  the  net  flux  (flow 
rate)  of  X  out  of  the  control  volume.   (i.e.,  what 
goes  out  minus  what  comes  in), 

X   can  be  momentum  ->  x  =  V    (used  to  transform  Newton's  2nd  law) 


r2 

X      can   be    energy   ->   x   =    e   =   u   +  ^ —  +  fe 
^-  2g^         g 


V"   +  gZ 


c    °c 
(used  to  transform  1st  law  of  thermo) 

X  can  be  mass  itself  -*■  x  -  1  (used  to  transform  continuity  eq.) 

or  other.  


Example:   This  may  help  to  sort  out  what  is  and  is  not  a 

Material  Derivative  (or  the  change  of  an  extensive 
property  within  a  control  mass) 

For  a  Control  Mass: 

dQ         dW         dE 
dt     ^r   dt         dt 

\^^  t 

These  are  NOT  Material         ^-This  is_  a  Material 
Derivatives  as  they  des-        Derivative, 
cribe  the  rate  of  transfer 
of  energy  across  a  surface . 
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Thus  we  may  summarize  the  two  approaches  to  a  flow  problem: 

1.)   CONTROL  MASS  -  wherein  a  given  mass  is  observed. 
This  is  called  a  Closed  System.   (i.e.  ,  no  new 
matter  can  enter  or  leave  the  system) 

2.)   CONTROL  VOLUME  -  wherein  a  particular  fixed  volume 
in  space  is  observed.   This  is  called  an 
Open  System. 

As  an  example  let  X  (the  extensive  property)  be  Mass  (M). 

In  the  transformation  equation  relating  the  Material  Derivative 

of  an  extensive  property  (X)to  the  Control  Volume,  in 

Equation  *1.28,  X  is  mass,  hence  x  is  mass  per  unit  mass, 

or  x  =  1.   Thus  Equation  *1.28  becomes: 


dM 
dt 


3t 


III 

c  .  v. 


p  dv 


+  J  J         PV. 


c  .  s 


dA 


dM 
dt 


0 

"material 
derivative 


SO 


at 


III   P*v  \JJ 


C.  V. 


c  .  s . 


pV.dA 


This  is  the  CONTINUITY  EQUATION  for  a  control  volume.  There 
are  other  forms  of  this  equation  for  special  conditions  such 
as  steady  flow,  etc.  which  we  will  now  develop. 
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1.3   Control  Mass  versus  Control  Volume 


COMMENT  - 

-  The 

concept 

developed 

in  1.2 

is 

presei 

ated 

in 

an  abbreviate 

d  form  by 

the  Text. 

Compare 

these 

two 

parallel 

deve 

lopments 

as 

this 

important 

concept 

must 

be  thoroughly 

understo 

od. 

a)  Read  section  1.6  of  the  Basic  text. 

b)  Can  you  interpret  the  conversion  from  the  consideration 
of  a  fixed  MASS  system  and  the  properties  associated 
with  it  to  a  CONTROL  VOLUME  approach?   Which  is  more 
convenient  and  WHY? 

c)  Do  you  understand  the  significance  of  each  term  in 
Equation  (1.7)?   What  does  each  term  mean  (Explain  in 
words ) . 

d)  Show  which  term  will  drop  out  of  Equation  (1.7)  in 
any  steady  flow  situation  and  explain  why. 


Note  1:   Equation  numbers  without  an  asterisk  refer  to   * 
equations  in  the  text  (JOHN). 

Equations  with  an  asterisk  (e.g.,   *1.23)  refer 
to  equations  in  the  study  guides. 
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1.4   CONSERVATION  OF  MASS 


a.)   Read  Section  1.7  In  the  text 


COMMENT  - 

In  this 

course  we 

shall 

assume 

that  Matt 

er 

and  Energy 

are  eac 

h  separately  conserved 

,  as  we 

shall  not 

consider 

transformations 

such 

as  occur 

in  nuclear 

reactions . 



b.)   What  is  the  Principal  of  Conservation  of  Mass,  in 
words? 

c.)   Now  skip  ahead  and  read  Section  1.12.   The  real 

importance  of  the  One-Dimensional  Flow  assumption 
is  that  it  permits  us  to  easily  obtain  approximate 
solutions  (with  acceptable  engineering  accuracy) 
to  problems  whose  exact  solutions  would  be  extremely 
difficult.   Any  problem  solution  based  on  a  simpli- 
fying assumption  will  yield  results  which  vary  to 
some  extent  from  the  real-world  problem.   The  decision 
to  utilize  a  simplification  depends  on  its  correlation 
to  experimental  results,  the  availability  and 
cost  of  alternate  methods  of  solution  and  the 
accuracy  required  of  the  solution. 

d.)   Now  return  to  pg.  8  of  the  text  and  determine  in 

detail  the  restrictions  imposed  in  progressing  from 
Eq.  (1.9)  to  Eq.  (1.10) . 
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NOTE  WELL: 

The  condition 

of  ONE-DIMENSIONAL  Flow 

imposes 

NO 

restrictions 

on  density,  area,  velocity, 

(or 

any 

other  fluid  prop 

erty)  from  section 

to 

section 

in 

the  direction 

of  flow,  but  does 

say  that 

for 

a  given  location  at 

a  given  instant  of 

time , 

the 

flow  variables  are  assumed  constant  over 

the 

entire 

cross-section 

• 

e.)   Equation  (1.10)  is  often  written  pAV  =  constant, 

in  which  case  the  V  is  assumed  to  be  the  component 
of  velocity  perpendicular  to  the  area.   What  quantity 
is  the  constant? 

HINT:   Work  out  the  Units 

f . )   Is  it  correct  to  say  that  the  continuity  equation 
can  be  stated  as  "The  flow  rate  into  a  control 
volume  equals  the  flow  rate  out  of  the  control 
volume?"   Is  this  mass  flow  rate  or  volumetric  flow 
rate?   What  condition  must  exist  before  the  above 
statement  is  true?   Follow  through  examples  1.1, 
1.2.   Can  you  LIST  the  assumptions  made  for  each 
of  these  problems? 

g. )   It  is  more  convenient  for  the  case  of  multiple 
entry  and  exit  situations  such  as 


13 


to  express  the  continuity  equation  for  one-dimen- 
sional, steady  flow  as:  ]T)pAV  =  0 
h. )   An  alternate  form  of  the  continuity  equation  (for 
one-dimensional  steady  flow)  is  obtained  by  first 
taking  the  natural  logarithm  of  [pAV  =  constant] 
and  then  differentiating  the  result.   Try  it! 
You  should  get : 

dp_  +  dV  +  dA  =  0  (*1.4l) 

p     V    A 

This  form  is  useful  in  interpreting  the  changes 
that  must  occur  as  fluid  flows  through  a  duct , 
channel  or  stream  tube.   It  indicates  that  if 
mass  is  to  be  conserved,  the  changes  in  density, 
velocity  and  cross-sectional  area  must  compensate 
for  one  another.   For  example,  if  the  area  is 
constant  (dA  =  0),  then  any  increase  in  velocity 
must  be  accompanied  by  a  corresponding  decrease 
in  density.   This  form  of  the  continuity  equation 
will  also  be  useful  in  future  derivations. 


1\H 


1.5   Conservation  of  Energy 

a.)   Read  Section  1.9  of  the  text. 

Section  1.8  will  be  considered  after  the  ENERGY  con- 
cepts are  firm. 

b.)   Do  you  recall  the  essential  difference  between  a 
PATH  FUNCTION  and  a  STATE  FUNCTION?   The  key  lies 
in  the  definition  of  PATH  functions  which  are  only 
defined  in  crossing  a  system  boundary . "Heat "  and 
"pressure"  are  examples  respectively  of  a  PATH 
and  a  STATE  function.   List  as  many  of  each  as  you 
can  recall. 

PATH  FUNCTIONS  STATE  FUNCTIONS 


c.)   What  notation  is  used  to  represent  the  infini- 
tesimal change  in  a  PATH  function  by  JOHN?   How 
about  the  notation  for  an  infinitesimal  change  in 
a  STATE  FUNCTION?   Why  the  difference? 
Example:   6Q  =  6W  +  dE 

Note  that  in  equations  1.16  and  subsequent  the 
expression  dQ/dt  uses  the  same  differential  symbol 
for  Heat  (Q),  a  PATH  function  and  for  Energy  (E) 
a  STATE  function.   This  is  valid  because  we  are 
now  dealing  with  a  control  volume,  so  the  change 
in  Heat  (or  work)  is  now  a  function  of  time  only. 
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d.)   Two  important  concepts  included  in  the  development 
of  equation  (1.19)  are: 

1)  FLOW  WORK 

2)  ENTHALPY 

Flow  work  is  the  work  done  to  push  the  fluid  into 
(or  out  of)  the  system  (i.e.,  CONTROL  VOLUME). 
Since  it  is  often  difficult  to  measure  Flow  work 
being  done  on  the  system  (whereas  the  Shaft  work 
and  other  types  of  work  can  usually  be  measured), 
it  is  advantageous  to  combine  Flow  work  with  specific 
Internal  energy  (a  thermodynamic  property)  and 
thus  define  ENTHALPY  (h). 
e.  )   The  formulation  of  the  general  ENERGY  Equation  as 
expressed  in  equation  (1.19)  can  often  be  simpli- 
fied to  facilitate  obtaining  a  solution  to  practical 
problems  in  engineering.   Consider  the  case  of 
STEADY  FLOW.   We  can  summarize  the  effects  of 
STEADY  FLOW  by  two  conditions: 

1)   MASS  FLOW  IN  equals  the  MASS  FLOW  OUT  of 

the  CONTROL  VOLUME.   Why  is  this  so?   Does 

this  define  steady  flow? 


2)       a(anything)    =    0 
at 


Note  carefully  the  use 
of  partial  versus  the 
ordinary  derivative 


This  assumption  causes  the  first  term  of  eq.  1.19 
to  vanish.   Similarly,  as  in  examples  1.5  and  1.6, 
real  problems  can  be  readily  treated  by  judicious 
definition  of  the  control  volume,  which  makes  the 
surface  integral  in  eq.  (1.19)  extremely  easy  to 
evaluate . 

^6 


f . )   The  general  ENERGY  equation  may  be  viewed  as 

analogous  to  a  bank  account.   Each  of  the  terms 
represents  a  certain  type  of  transaction,  which 
must  be  balanced  against  each  other  to  balance 
the  books   in  accordance  with  the  Law  of  Conserva- 
tion of  ENERGY.   To  consolidate  your  understanding 
of  the  effect  of  each  term,  next  to  each  of  the 
statements  below,  write  the  term  from  the  GENERAL 
ENERGY  EQUATION,  Equation  (1.19),  whose  function 
is  described.   (-+- You  may  have  to  separate  certain 
terms ) 

1)   The  time  rate  of  change  of  total  energy 
within  the  control  volume 


2)   The  net  efflux  of 
KINETIC  ENERGY 


3)   The  net  efflux  of  POTENTIAL  ENERGY 
associated  with  the  earth's 
gravitational  field.  


4)   The  net  efflux  of  flow  work 
and  internal  energy  


5)   The  rate  at  which  heat  enters 
the  control  volume 


6)   The  rate  at  which  work  (other  than  flow  work) 
leaves  the  control  volume 
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g)   From  Thermodynamic  considerations  entropy  can  be  broken 
down  as  follows: 

ds  =  dse  +  ds±  (*1.51) 

where :  ^  "ds  "  represents  entropy  change  caused  by 
External  Heat  Transfer  (or  HEAT  crossing 
the  boundary).   It  is  assumed  that  this  Heat 
is  transferred  Reverslbly . 

"ds  "  represents  the  Entropy  change  caused 
by  all  the  Internal  IRREVERSIBLE  effects. 

THUS    1)  ds   =  tft      (This  can  be  +  or  -  depending 

on  the  sign  of  SQ) 

2)  ds.  is  always  positive! 

Now  recall  the  "INEQUALITY  of  CLAUSIUS."   Develop  this 
from  *1.51. 

^fiO  (*1.52) 

From  your  Review  Questions  you  can  recall 

ds£  =  0  implies  ADIABATIC 

ds.  =  0  implies  REVERSIBLE 
and   ds  =  0  means  ISENTROPIC ,  therefore  if  a  situation 
both  ADIABATIC  and  Reversible  it  is  ISENTROPIC. 

h)   To  simplify  the  application  of  the  energy  equation  (1.19) 
to  many  engineering  applications  we  will  again  apply  the 
assumption  of  ONE-DIMENSIONAL  FLOW. 
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Now  for  ONE-DIMENSIONAL  steady  flow  a  differential 
control  volume  may  be  conceptualized  as  the  dotted  area 
in  Figure  *1.51. 


Pionv) 
V 

A 
u 

P 


P+  dp  (or.-  v+dv) 

V+dV 

FLOW  Figure    *1.51 

A+  dA 

i  u+  du 
'p  +  dp 


Thus  we  may  apply  the  Energy  equation  (1.19)  to  this 
control  volume,  and  divide  out  the  flow  rate  with  the 
following  result: 


6q  =  [ (p+dp) (v+dv)  -  pv]  +  6w 


shaft 


+  [(u+du)  +  (v:dv)  + g(z+dz)]  -  [u  +  ^l  +  ££] 


2g 


2SC    Sc 


(*1.53) 


If  we  expand  this  and  neglect  terms  of  higher  than  linear 
order  we  now  have  the  Energy  Equation  in  differential  form: 


j  /   \  .  ,        .  j     d  ( V  )  .  gd  z 

6q  =  d(pv)  +  5w  ,   „,  +  du  +  ^ +  fe — 

shaft        2gn     g^, 


(*1.54) 


Fill  in  the  intermediate  steps  for  exercise. 
How  would  you  define  6q  and  6w? 
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i)   The  stage  is  now  set  to  develop  the  "PRESSURE  ENERGY 
EQUATION. "   From  the  Thermodynamic  relation 

T  ds  =  dh  -  vdp  (*1.55) 

we  can  now  substitute  Equation  *1.51  to  obtain 

dh   =   Tds^  +  Tds.  +  ^  (*1.56) 

e      1    p 

Nov;  by  substituting  from  equation  *1.56  for  dh  in  place 
of   du  +  d(Pv)   in  equation  *1.5^  and  recognizing  that 
Tds   =  5q,  we  can  obtain  the  useful  form  of  the  energy 
equation  often  called  the  PRESSURE  ENERGY  equation. 


j )   The  PRESSURE  ENERGY  EQUATION  is  the  source  of  a  familiar 
relation  from  FLUID  MECHANICS.   If  we  consider  the  flow 
as  REVERSIBLE  (ds .  =  0)  and  No  shaft  work  crosses  the 

boundary  (sw    „  =  0)  and  if  the  fluid  is  considered 

s nai  u 

incompressible  (p  =(t  )  then  we  can  integrate  the  Pressure 
Energy  Equation  between  two  points  and  obtain  BERNOULLI ' s 
Equation. 

2  2 

-2  +   1   +  — 1   =   -2  +   2   +  __2  (*1.58) 

P    2gc    gc       p    2gc    gc 

Summarize  the  assumptions  made  for  Bernoulli's  Equation 
for  your  own  review. 

1.   

2-   

3.   

4.   

5.   

so 


1.6  -  MOMENTUM 

a.  Prior  to  the  development  of  the  equation(s)  describing 
momentum,  consider  the  concept  of  momentum  —  could  you 
describe  momentum  to  a  'non-engineering  type'?   Reflect 
on  your  own  understanding  and  if  you  feel  uncertain 
about  this  concept,  break  out  your  Physics  text  and 
refresh  your  understanding  of  'MOMENTUM'  and  'IMPULSE'. 

b.  Study  Section  1.8  in  the  text,  including  examples  1.3 
and  1.4.   Are  we  considering  angular  momentum? 

How  about  acceleration  effects? 

c.  Since  neither  angular  momentum  nor  any  acceleration 
effects  are  considered,  we  find  that,  from  the  text, 
equation  (1.13)  is  in  reality  valid  only  under  quite 
restricted  conditions.   Explain  in  your  own  words  the 
restriction  related  to  each  of  the  below  listed  concepts 

a)  Inertial  Reference 

b)  No  Angular  Momentum 

c)  Control  Volume  (Fixed  or  Constant  Velocity) 

Can  you  take  equation  (1.13)  and  "cross"  both  sides 
with  the  radius  and  come  up  with  a  similar  relationship 
to  deal  with  angular  momentum? 

d.  In  Section  1.8  of  the  text,  the  material  derivative  of 
the  extensive  property  momentum  is  being  transformed  to 
a  control  volume  consideration,  which  is  then  related 
to  the  forces  acting  on  the  control  volume  by  Newton's 
Second  Law  of  Motion.   We  can  apply  to  equation  (1.13) 
in  the  text,  the  simplifying  assumption  of  Steady  and 
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One-dimensional  flow.   For  One-dimensional  flow,  the  second 
term  on  the  right  side 


m 


out 


Figure  *1.6-1 

which  describes  the  net  amount  of  momentum  per  unit  mass 
crossing  the  control  boundary  simplifies  as  follows: 


// 

c  .  s. 


t-s-x      mV 


g, 


p(?.dA)   = 


out 


mV 


(•1.6,1) 


m 


From  Section  1.5  we  can  recall  the  two  useful  effects  of 
Steady  Flow: 


at 


(anything)   =   0 


and:    ^Mass  flow  in 


2_,Mass  flow  out   =   (f  =   m 
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Therefore  Equation  (1.13)  becomes 


*— '  on 


m  ,—  _>. 

fl  .  ,  =  —  (V   ,  -  V.  ) 
fluid    g     out     in 


(•1.6,2) 


e.   Having  obtained  a  relation  which  is  now  expressed  in 

terms  of  the  mass  flow  rate  and  velocities  in  and  out  of 
the  control  volume,  we  need  now  to  examine  the  Forces 
acting  on  the  control  volume. 


FLOW         rz 


w 

Figure    *1.6-2 


In  Figure  *1.6-2  we  see  summarized  the  various  forces 

acting  on  the  control  volume,  where: 

W  =  Weight  of  Fluid  in  the  control  volume 

F  =  Upstream  pressure  forces 

Pi 

F  =  Downstream  pressure  forces 

P2 

F  =  Pressure  forces  exerted  by  the  walls 
Pw 


=  Friction  forces  exerted  by  the  wall 


w 


It  is  often  common  practice  to  group  the  last  two  together 
and  refer  to  them  as  'Enclosure  Forces'.   Where  necessary 


53 


we  can  also  include  other  forces  such  as  electrical, 
magnetic  or  other  types  of  field  forces  which  for 
simplicity  we  will  not  consider  here. 

We  establish  co-ordinates  and  designate  the  various 
forces  by  subscript  acting  on  a  differential  control 
volume  of  length  dx. 


dVA 


incoming 
properties 


P 


,  V^V 


p+dp 


W 


NOTE  WELL 


z  and  x  are  not 

orthogonal 

dimensions 


outgoing 
properties 


Figure  *1.6-3 


Since  this  is  steady,  ONE-DIMENSIONAL  flow  in  the 
x-direction,  the  sum  of  Forces  in  any  other  direction 
must  equal  zero  from  equation  (*1.6-2). 
f.   We  can  nov;  summarize  the  x-components  of  the  Forces. 
Note  we  will  establish  the  sign  convention  that  force 
components  in  the  upstream  direction  are  positive . 
Upstream  Pressure  Forces   (T) 


Fn   =  PA 
Ix 


(•1.6,3) 


54 


Downstream  Pressure  Forces  (2) 

F2x   =   -  (P  +  dP)(A  +  dA)  (*1.694) 

Expanding  and  Neglecting  higher  ordered  terms: 

F2x  =   "  tpA  +  PdA  +  Adp]  -(*1.6,5) 

Wall  Pressure  Forces  (3) 


F 


3x 


(P  +  ^)(f dL)  sin  e  (*1.6,6) 


where  ^    is  the  Mean  Wetted  Perimeter 
and    e  is  the  divergence  angle  of  the  wall 
if  we  recognize  that 

( f dL)  sin  9   =   dA  (*1.6,7) 

and  we  expand  and  neglect  higher  ordered  terms. 


F0   =   PdA  (*1.6,8) 


Friction  Forces   (V) 


If  we  define  if   as  the  mean  shear  stress  alonj 
the  wall 


F 


HX     =  -  [xf  (fdL)  cos  e  ]  (*l.6,9) 


recognizing  that   dL  cos  e   =   dx 

F„    =   -  t,  F  dx  (*1.6,10) 

4x        f  * 
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Gravity  Force   (5 


F5x   =   -  [  (A  +  4£)  dx  ]  [p  +  ^-)  ]  J*-   cos  (f>  (*1.6,11) 

where  <j>  is  the  angle  of  inclination  of  the  x  axis 
but   dx  cos  <J)  =  dz .   Expanding  and  neglecting  higher 
order  terms: 


F5x  =   "  Ap  g   dz  (n.6,12) 


g.   Having  shown  that  the  right  side  of  (*1.6,2)  is: 

—  [V   .  -  V.  ]  (*1.6,2a) 

g   L  out    mJ  » 

&c 

The  x  component  of  this  becomes: 

-^-  [V  +  dV  -  V]  (*1.6,2b) 

sc 

or 

^  dV 
sc 

You  should  equate  the  forces  to  the  right  side,  cancel 
terms  and  obtain  the  following  result : 

AdP  +  Tf  f  dx  +  PA  -£-  dz  +  p  ^-     =      0  (*1.6S13) 

sc         gc 

To  evaluate  the  Shear  term  we  will  introduce  a  non- 
dimensional  friction  co-efficient  (f)  which  will  relate 
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shear  force  to  the  dynamic  pressure.   (Sometimes  the 
friction  factor  is  defined  without  the  factor  of  4.) 


»*f 

PV2 
2Sc 

(*1.6,l4) 


We  also  introduce  the  concept  of  an  'Equivalent  diameter'. 


F 


DE  =    TT  (n.6,15) 


Then 


p  -  ■  M 

*     DE  (*1. 6,15a) 

and  substituting  these  expressions  into  Equation  (*1.6,13), 
divide  the  equation  by  density  and  we  obtain  a  useful 
form  of  the  Momentum  Equation  in  the  direction  of  fluid 
flow  for  steady,  one-dimensional  flow. 

dP      fV2   ,    ,   g  ,      VdV     n      r*i    a    iz\ 

—  +   — -* —  dx   +  -**-   dz   +   =   0      (*1.6,16) 

P      DT?2g^         g„         g^ 
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1.7  -  THERMODYNAMIC  CONSIDERATIONS  and  EQUATION  of  STATE 

a.   Read  Section  1.10.   From  your  basic  Thermodynamics  course 

you  may  recall  that  entropy  changes  can  fall  into  one  of 

two  different  categories: 


b. 


ds 


ds.  +  ds^ 
i     e 


where  ds.  is  the  change  in  entropy  due  to  irreversible 

effects  (such  as  friction) 

and  ds   (  -^   )  is  the  change  in  entropy  due  to  heat  transfer 


Comment  on  ds . 

ds-  is  always  POSITIVE 

since  all 

irreversible 

effects 

generate  entropy. 

But 

if,  in 

a  p 

rocess,  the  c 

hanges 

in  pressure  temper 

'ature 

,  etc . 

are 

vanishingly 

small , 

then  this  process 

may  b 

e  treat 

ed 

as  reversible 

'  > 

hence  ds .  =0. 

Comment  on  ds 
e 

ds  can  be  positive  or  negative  (depending  on  the 
direction  of  the  heat  transfer)  and  is  zero  if  no 
heat  transfer  occurs. 


Read  Section  1.11.   This  section  reviews  some  vital 
concepts  and  relations  which  will  be  vital  in  many 
applications.   Example  1.7  provides  an  excellent  opportu- 
nity to  follow  through  a  problem  using  a  systematic 
problem  solution  technique. 
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A  recommended  method  for  approaching  a  problem  is  provided 
below.   Although  each  of  these  factors  may  not  figure  in 
every  problem,  as  was  the  case  in  example  1.7,  they  should 
help  you  organize  a  problem  solution. 

1)  Identify  (and  sketch)  CONTROL  VOLUME,  the  dotted 
area  in  Example  1.7. 

2)  Identify  FORCES  acting  on  the  matter  inside  the 
CONTROL  VOLUME. 

3)  Identify  where  the  fluid  enters  and  leaves  the 
CONTROL  VOLUME. 

4)  Identify  any  energy  (Q  and  W)  being  transferred 
across  the  control  surface. 

5)  ESTABLISH  A  CO-ORDINATE  SYSTEM.   DON'T  FORGET  SIGN 
CONVENTIONS. 

6)  BE  ESPECIALLY  CAREFUL  WITH  THE  SIGNS  OF  VECTOR 
QUANTITIES  SUCH  AS  V  and  F. 

7)  EVALUATE  WHAT  IS  KNOWN  and  WHAT  IS  UNKNOWN.   DECIDE 
WHAT  RELATIONSHIPS  WILL  BE  USEFUL. 


The 

Basic  Relationships  with 

which  we  will  analyze 

Gas 

Dynamics  will  be 

i  • 

1) 

CONTINUITY 

(Derived 

from  Conservation  of  Mass) 

2) 

ENERGY 

(Derived 

from  Conservation  of  Energy) 

3) 

MOMENTUM 

(Derived 

from  NEWTON'S  SECOND  LAW) 

4) 

An  EQUATION 

of  STATE 

P  =  $ 
p  =  pRT 
others 

5) 

THERMODYNAMIC  CONSIDERATIONS 

Review  Section  1.13. 
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UNIT  2  -  STAGNATION  CONCEPTS   WAVE  PROPAGATION  and  MACH  NUMBERS 

-  OBJECTIVES  - 

After  successfully  completing  this  unit  you  should  be  able  to: 

1.  Explain  the  stagnation  state  concept,  how  it  is  achieved 
and  utilized. 

2.  Compare  the  concept  of  a  STAGNATION  property  with  that  of 
a  STATIC  property. 

3.  Draw  a  T  -  S  diagram  representing  a  flow  system  and  indicate 
static  and  stagnation  points  for  an  arbitrary  section. 

4.  Define  STAGNATION  ENTHALPY. 

5.  Introduce  the  stagnation  concept  into  the  "Pressure  Energy" 
equation  and  derive  the  "Stagnation  Pressure  Energy" 
equation. 

6.  Explain  how  sound  is  propagated  through  any  medium  (solid, 
liquid  or  gas) . 

7.  Define  SONIC  VELOCITY. 

8.  Utilize  a  control  volume  analysis  to  derive  the  general 
expression  for  the  velocity  of  wave  propagation  in  an 
arbitrary  medium,  starting  with  the  "Continuity"  and 
"Momentum"  equations  for  steady,  one  dimensional  flow. 

9.  Recall  the  relations  for: 

a.  Speed  of  Sound  in  an  arbitrary  medium. 

b.  Speed  of  Sound  in  a  Perfect  Gas. 

c.  Mach  Number. 

10.  Discuss  the  propagation  of  signal  waves  from  a  moving  body 
in  a  fluid,  and  explain  what  is  meant  by  'Zone  of  Action', 
'Zone  of  Silence',  'Mach  Cone',  and  'Mach  Angle'. 

11.  Graphically  portray  the  simplified  flow  over  a  wedge 
shaped  body  for  SUBSONIC  and  SUPERSONIC  flow. 

12.  Express  the  basic  equations  in  terms  of  Mach  number  for 
a  Perfect  Gas  by  showing  that : 

a.   the  "Continuity"  equation  can  be  written  as 

/ygc 

m      =      PAM   y-jfir 
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b.   the  "Energy"  equation  can  be  written  as 


(1  +  ^^  M2)dh   +  '*—  hdM2   =   6q  -  6w 


c.   the  "Momentum"  equation  can  be  written  as 


IP   +   X  dM2   +   Y  M2  dT   +  I  M  fdx   =   0 

P      2  u        2     T      2    D„ 

hi 
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UNIT  2  -  STUDY  GUIDE 


2 . 1   The  Stagnation  Concept 

a.  We  introduce  the  concept  of  a  reference  state  defined  as 
that  thermodynamic  state  which  would  exist  if  the  fluid 
were  brought  to  zero  velocity  and  zero  potential.   To 
yield  a  consistent  reference  state  we  must  qualify  how 
this  "stagnation  process"  would  be  accomplished.   The 
stagnation  state  must  be  reached: 

(a)   without  any  energy  exchange   (Q  =  V,T  =  0) 
•  (b)   without  losses. 
By  virtue  of  (a),  dSe  =  0,  and  from  (b),  dS-j_  =  0.   Thus 
the  stagnation  process  is  isentropic. 

b.  We  can  imagine  the  following  example  of  actually  carrying 
out  the  stagnation  process.   Consider  fluid  which  is 
flowing  under  condition   1   in  Figure  *2.1  (these  condi- 
tions are  referred  to  as  the  "static"  conditions).   At 
(2)  the  fluid  has  been  brought  to  zero  velocity  and  zero 
potential  under  the  above  restrictions. 


I?  ,  X  ,  T,  *  0 


STREAMTUBE 


Figure    *2.1 


ARBITRARY       z  =0 


c.   Apply  the  "energy"  equation  to  the  streamtube  and  you  have: 

0  ,0     0       0 


(*2.1,1) 


which  simplifies  to 
Vl2    gzl 

h-,   +   7T±-  + 

1    2g^    g^ 


h, 


(•2.1,2) 


But  condition  (2)   represents  the  "stagnation  state" 
corresponding  to  the  "static  state"  Q.   We  thus  call 
hp  the  stagnation  or  total  enthalpy  corresponding  to 
state  (T)  and  designate  it  as  h^  . 
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Thus 


=   h. 


2g, 


gz. 


(*2.1,3) 


Or,  for  any  static  state  we  have  in  general: 


=   h 


V' 


2g, 


(•2.1,5) 


The  introduction  of  the  stagnation  (or  total)  enthalpy 
makes  it  possible  to  write  equations  in  a  simplified 
form.   For  example,  the  steady  flow  energy  equation 
becomes : 


dq   =   dw   +  dh, 
1       s     t 


(*2.1,6) 


or 


\   +  ql-2 


=   h, 


+  w 


(*2.1,7) 


'1-2 


You  should  note  that  the  stagnation  state  is  a  reference 
state  and  may  or  may  not  actually  exist  in  the  flow 
system.   Also,  in  general,  each  point  in  a  flow  system 
has  a  different  stagnation  state  as  shown  below: 


T 


® 


Stagnation 
States 


States 


Physical  System 


© 


Thermodynamic 
State  Plot 


Figure  *2.2 
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Also  one  must  realize  that  when  the  frame  of  reference 
is  changed  then  stagnation  conditions  change,  although 
the  static  conditions  remain  the  same.   (Static  properties 
are  defined  as  those  that  would  be  measured  if  the 
measuring  devices  move  with  the  fluid.) 

Consider  still  air  and  the  earth  as  a  reference  frame. 


V  =  0 


P  =  14. 7  psia 
Pt  =  14. 7  psia 


T  =  520°R 
Tt  =  520°R 


Figure  *2 . 3 


^   \\\\\\\  \  \\\\\\  \ 


In  this  case,  since  the  velocity  is  zero  (with  respect  to 
the  frame  of  reference),  the  static  and  stagnation 
conditions  are  the  same. 

Now  let's  change  the  frame  of  reference  by  flying  through 
the  air  on  a  missile  at  600  ft/sec.   As  we  look  forward  it 
appears  that  the  air  is  coming  at  us  at  600  ft /sec .   Its 
static  pressure  and  temperature  are  unchanged  at  14.7  psia 
and  520°  R  respectively. 


V  =  600  ft/sec 

P  =  14.7  psia  1   P1 

T  =  520°  R  ?   T, 


__  ?  Figure  *2  .  4 


g. 


But  now  the  air  has  a  velocity  (with  respect  to  the  frame 
of  reference)  and  thus  the  static  and  stagnation  condi- 
tions are  different . 

You  will  soon  learn  how  to  compute  the  stagnation  condi- 
tions.  Incidentally,  is  there  any  place  in  this  last 
system  where  the  stagnation  conditions  actually  exist? 
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2 . 2   The  Stagnation  Pressure  Energy  Equation 


a.   Consider  two  stations  in  a  flow  system  that  are  very 

close  together  and  with  thermodynamic  states  differentially 
separated  as  shown  in  Figure  *2.5«   Also  shown  are  the 
corresponding  stagnation  states  for  these  two  stations. 


T 


o  o 


Figure  *2.5 


-o 


TdT 


2 


We  may  write  the  following  property  relation  between 
points  1  and  2: 


Tds   =   dh  -  vdp 


(•2.2,1) 


Since  this  is  a  valid  property  relation  we  may  also  write 
this  between  points  1,  and  2  . 


Ttdst 


dht  "  Vtdpt 


(*2.2,2) 


However,     ds,  =  ds 


and 


ds   =  ds   +  ds. 


Thus  we  may  write: 


Tt(dse  +  ds.)   =   dht  -  vtdpt 


(•2.2,3) 


Recall  the  "energy"  equation  written  in  the  form: 


5q   =   5w   +  dht 


(*2.2,4) 
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Substitute   equation   *2.2,3   into   equation    *2.2,4    and   you 
should   obtain: 


3q      =      5v7s    +   Tt(dse    +   dsi)    +   vtdpt 


(•2.2,5) 


Also  recall  that:   8q  =  Tds 


(•2.2,6) 


By  substituting  equation  *2.2,6  into  *2.2,5  and  noting 
that  - 

1 
t     pt 

you  should  obtain  the  following  equation  which  is  called 
the  "Stagnation  Pressure  Energy"  equation: 


dP 


t  A 


8w   +  ds  (T,-T)  +  T.ds.   =   0 
s     e   t       t   i 


(*2.2,7) 


b.   Consider  what  happens  when 

(a)   There  is  no  work  transfer 


8ws  =  0 


(b)   There  is  no  heat  transfer 


dsQ  =  0 
e 


(c)   There  are  no  losses 

dP 
Under  these  conditions   — 


t 


Pt 


=  0 


ds±  =  0 


or   dPt  =  0 


or    P,  =  constant 


(«2.2,8) 
(•2.2,9) 


Note  that  in  general  the  total  pressure  will  not  remain 
constant;  only  under  a  special  set  of  circumstances  will 
equation  *2.2,9  hold  true. 
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2.3   Velocity  of  Sound 


Read  Sections  2.1,  2.2  and  2 . 3  in  the  text. 

Recall  from  your  previous  work  in  fluid  mechanics,  the 
definition  of  a  stream  line.   What,  if  any,  properties 
reamin  constant  along  a  stream  line?   What  significance 
does  the  distance  between  stream  lines  have  in  an 
incompressible  medium? 

Consider  carefully  the  implications  of  section  2.2  with 
regard  to  the  velocity  of  sound  in  a  compressible  medium. 
It  is  important  to  note  that  a  sound  wave  results  from  an 
inf initesimally  small  impulse.   The  text  uses  a  moving 
control  volume  which  encloses  the  wave  to  analyze  this 
unsteady  situation  as  a  steady  flow  problem.   We  can  view 
this  process  as  a  simple  change  of  reference  from  that  of 
a  fixed  observer  viewing  the  moving  wave  to  that  of  an 
observer  who  is  "riding"  the  wave.   The  moving  observer 
experiences  the  fluid  flowing  from  right  to  left  in 
Figure  2.5  at  a  velocity  of  sound  (a)  upstream  and 
velocity  (V-a)  downstream.   This  does  not  change  the 
static  conditions.   What  about  the  stagnation  conditions? 

Follow  through  Example  2.1.   The  author  of  our  text 
introduces  the  relationship  of  'compressibilities'  which 
may  be  unfamiliar  to  you.   What  other  ratio  is  expressed 
by  (y)?   We  need  not  concern  ourselves  at  this  stage  with 
situations  in  which  air  does  not  behave  approximately  as 
a  perfect  gas,  other  than  to  examine  each  problem  with 
this  in  mind.   In  general,  for  flow  situations  where 
extremely  high  pressures  and/or  extremely  low  temperatures 
are  encountered  the  'Perfect  Gas'  approximation  is  not 
valid.   Otherwise  we  may  treat  air  and  other  gasses  (if 
no  appreciable  molecular  interaction  is  taking  place)  as 
a  perfect  gas. 

Read  sections  2.4  and  2.5  in  the  text.   In  Figure  2.10  in 
the  text,  the  circles  indicate  the  location  of  the  original 
sound  wave  emitted  while  at  point  '0'.   Figure  *2.6  in  the 
Study  Guide  shows  the  location  of  the  sound  wave  initially 
emitted  at  times,  t  =  0,  1,  &  2  as  they  appear  at  time, 


t  =  3. 


emitted  at  t=0 
emitted  at  t=l 


1 
V  =  p-  a 


Figure  *2.6 


emitted  at  t=2 
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f.  Notice  the  distance  between  successive  sound  waves  is 
less  in  the  direction  of  travel  of  the  point  projectile. 
Now  compare  Figure  *2.6  to  Figure  2.11  in  the  text. 

In  Figure  2.11  the  point  is  overtaking  the  sound  waves 
it  emits  as  it  moves.   The  waves  are  left  behind  creating 
a  situation  in  a  sense  like  the  wave  created  by  dragging 
the  point  of  a  stick  rapidly  across  the  surface  of  a 
still  body  of  water.   The  Mach  Wave  (or  Mach  Cone) 
produced  by  a  point  source  is  an  infinitely  weak  pressure 
wave,  but  the  addition  of  a  multitude  of  such  waves 
produced  by  a  finite  body  results  in  the  shock  shown  in 
Figure  2.13  in  the  text. 

g.  It  should  be  emphasized  again  that  both  velocity  and 
speed  of  sound  are  local  conditions. 


air 


Now  consider  the  problem  below 


Q) 


Note:   No  information 

may  be  assumed  about- 

pressure,  area, 

temperature,  etc.- 

changes  from  1  to  2. 

V1    =  500  ft/sec     V2  =  1000  ft/sec 


M1  =  .5 


M    =   9 

n2 


Therefore : 


a1   =  1000  ft/sec   a2  =  ? 


Is  enough  information  given  to  compute  a2 
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2.k    -   Basic  Equations  in  Terms  of  Mach  Number 

a.  The  importance  of  Mach  Number  as  a  parameter  which 
controls  the  behavior  of  fluid  flow  suggests  the  value 
of  expressing  the  basic  equations  developed  in  Unit  1 
in  terms  of  Mach  Number  as  well  as  velocity.   This  can 
be  done  for  steady,  one  dimensional  flow  quite  readily, 
for  the  case  of  perfect  gases. 

b.  For  the  CONTINUITY  equation  — 

Equation  (1.10)  from  the  text  may  be  expressed  for  steady 
flow  as : 

m  =   PAV  =  <£  («2.4,1) 

where  (m)  is  mass  flow  rate. 

From  the  Perfect  Gas  Lav;: 


P   =  ~  (*2.4,2) 

and  the  definition  of  Mach  Number: 

V  =   Ma  (*2.4,3) 

Now  recall  the  expression  for  sonic  velocity  in  a  Perfect 
Gas: 

a   =  V7i~RT  (*2.M) 

where  y,  g   and  R  are  defined  in  Section  1.11  of 

the  text. 

We  thus  have: 


m   =   PAV   =   ^i^/y^RY   =   PAM\-^        (*2.4,5) 

Thus,  for  the  steady,  one  dimensional  flow  of  a  Perfect 
Gas,  the  "Continuity"  equation  becomes: 


PAMy^   =  <£  (*2.4,6) 
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c.   Differentiating  the  "Energy"  equation,  Equation  (1.19)  in 
the  text,  and  neglecting  the  potential  term,  we  can  show 
for  steady  flow  that: 


V2 
d[h  +  ^-]   =   6q  -  <5w  (*2.4,8) 

6c 

You  saw  that  for  a  perfect  gas: 

a2   =   YgcRT  (*2.4,9) 

hence 

2 
d[h  +  M  (^RT)]   =   6q  -  5ws  (*2.4,10) 

We  can  re-write  this  expression  as: 

d[h  +  (  Xzl  )   ^lljW]   =   5q  _  5Wg  (*2.H,11) 

From  section  1.11  recall  that  for  a  perfect  gas: 

c 
(_I_)   =   _£  (*2  4  12) 

and 

dh   =   c  dT  (*2.*»,13) 

p  ' 

Show  that : 

d[h(l  +  *~  M2)]   =   6q  -  6wc  (*2.4,lH) 

Now  use  the  product  rule  for  differentials  to  show  that : 


(1  +  ^~  M2)  dh   +   *~  hdM2   =   «Sq  -  6w£       (*2.4,15) 


70 


d.   We  derived  In  Unit  1.6  of  the  Study  Guide  an  expression 
for  the  "Momentum"  equation  in  steady,  one  dimensional 
flow: 

dP  fV  g^z  VdV 

0£     +  ^1 —  dx      +      a_     +      El     =      o  (*2.4,16) 

p  DE2gc  gc  gc 

Divide  through  by  P,  multiply  each  term  by  p,  neglect  the 
dz  term,  and  rewrite  the  last  term  to  yield: 

d£  +  pfXl_  +  P^f_  =     o  C»2.4,17) 

P      PDE2gc      P2gc  J 

Now,  substitute  pRT  for  P,  multiply  the  second  and  third 
terms  by  y/y  and  regroup  the  second  term  to  obtain: 

I?   +   X  M2  JL  dx   +   ^dV    =   o  (*2.4  18) 

P      2  '   D™   x      2yg  RT  l  ^  ,1Cj 

E  °C 

The  third  term  can  be  modified  by  considering  the  product 
rules  for  the  differential  of  the  two  variables  M  and  T: 

I^  =   X^(a!«ii  .   I  etl2  +  I  H2  f        V2.H.19) 
a  a 


So  now  the  "Momentum"  equation  may  be  written: 


dP   +   I  dM2   +   IM2M   +  1  M2(f  dx}   =   Q   (*2>4a20) 

E 
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UNIT  3  -  ISENTROPIC  FLOW  -  OBJECTIVES 


The  student  shall  be  able  to: 

1.  Define  an  ISENTROPIC  PROCESS  and  explain  the  relation- 
ship among  reversible,  adiabatic  and  isentropic  processes 

2.  Show  graphically  how  pressure,  density  and  area  vary  in 
steady,  one-dimensional,  isentropic  flow  as  Mach  Number 
ranges  from  zero  to  supersonic  values. 

3.  Explain  the  difference  between  STATIC  and  STAGNATION 
properties.   Given  all  the  static  properties,  compute 
the  stagnation  properties. 

M .   Describe  what  is  meant  by  a  "choked"  flov;  passage. 

5.  Compare  the  function  of  a  NOZZLE  and  a  DIFFUSER.   Sketch 
physical  devices  that  perform  as  each  for  subsonic 

and  supersonic  flow. 

6.  Simplify  the  basic  equations  for  CONTINUITY,  ENERGY 
and  MOMENTUM  to  relate  differential  changes  in  density, 
pressure  and  volume  to  a  differential  change  in  area 
for  steady,  one-dimensional  flow  through  a  varying  area 
passage . 

7.  Define  the  reference  condition,  (*)  and  the  properties 
associated  with  it.   (i.e.,  A*,  P*,  T* ,  p*,  etc.) 

8.  Define  stagnation  temperature  (T^ )  of  a  perfect  gas  in 
terms  of  temperature  (T) ,  Mach  Number  (M)  and  ratio  of 
specific  heats  (y). 

9.  Define  stagnation  pressure  (Pt )  of  a  perfect  gas  in 
terms  of  temperature  (T),  Mach  Number  (M)  and  ratio  of 
specific  heats  (y)- 

10.  Define  stagnation  enthalpy  (ht )  in  terms  of  enthalpy  (h) 
and  velocity  (V) .  Put  in  terms  of  M  and  y  for  a  perfect 
gas. 

11.  Depict  graphically  the  general  relationships  of  velocity, 
density  and  area  in  varying  area  adiabatic  flow  where 
the  Mach  Number  varies  from  zero  to  supersonic  values. 

12.  Express  the  loss  relationship  in  adiabatic  flow  (Asj_)  as 
a  function  of  stagnation  pressures  (P-^)  or  reference 
areas  (A*)  between  two  points  in  the  flow. 
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13.  Cite  the  necessary  conditions  for  the  reference  area 

A*  to  remain  constant  for  all  points  in  the  flow  field. 

14.  Derive  the  working  equations  for  a  perfect  gas  relating 
property  ratios  between  two  points,  in  both  adiabatic 
and  isentropic  flows,  as  a  function  of  Mach  Number  (M) , 
ratio  of  specific  heats  (y)  and  change  in  entropy  (As). 

15.  State  and  interpret  the  relation  between  stagnation 
pressure  (Pt)  and  the  reference  area  (A*)  for  the 
process  between  two  points  in  adiabatic  flow. 

16.  Utilize  the  adiabatic  and  isentropic  flow  relations 

and  the  isentropic  tables  to  solve  typical  flow  problems 
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STUDY  GUIDE  -  UNIT  3 
ADIABATIC  AND  ISENTROPIC  FLOW 


N 


Engineering  problems  typically  involve  many  factors  which 
affect  the  fluid  flow  conditions  and  frequently  one  of 
these  factors  is  predominant.   In  this  case  it  is 
possible  to  consider  only  the  controlling  factor  and 
develop  a  simple  solution  to  the  problem.   These 
approximate  solutions  can  be  applied  with  great  confidence 
to  many  engineering  situations.   To  the  engineering 
student,  this  approach  also  provides  insight  concerning 
the  significance,  magnitude  and  relative  importance 
of  each  of  the  factors  considered,  as  the  basic  state 
variables  are  changed. 


3.1  -  Variable  Area  Flow  in  General 

a.  We  will  now  consider  flow  in  a  variable  area  channel 
without  heat  transfer,  or  "adiabatic"  flov;.   Many  real 
life  situations  fit  into  this  category.   We  have 
learned  that  for  many  situations  air  and  other  gases 
flowing  in  a  continuous  field  obey  the  Perfect  Gas  Law 
reasonably  well.   Since  this  is  the  case,  we  will 
develop  working  equations  using  the  'perfect  gas' 
assumption.   Additionally,  we  are  motivated,  for  the 
sake  of  simplicity,  to  continue  to  assume  steady,  one- 
dimensional  flow  with  no  work  (that  is,  shaft  work) 
done  on  or  by  the  system. 

If  we  consider  the  adiabatic  flow  as  it  occurs  in  real 
situations,  some  accounting  must  be  made  for  losses 
which  occur  as  a  result  of  friction  and  other  irreversi- 
ble factors.   When  these  factors  are  considered  negli- 
gible and  the  flov;  is  taken  to  be  reversible,  the  flov; 
meets  the  criteria  for  isentropic  flov;.   This  concept 
was  covered  in  brief  in  the  Review  Questions.   We  may 
view  the  isentropic  flov;  case  as  the  'ideal'  or 
'Standard'  with  which  to  measure  all  adiabatic  flov;. 

b.  Read  sections  3.1  and  3 . 2  in  the  text. 

c.  What  common  simplification  is  used  to  obtain  equations 
3.1  and  3.2  from  the  basic  equations?   When  is  it 
justifiable  to  neglect  the  higher  ordered  terms  and 
assume  that  these  versions  of  the  basic  equations 
(3.1,  3.2,  etc.)  are  correct? 
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d.   To  gain  some  insight  into  the  manner  in  which  the  flow 

variables  are  changing  when  a  fluid  encounters  a  variable 
area,  it  will  be  convenient  to  restrict  the  analysis  to 
flow  with  no  losses.   We  shall  start  with  the  energy 
equation: 

S.q   =   6wB  +  dht  (*3.1,D 

Since  we  have  made  the  assumptions  of  no  heat  transfer 
or  shaft  work: 

dht  =  0  (•3-1,2) 

Recalling  the  definition  of  stagnation  enthalpy  (neglecting 
potential  energy  change): 

ht   -  h  +  £  ■   (*3.1,3) 

Differentiating: 

,,     ,.  ,  VdV 

dh,  =  dh  +  

&c 
but: 


dht  =  0 


so  : 


dh   =   - 


XIX  (*3.1,i*) 


Now,  we  may  start  with  equation  (1.22),  from  the  Second 
Law  of  Thermodynamics: 

Tds   =   dh  -  ^  (*3.1,5) 

p 


Since  this  process  is  both  reversible  and  adiabatic,  it 
is  isentropic .   Thus  ds  =  0. 
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So  finally: 


dh   = 


dP 


(•3.1,6) 


Considering  the  equations  (*3.1,^)  and  (*3.1,6)  we  obtain 


VdV 


dP 
p 


or 


dV  =   - 


gcdP 
pV 


(•3.1,7) 


We  introduce  this  into  equation  3-1  and  the  differential 
form  of  the  continuity  equation  becomes: 


dp_ 
P 


dA 
A 


g  dP 

;v2 


=  o 


(*3.1,8) 


which  can  be  solved  for 


dP   _   Vf_   r  dp_   dA  -, 
P      go     P     A  J 


The  definition  of  sonic  velocity 


(*3.1,9) 


•  11 

'C  5p 

s=constant 
can  be  written  for  our  case  as: 


dP 
'c  dp 


(2.3) 


because  this  flow  is  isentropic 


Now,  rearranging  equation  2.3 

a2 

dP   =   —  dp 


(*3.1,10) 
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Substituting  this  expression  for  dP  into  equation  *3.1,9 
yields : 

T  "  4f  V  +  f]  (•3.1.11) 

a 


From  the  definition  of  Mach  Number 

2 
a 

show  that : 

¥■  -  C  -^y  ]  X  (*3.1,12) 


If  we  now  substitute  equation  *3.1,12  into  equation  3.1 
in  the  the  text  you  can  show  that: 


*v  =  -  [   1   ]  <£  (-3.1,13) 

v  1-JYT    A 

But  we  recall  that: 

-g 

dV  =  — =£.  dP  (*3.1,7) 

p  v 

or 

ffi  =  -  -%  dP  (-3.1,7a) 

JTT 

Substitute  this  expression  for  y  into  equation  *3-l»13 
and  solve  for  dP: 


dP   ■   #  c-^f  (3-5) 

Bc     1-M 
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We  now  collect  the  principal  relations  which  will  be 
referred  to  in  the  following  sections. 


sc    1-NT    A 


(3.5) 


(•3.1,12) 


(•3.1,13) 


e.   Consider  now  what  must  be  happening  when  fluid  flows 
through  the  variable  area  channel.   Fill  in  the  blanks 
below.   Assume  that  the  pressure  is  always  decreasing. 
Thus  dP  is  negative. 


From  equation  3.5  you  see  that  if  M  <  1,  dA  must  be  

indicating  that  the  area  is  ,  whereas  if  M>  1, 

dA  must  be  and  the  area  is 


From  here  you  can  move  to  equation  *3-l,12.   If  M <  1  and 

dA  then  dp  must  be  .   If  M  >  1  and 

dA  then  dp  must  be . 


Looking  at  equation  *3-l313  reveals  that  if  M<  1  and 

dA  then  dV  must  be  meaning  that  velocity 

then  dV 


is 


must  be 


,  whereas  if  M  >  1  and  dA 
and  velocity  is 
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f.   We  can  summarize  the  effects  and  variation  of  velocity, 
density  and  area  as  Mach  Number  varies  from  zero  to 
supersonic  values.   Refer  to  Figure  *3«1  and  recall  the 
continuity  equation: 


dp_   dA   dV  =   n 
p     A    V 


(*1.h1) 


P,A,V 


Mrl.O    M>  1 


Summary  of  variation  of  p ,  V  and  A. 


M  <<  1 

M  =  1 

M  >  1 

1.  Density^  constant 

2.  Area 

compensates 
Velocity 

1.  Density 

compensates 
Velocity 

2.  dA  =  0 

1.  Density 

compensates 
Area 

and 
Velocity 

g.   Now  read  sections  3-3  and  3-*'  in  the  text.   Consider 
carefully  the  results  pictured  in  Figures  3.3,  3.^  and 
3.5.   The  situation  depicted  in  Figure  3. 5  in  the  text 
can  only  exist  if  the  proper  pressure  conditions  exist 
throughout  the  nozzle. 
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h.   A  'nozzle'  is  a  device  which  converts  enthalpy  (or  pressure 
energy,  for  the  case  of  incompressible  fluids)  into 
kinetic  energy.   A  'diffuser'  is  a  device  which  converts 
kinetic  energy  into  enthalpy.   This  distinction  is  often 
overlooked  and  a  device  is  labeled  a  "nozzle"  simply 
because  it  is  a  convergent  section,  with  no  regard  for 
the  flow  conditions. 


3.2  -  Computation  of  Stagnation  Conditions  for  a  Perfect  Gas 

a.  Read  the  first  part  of  Section  3.^  (from  page  39  to  the 
middle  of  page  4l). 

b.  Particularly  note  equations  3-6  and  3- 7.   These  will  be 
used  frequently  and  should  be  learned. 


3.3  -  Adiabatic  Plow  for  a  Perfect  Gas 


a.  We  now  return  to  the  problem  of  adiabatic  varying  area 
flow  and  include  any  losses  that  might  exist.   To  make 
this  problem  tractable  we  will  assume  that  the  fluid 
is  a  perfect  gas. 

b.  From  earlier  work  recall  the  following  relations  which 
hold  true  for  a  perfect  gas : 


Ah   =   c  AT  (1.25) 

c   -  c    =   R  (1.26) 

P    v 


c 
v  f     Prom  section 


1.11  in  text 


c 


P       Y-l 


±r   R 


From  3 . 1  we  showed  that 

&q   =   SwQ  +  dh,  (*3.1,D 


s    ~"t 
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So  for  adiabatic  flow  with  no  work 


dht   =   0 


(•3.1,2) 


or 


ht  -  <t 


This  is  true  for  any  fluid  under  conditions  where  Q  =  W  =  0 
From  equation  1.25,  if  one  considers  a  perfect  gas, 
the  stagnation  temperature  (T^.)  is  also  constant. 


c.   Review  the  relations  we  already  know  to  determine  which 
ones  may  be  useful  in  analyzing  adiabatic  flow.   First 
consider  the  continuity  equation  for  steady,  one- 
dimensional  flow: 


m 


=   pAV   =  (J 


(*3.3,D 


Examine  the  flow  between  two  arbitrary  points  in  a  flow 
system  such  as  in  Figure  *3-2. 


FLOW 


v         \ 

V                \                 V                   V                \ 

/ 

Figure  *3-2 


Equation  *3»3,1  may  be  written 


PlAlVl   =   p2A2V2 


(•3.3,2) 


or 


A- 


PlVl 

p2v2 


(•3.3,3) 
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From  the  definition  of  Mach  Number: 

V  =  Ma  (*3.3,*D 


and  the  expression  for  the  speed  of  sound  in  a  perfect 
gas: 


a   =   VYgcRT  (*3.3,5) 

and  recalling  the  equation  of  state  for  a  perfect  gas: 

P   =   PRT  (*3.3,6) 

we  can  re-write  equation  *3.3>3  as: 


A2     P1M1    T2   1/2 

%  =  ^?2   ^  ^l  ^  (*3.3,7) 


Work  out  the  intermediate  steps  from  *  3 • 3  >  3  to  *3«3>7  in 
order  to  corroborate  this  expression. 


d.   We  will  now  turn  to  the  energy  equation.   We  have  seen 
that: 

ht   =  h  +  X-  (*3.3,8) 

Hence  the  energy  equation  may  be  written  between  two 
points  as: 


K        +  q-,_>P   =  w      +  h  (*3.3,9) 

t1      l->2      s1_>2      t2 


Since   T,   =  T,    in  a  perfect  gas  under  adiabatic  flow 
t1  t2 

conditions  with  no  shaft  work,  you  can  now  re-write  equation 
1.25  from  the  text,  between  any  arbitrary  point  in  the 
system  and  stagnation  conditions.   You  should  get: 
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ht  -  h 


c   [T,  -  T] 
p    t 


(*3.3,10) 


Now  using  equation  *3-3,8  we  can  show  that 


Tt   =   T  [1  +  *|i  M2] 


(3.6) 


The  development  of  this  equation  is  found  in  section 
3.^  in  the  text. 

Since  Tt  is  a  constant  between  any  two  points  in  the  flow, 
we  can  now  write: 


-1  „  2 


-1  „  2 


Tx  [1  +  l=i  M^]   =   T2  [1  +  ^  M2^] 


(•3.3,11) 


or 


T. 


1  +  l=i  M^ 


1  +  I=i  M22 


(•3.3,12) 


e.   Prom  the  stagnation  energy  equation  for  adiabatic  flow 
of  a  perfect  gas  with  no  work  : 


dP, 


ds 


R 


(•3.3,13) 


If  we  integrate  both  sides  between  two  points  (1  &  2), 
we  get : 


In 


t2 

'tl 


-  (s  .   -  s  . 

X2     Xl 
R 


(•3. 3,1*0 


Since   dse  =  0,  Asj_  =  As   so,  taking  the  anti-log  of 
both  sides  of  equation  *3.3,1^  yields: 


(•3.3,15) 


This  expression  is  of  particular  value 
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f.   Stagnation  pressure  is  defined  as  the  pressure  obtained 
if  the  flow  is  brought  to  rest  isentropically . 

Between  two  points  we  can  recall  from  section  1  (for  an 
isentropic  process): 


(  Tl 


(1.30) 


Let  condition  2  be  the  stagnation  state  and  condition  1 
be  the  static  state.   Use  equations  1.30  and  3.6  to  show 
that: 


X 


-£«(!  +  X^mV"1 


(3.7) 


Thus  between  two  points: 


t2 


tl 


P2(l  +  izi  M22)Y-1 

P1(l  +  3Czl  M12)Y^I 


=   e 


As. 
R 


(•3.3,16) 


Rearrange  this  equation  to  show  that 


1    +   ^~  M22 

Y-l 

1    +   X— ^  M 

1           2         1     J 

AS 
'  R 


(*3.3,17) 


We  can  now  utilize  equation  *3.3>17  to  obtain  a  relation 
which  expresses  the  change  in  area  in  terms  of  Mach 
Number  and  a  loss  term. 

From  continuity  show  that : 


A, 
T. 


¥i 

P2M2 


1/2 


(*3.3,18) 


8H 


You  can  now  substitute  from  equations  *3.3,17  and  *3.3,12 
to  obtain: 


A. 

T. 


Y-l      2  -i 


+  X7^  ] 


1   J 


Y-l  /  As 


1/2 


R  \l  _1 


i  +  1~  i^2 

Ll  +  ^M22  j 


(•3.3,19) 


By  combining  terms,  equation  *3-3,19  can  be  written: 


A. 


M. 

K 


1  +  — —  M 
X     2    2 

L  -1     2   Hl 


(•3.3,20) 


3.4  -  Isentropic  Tables 


Read  the  remainder  of  section  3.^  in  the  text.   Note  the 
introduction  of  a  new  reference  condition  denoted  by  *. 
At  this  location  the  Mach  Number  is  unity  and  this 
condition  can  be  reached  by  many  processes.   This  is  a 
hypothetical  reference  condition  which  may  or  may  not 
exist  in  the  system  (as  was  the  case  for  the  stagnation 
condition) . 

One  can  now  simplify  many  general  equations  that  have  been 
developed  by  assuming  no  losses. 

Consider  an  isentropic  flow  case  viewed  in  the  T —  S  plane 


T 


i 

r 
i 

i 

*  i 


stagnation  point 


an  existing  subsonic 
flow  condition 

M  =  1.0 


Figure  *3-3 
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By  utilizing  the  expressions  for  area,  pressure  and 
temperature  ratios  and  letting  one  state  point  be  a 
reference  condition,  if  we  know  what  the  fluid  is  (i.e., 
the  value  of  y)  then  we  can  construct  tables  which  relate 
these  ratios  to  various  values  of  Mach  Number  for  the 
isentropic  case. 

As  an  example,  let  2  be  an  arbitrary  point  in  the  flow 
system  and  let  1  be  the  stagnation  reference  state.   Then 


T. 


=   T   (hence  M0  =  M  =  any  value) 


T1   =   Tt  (hence  M  =  0) 


and  equation  *3»3,12  can  be  written  as: 


_T_ 
1\ 


1  +  *_I  m2 


=   f(M,y) 


(•3.4,1) 


As  another  example,  let  2  be  an  arbitrary  point  in  the 
flow  system  and  let  1  be  the  reference  state  *. 


Then: 


A2   =   A   (hence  M- 


M  =  any  value) 


A1      =   A*  (hence  IV^  =  1) 


Remembering  that  we  are  considering  an  isentropic  process, 
equation  *3«3>20  can  be  written  as: 


A 


1 

M 


1  +  ^|i  M2 

7+1 


Y  +  l 
2(y-D 


=   f(M,Y) 


(•3-4,2) 


We  see  that  for  isentropic  processes  these  important 
ratios  are  simply  functions  of  Mach  Number  and  gamma  (y) 
which  lead  to  tabulation  of  them  in  Appendix  A. 

Evaluate  equation  *3.4,1  for  an  arbitrary  value  of  M>  1 
and  an  arbitrary  value  of  M<1.   Compare  your  answers 
with  the  values  for  T/Tt  in  the  appendix.   Notice  that 
the  tables  are  dimensionless  and  are  always  expressed  in 
terms  of  a  ratio  of  pressures,  temperatures,  etc. 
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3.5  -  Relations  Between  P-^  and  A* 


a. 


We  can  now  relate  the  stagnation  pressure  (P+-)  and  the 
sonic  reference  area  (A*).   Take  A]_  as  A^*  and  take 
A2  as  A2*  in  equation  *3.3,20.   By  definition  IY^  =  M«  =  1 
This  expression  reduces  to 


A  * 
A2 

AT*" 


=   e 


AS 
R 


(•3.5,1) 


We  saw  that : 


t2 
tl 


=   e 


AS 

R 


(•3.3,15) 


Thus 


(£)(£)■  ( 


AS 
R 


AS 
'  R 


=   1 


(*3.5,2) 


or 


P    A  * 
rt2  R2 


P   A  * 
^tlAl 


constant 


(•3.5,3) 


Note  that  equations  *3-5,2  and  *3«5,3  apply  to  any 
adiabatic  situation  of  flow  of  a  perfect  gas  where 


ws  =  0. 


b.   It  should  be  noted  (and  we  can  verify  from  equation  *3-3,20) 
that  for  two  different  points  in  the  flow  field  (1  &  2), 
A]_*  is  only  equal  to  A2*  if  the  flow  is  isentropic  between 
points  1  &  2. 


Now  read  sections  3-5,  3.6  and  3. 7  in  the  text 
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UNIT  4  -  NORMAL  and  OBLIQUE  SHOCKS  -  Objectives 


The  student  shall  be  able  to: 

1.  List  the  assumptions  used  to  analyze  a  standing  normal 
shock. 

2.  Sketch  a  "shock  process"  on  a  T-S  diagram,  noting  the 
pertinent  features  of  the  phenomenom. 

3.  Explain  why  an  "expansion  shock  wave"  cannot  exist. 

4.  Recognize  the  necessary  and  sufficient  conditions  to 
solve  "normal  shock"  problems. 

5.  Derive  the  expression  which  describes  the  properties 
on  the  downstream  side  of  a  standing  normal  shock  in 
terms  of  properties  on  the  upstream  side. 

6.  Explain  how  shock  tables  may  be  developed  from  the 
relations  describing  the  properties  on  the  downstream 
side  of  a  standing  normal  shock  in  terms  of  the  properties 
on  the  upstream  side. 

7.  Demonstrate  the  ability  to  solve  typical  normal  shock 
problems  by  use  of  the  tables  and/or  normal  shock 
equations . 

8.  Explain  how  an  oblique  shock  can  be  described  by  the 
superposition  of  a  normal  shock  and  another  flow  field. 

9.  Identify  which  properties  remain  constant  and  which 
change  when  a  uniform  velocity  is  superimposed  on  the 
flow  field. 

10.  Show  (by  diagrams)  how  the  "shock  angle"  and  the 
"deflection  angle"  are  defined. 

11.  Describe  the  general  results  of  an  oblique  shock  analysis 
in  terms  of  a  diagram  such  as  "shock  angle"  versus 
"Mach  number"  for  various  deflection  angles. 

12.  Distinguish  between  weak  and  strong  shocks.   Know  what 
conditions  cause  each  to  form. 

13.  Describe  the  conditions  which  cause  a  detached  shock 
to  form. 

14.  Solve  typical  problems  involving  oblique  shocks  (such 
as  air  inlets,  nozzle  outlets,  wedges,  etc.) 
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15.   Describe  the  treatment  of  a  moving  normal  shock  wave 
so  as  to  apply  the  relations  developed  for  the 
standing  normal  shock. 
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Unit  4  -  NORMAL  and  OBLIQUE  SHOCKS 


k . 1  -  The  Shock  Mechanism 

a.  To  this  point  we  have  treated  fluid  flow  problems  which 
occured  in  a  continuous  medium;  that  is,  wherein  no  discon- 
tinuities existed  within  the  control  volume  under  consider- 
ation.  We  now  turn  to  consideration  of  a  discontinuous 
process  -  the  "shock".   A  shock  wave  characteristically 
occurs  within  a  very  thin,  but  finite  volume.   The  thick- 
ness of  a  shock  is  on  the  order  of  10~5  inches.   Due  to 

the  complex  interactions  involved,  analysis  of  the  pro- 
perties within  the  shock  are  beyond  the  scope  of  this  course. 
We  will  analyze  the  change  in  fluid  properties  across  the 
shock  by  investigating  the  behavior  on  both  sides  of  a 
shock  wave. 

b.  Now  read  Sections  4.1  and  4.2  in  the  text. 

c.  We  will  first  consider  a  "standing  (i.e.  NON-MOVING) 
normal  shock";  that  is,  a  standing  shock  which  is  perpendicular 
to  flow.   With  the  principles  developed  in  this  analysis 

we  will  then  investigate  the  possibility  of  the  existence 
of  an  expansion  shock.   Then  moving  normal  shocks,  conver- 
gent-divergent nozzle  operation,  and  oblique  shocks  will 
be  covered. 
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4.2-  Standing  Normal  Shock 


a.   Consider  a  flow  situation  as  shown  in  Figure  *4.1.   As 

with  previous  analyses  we  first  establish  a  control  volume, 


shock 


FLOW 


*■   x 


Figure  *4.1 


then  apply  the  basic  equations  to  the  properties  across 
the  control  volume. 


b.   The  following  assumptions  will  be  made: 

1)  Since  the  shock  is  very,  very  thin,  we  may  consider 
the  control  volume  equally  thin  and  assume  that 
there  is  no  surface  along  the  wall  to  create  friction 

2)  Due  to  the  thin  control  volume,  the  cross-sectional 
area  on  both  sides  of  the  shock  may  be  considered  to 
be  essentially  equal. 

3)  We  will  continue  to  assume  steady,  one-dimensional 
flow  with  negligible  potential  changes. 


c.   Consider  first  the  continuity  equation 
steady,  one-dimensional  flow: 

m   =   constant 


Since  this  is 


(•4.2,1) 


or 
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P1A1V1   ■   P2A2V2 


but 


Al      A2 


so 


p  V   =   («4.2,2) 


d.   Now  recall  the  energy  equation  for  steady,  one-dimensional 
flow  between  two  points: 


htx  +  Im   =   ht2  +  w  (*2.1,7) 

Sl->2 

We  can  consider  this  an  adiabatic  process  with  no  shaft  work 
being  done,  so   q,  0      =   w      =   0   ,  or 

l->2 


I 


h,    =   (»4.2,3) 

t]_  9 


e.   The  x-component  of  the  momentum  equation  for  steady, 
one-dimensional  flow  (Equation  1.13)  becomes: 

*— ^        c.v.      c 


Integrating  the  right  side  of  equation  *4.2,2  yields 


/^_  v   (V.dA)  =  —  [V0  -  Vn] 
J&c      x  gc    2     1 

We  can  also  express  the  summation  of  forces  in  the  x- 
direction  as 


XX   ■   P1A1  -  P2A2   "   (P1  "  VA 
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So  equation  *4  .  2 , 4  becomes 


(P   -  P2)A   -  JL  [V2  -  V,] 

&C 


(•4.2,5) 


but   m  =   pAV   so  show  that: 
2 


Pl  + 


»1V1 


=   P2  + 


p2v2' 


(*4.2,6) 


f.   To  this  point  we  have  been  dealing  with  a  general  fluid 
By  introducing  the  Perfect  Gas  Law 


P   =   PRT 


(1.23) 


and  recalling  the  expression  for  the  speed  of  sound  in 
a  Perfect  Gas  (Equation  *3.3,zO: 


V  =   Ma   =   M  \/  yg  RT 


(•4.2,7) 


We  can  now  express  Equation  *^.2,2  in  terms  of  Mach  number 


P  M 
rl11l 


P2M2 


V*i      V*i 


(*4.2,8) 


Obtain  this  expression  yourself  for  drill. 

;umed 


Similarly,  recalling  from  Equation  (3.6) 


T    =   T[l  +  ^  M2] 


(3.6) 


and  since 


ht  ■  cPTt 
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from  Section  1.1  in  the  text,  the  energy  equation 
(Equation  *4.2,3)  can  be  written 


T1[l  +  *-±   M!2]   =   T2^1  +  V"  M22] 


(•1.2,9) 


g.   Similarly,  using  the  definition  of  Mach  Number,  the 

Perfect  Gas  Law  and  the  expression  for  speed  of  sound  in 
a  perfect  gas,  show  that  the  Momentum  Equation  (Equation 
(4.2,6)  can  be  written: 


P1Cl  + 


ym12] 


P2[l  +  YM2  ] 


(*4.2,10) 


h.   We  have  in  effect  seven  variables 


Y,  P1,  M1,  T1,  P2,  M2  and  T2 


and  three  governing  equations,  nos.  *4 . 2 , 8 ,  *4.2,9,  and 
*4.2,10.   Given  any  4  of  these  we  should  be  able  to  solve 
for  the  remaining  three  variables. 


i.   We  now  find  it  possible  to  combine  the  three  basic  equations 
so  as  to  derive  an  expression  for  M2  in  terms  of  M,  and  7  . 
Rearrange  the  continuity  and  energy  equations  as  follows: 


P2M2 


(*4.2,8a) 


1/2 


T. 


i  +  i 


-l 


M 


2 


1     2    1 


(*^.2,9a) 


and  the  momentum  equation 


\  1  +  YM2' 


1  +  yM- 


(*4.2,10a) 
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Now  combine  these  three  and  show  that: 


1  +  YM2 


2  -i 


l  +  ym 


l  J 


M. 


M. 


1  +  Y  ^  M  2 
2   U2 

1  +  ^  Ml2  . 


1/2 


(*4.2,11) 


Solve  for  Mp  and  obtain 


M, 


2     2 

Mi  +7=T 
^)Mi2  "  x 


(*4. 2,11a) 


So  for  a  flow  situation  such  as  in  Figure  *4.13  knowing 
the  conditions  at  (1)  (i.e.,  before  the  shock)  and  solving 
for  M2  we  can  utilize  equations  *4.2,10a;  *4.2,7;  *4.2,6; 
*3-3,15;  *4.2,9a;  and  equations  3.6  and  3.7  from  the  text 
to  solve  for  all  other  properties. 

REPEAT!   If  the  properties  are  known  ahead  of  the 
shock  (condition  1)  then  ALL  conditions 
after  the  shock  (condition  2)  can  easily 
be  found. 

For  a  given  fluid  (i.e.,  y  Is  known)  since  the  property 
ratios  are  functions  of  y  and  M]_  only ,  it  is  possible  to 
pre-compute  the  ratios  of  the  fluid  properties,  and 
tabularize  them.   Appendix  B  in  the  text  is  such  a  table 
Take  M^  to  be  2.0  and  solve  for  M2 ,  and  the  property 
ratios,  using  the  relations  we  have  developed.   Compare 
your  answers  with  the  values  in  the  Normal  Shock  table. 

Notice  that  it  is  possible  to  use  a  subsonic  value  of  IVh 
and  obtain  a  MATHEMATICALLY  correct  answer  for  M2  which 
is  INVALID.   Why? 


HINT: 


Take  M-,  to  be  .577  and  solve  for 


t2 


tl 


Now  what  sort  of  value  must  As  (in  this  case 
Asi)  have  to  satisfy  Equation  *3-3,15?   This 
proves  that  an  "expansion  shock"  is  NOT 
POSSIBLE. 
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k.   Examine  Appendix  B  in  the  text.   Note  that  as  IVL  increases 
M2  decreases  and 

P2      p2  T 

p—  ,   —    and   fp—    are  all  >  1.0 

rl     pl  ll 

Therefore  a  shock  is  always  a  compression  process,  with 
?2   always  greater  than  P^.   Notice  that  Pt2  is  always 
less  than  P-^  (though  for  values  of  M^_  near  1.0  only 
slightly) . 

Thus  a  shock  process  is  IRREVERSIBLE. 
1.   Read  Section  4.3  in  the  text  and  study  the  examples. 
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4.3-  Moving  Normal  Shocks 


T   =   70  ^F 


V   =  400  fps- 


P   =   10  psig 


Figure  *4 . 2 


valve  open 


V  =  400  fps 


■&*■ 


s 


dir 

V  =  0 


valve  closed 

Figure  *4.2a 


T   =   ? 
P   =   ? 


a.   Now  consider  a  situation  as  pictured  in  Figure  *4.2  where 
air  is  flowing  in  a  pipe  at  the  conditions  shown  and  the 
valve  at  the  end  is  suddenly  closed.   This  will  cause  a 
shock  wave  to  propagate  back  through  the  duct  as  pictured 
in  *4.2a.   We  do  not  know  the  speed  of  this  shock  wave. 
Is  this  a  steady  flow  problem?   If  you  said  "Yes", 


GO  DIRECTLY  TO  JAIL 
DO  NOT  PASS  ? 


and  while  you're  there,  re-read  your  handout  entitled 
Comments  on  Steady  Flow. 
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b.   We  can  change  this  problem  into  a  steady  flow  problem  by 
"hopping  aboard  the  shock  wave." 


V   =  400  +  V 
1  s 


P-,  =  10  psig 
T1  =  530  °R 


V„  =  V 


■*- 


p    =   9 

2 

T    =   9 

i2 


Figure  *4.2b 


We  now  have  a  steady  flow  standing  normal  shock  problem. 

The  velocity  upstream  of  the  shock  wave  will  be  

and  the  speed  downstream  will  be  .   Why  can  we 

not  derive  a  so-called  closed  form  solution  for  Mp? 

c.   A  "trial  and  error"  process  is  required  to  solve  this 
problem.   Let  us  assume  that  IVL  =  1.3 


a. 


(l.M)  (32.2)  (53.3)  (530)     =   1128  ft/sec 


Thus  V^  =  M1a1  =  (1.3)  (1128)   =   1467  ft/sec 


and  Vs  =  V1  -  400   =   1067  ft/sec 


T 
2 
From  the  shock  tables  obtain  M0  =  .786  and  =-  =  1.191 


1 


Thus  T2  =  1.191  (530)   =   632  °R 
a2   =   1233  ft/sec 


and  V2  =  M2a2   =   970  ft/sec 
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But  V2  must  equal  Vs ,  which  it  does  not!   Thus,  our 
assumed  Mach  Number  was  in  error.   A  few  trials  will 
converge  on  an  assumption  of  Mj_   =  1.235-   Follow  these 
calculations  through  and  show  that  V2  =  Vs  =  992  ft/sec. 


What  would  the  pressure  be  behind  the  shock  wave?   How 

do  these  values  compare  with  those  in  the  actual  "unsteady" 

problem  posed  in  Figure  *4.2a? 


d.   Read  Sections  4.4,  4.5  and  4.6  in  the  text. 
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*J .  4  -  Convergent-Divergent  Nozzle  Operation 


a.   We  recall  that  a  convergent  section  acting  as  a  nozzle 
will  become  "choked"  at  M  =  1 . 


If  the  pressure  downstream  of  such  a  nozzle  is 
decreased  the  flow  will  accelerate  until  M  =  1 
is  reached  at  the  exit  (see  Figure  *4.3). 


UPSTREAM 


DOWNSTREAM 


'EXIT 


Figure  *4 . 3 


Further  decrease  in  pressure  outside  the  nozzle 
exit  will  not  induce  a  further  decrease  in 
pressure  inside  the  nozzle.   The  nozzle  is  said 
to  be  "choked". 


We  also  recall  (from  Figure  *3.1)  that  in  order  for  a 
device  to  function  as  a  'nozzle'  above  M  =  1.0  (i.e., 
continue  to  accelerate  the  flow,  thus  converting  more 
enthalpy  into  kinetic  energy)  the  cross-sectional  area 
must  begin  to  increase .   Thus  we  can  presage  the  existence 
of  a  convergent-divergent  section  (with  M  =  1.0  at  the 
throat)  which  functions  as  a  nozzle  throughout  its  entire 
length.   This  is  called  a  DeLaval  nozzle. 
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De  Laval  Nozzle 


First 
Critical 


Third 
Critical 


Figure  *H.h 


Entrance  Throat 


Exit 


The  convergent  section  of  the  C-D  nozzle  offers  little 
design  difficulty,  but  the  divergent  section  must  be 
carefully  designed  to  reduce  losses  due  to  separation, 
wave  interaction  and  other  fairly  complex  phenomena. 


b.   Below  the  diagram  of  the  C-D  nozzle  in  Figure  *4.4  is 

depicted  a  plot  of  static  pressure  ratio  versus  location 
in  the  nozzle.   V/e  can  envision  two  possible  flow  condi- 
tions with  M  =  1  at  the  throat.   Each  of  these  conditions 
is  associated  with  a  pressure  value  at  the  exit  known  as 
a  "critical  point."   The  "First  Critical  Point"  represents 
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flow  which  reaches  M  =  1.0  at  the  throat,  but  is  subsonic 
in  both  the  convergent  and  the  divergent  sections  (dotted 
line).   The  "Third  Critical  Point"  represents  the  design 
point  (or  pressure  ratio)  of  the  C-D  nozzle.   Both  1st 
and  3rd  critical  points  represent  isentropic  flow  condi- 
tions.  Any  exit  pressure  above  1st  critical  will  result 
in  subsonic  flow  throughout  the  nozzle,  or  typical 
"venturi"  operation.   A  point  between  1st  and  3rd  critical 
will  induce  non-isentropic  flow  through  the  nozzle. 


c.   Since  both  1st  and  3rd  critical  represent  isentropic  flow, 
Appendix  A,  "Isentropic  Plow  Tables",  may  be  used  to 
determine  the  values  of  Mach  Number  at  first  and  third 
critical,  provided  we  know  the  area  ratio  between  the  exit 
and  the  throat . 

As  an  example  suppose: 

A   .  . 
BX^ =   3.0  (*4.4,1) 


throat 


We  know  that  at  the  throat  M  =  1.0  so  A,,    ,  =  A*,  or: 


A 
A* 


-   3.0  (*JJ.4,la) 

Prom  Appendix  A  in  the  text : 


for  first  critical  M  ,.  ~.197  w~     »  -973    (*4.4,2) 

6X11/  4- 

for  third  critical  M   .,  -  2 . 64  £-     =  .0^72   (*4 .4 ,3) 

GX1I/  x  «_ 

t 


Two  Notes  on  this  solution 

1)  Linear  interpolation  can  be  used  to  obtain 
better  accuracy. 

2)  The  small  figure  (-1)  indicates  the  power 
of  ten  with  which  the  tabular  value  must  be 
multiplied. 

It  is  also  interesting  to  compare  these  results  to  that 

of  a  converging  (only)  nozzle.   In  the  case  of  a  converging 

nozzle  with  y  =  1 . ^  we  need  a  pressure  ratio  of  .528. 
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(P   ,,  to  P.  .  . 
exit     inlet 


)  to  produce  sonic  velocity.   In  the  above 


converging-diverging  nozzle  we  have  sonic  velocity  (in  the 
throat)  with  a  pressure  ratio  of  only  .97. 


As  you  might  expect,  somewhere  between  "FIRST  BASE"  and 
"THIRD  BASE"  there  must  be  a  "SECOND  BASE".   Let's 
consider  what  happens  in  the  region  between  "first  critical" 
and  "third  critical." 


As  the  receiver  pressure  is  lowered  belox^  "first  critical," 
a  normal  shock  forms  just  downstream  of  the  throat  (see 
Figure  *4.5).   The  remainder  of  the  "nozzle"  is  now  acting 
as  a  diffuser  since  the  flow  is  subsonic  and  the  area  is 
increasing.   The  shock  will  locate  itself  in  a  position 
such  that  the  pressure  changes  that  occur  ahead  of  the 
shock,  across  the  shock,  and  downstream  of  the  shock  will 
produce  a  pressure  that  exactly  matches  the  outlet  pressure. 
In  other  words,  the  outlet  pressure  determines  the  location 
and  strength  of  the  shock.   As  the  pressure  is  further 
lowered,  the  shock  continues  to  move  toward  the  exit.   When 
the  shock  is  located  at  the  exit  plane  this  condition  is 
referred  to  as  the  2nd  critical  point.   If  the  receiver 
pressure  is  between  2nd  and  3rd  critical,  then  a  compres- 
sion takes  place  outside  the  nozzle.   This  is  called 
"overexpansion"  (i.e.,  the  flow  has  expanded  too  far  within 
the  nozzle).   If  the  receiver  pressure  is  below  3rd  critical 
then  an  expansion  takes  place  outside  the  nozzle.   This 
condition  is  called  "underexpansion . "   We  shall  investigate 
these  conditions  later  in  the  course. 


Entrance   Throat 


Exit 


First  Critical 


Second  Critical 


Third  Critical 


Normal  Shock 
Inside  Nozzle 


Figure  *4.5 
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f .   We  can  now  use  what  we  know  about  standing  normal  shocks 
to  find  the  conditions  at  '2nd  critical'.   Just  ahead 
(e  away",  as  the  mathematicians  are  wont  to  say)  of  the 
shock  at  2nd  critical,  the  conditions  are  identical  to 
those  at  3rd  critical  (or  design).   If  you  have  not  done 
so,  use  the  tables  to  confirm  the  exit  Mach  at  3rd  critical 
is  M  =  2.64  for  an  Area  ratio  of 


exit 


A 


=      3.0 


throat 


FLOW 


£6 


Shock  - 

2nd  critical 


Figure  *4.6 


exit 
> 

inlet 


tl 


If  we  assume  stagnation 
conditions  at  the  inlet 


Is  this  a  reasonable  assumption? 
From  the  Tables : 


t2 


P 


(*4.4,4) 


tl 


From  Isentropic  Flow  Tables  %   M  =  2.64 


From  Normal  Shock  Tables  @  M  =  2.64 
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Compute  this  and  verify  that 


exit 
inlet 


.376 


Now  suppose  we  subject  the  same  De  Laval  nozzle  to  a 
pressure  ratio  of  0.4,  whre  will  the  shock  form?   Is  this 
'underexpanded'  or  ' overexpanded '  Flow?   Think  about  it 
before  checking  the  solution  below. 


Solution 

Since  the  flow  is  between  1st  and  2nd  critical,  a  normal 
shock  has  formed  inside  the  nozzle  —  hence  the  terms 
"underexpanded"  and  "overexpanded"  do  not  apply  to  the 
flow. 


Figure  *4.7 


inlet 


exit 


We  have  now  enough  information  to  analyze  the  flow, 
including  losses. 


Known 


exit 
throat 


=   3.0 


exit 
inlet 


0.4 
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But:     P       =   p 
inlet      t 


inlet 


and     A.,    .  =   A*   ,  . 

throat      inlet 


since  we  assume  no 
losses  between  the 
inlet  and  the  throat. 


Now  recall  Equation  *3.5,3.   Therefore: 


A  *     P 
inlet   rt 


inlet 


A  *    P, 

exit    t 


exit 


=   Constant     (*4.4,4) 


Thus 


A 


exit 
Lthroat 


P  P., 
exit 


inlet 


'  t     t 

known  known 

i     i 

(3.0)  (0.4) 


A   ., 

exit 

A  *   . 
inlet 


exit 


t 


inlet 


A   ,. 

exit 

A  * 

exit 


exit 


_  exit_ 


1.2 


or: 


A   ,, 

exit 

A  * 
exit 


exit 


t 


exit 


=   1.2 


From  the  equations  developed  in  Unit  3: 


A* 


=   f(M,Y) 


(*3.4,2) 


4-      =   f(m,Y) 

r  t 


(3.7) 
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You  can  obtain  an  expression  for 


p-  jr-   =  f(M,7)  («4.4,5) 


John  does  not  tabulate  this  value  but  a  sample  table 
including  this  function  was  passed  out  in  class.   Excerpts 
from  a  similar  table  for  y  =  1.4  show: 

_A_  ?_ 
M  A*   Pt 


.47  1.205 

.48  1.178 

Interpolation  reveals  the  Mach  No.  at  the  exit  to  be  about 


M   .    &    .472 
exit  -     ' 


If  we  can  determine  the  stagnation  pressure  ratio 

P      /p        we  will  be  able  to  locate  the  shock 
exit'    inlet 


P  /P      \  /  P 

t   ..         t   ..  \   exit 
exit         exit 


Ptinlet     l?exlt  /lPtinlet'       ^859 

\ 

N  Given 

From  Isentropic  tables  for  M  =  .472 


Since  all  losses  occur  across  the  shock  we  can  say  that 


Pt  Pt 

exit  =   2   i  (*H.4,6) 

p^  Pt 

inlet         1  /   across  shock 
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From  the  shock  tables  find  M,  =  2.59  and  from  the 
isentropic  tables  find 


Asho_ck_  =   2>g69 


throat 


We  also  have  a  measure  of  the  loss  from 


P,  As 

--£   =    e"  "R" 

Pt 


As  we  recall  that   As  =  As.  here. 


h.   Read  Section  5.2  and  5.3  in  the  text. 


(•3.3,15) 
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4.5  -  Oblique  Shocks 


a.   Read  Section  6.1  in  the  text 


b.   Consider  a  Standing  Normal  Shock 


© 


V 


In 


V 


2n 


Recall:   using  Shock  as 
Reference 


Vn   >  Sonic 
In 

Vp   <  Sonic 


Figure  *4  .  8 


Now,  let  us  superimpose  a  velocity  to  V^n  &  v"2n  (say  Vt ) 
Note  that  this  is  the  equivalent  of  running  along  the 
shock  front .   Recall  that  by  doing  this  we  have  not 
changed  the  static  states  of  the  fluid,  but  we  have 
changed  the  stagnation  conditions.   We  then  have: 


V 


Since  V2n  <  Vln 
and  Vt   =  Vt 


Notice  that  V"2  has  turned 
towards  the  shock  front. 


Figure  *4 . 9 
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We  normally  picture  this  as  follows 


Shock 
0      Angle 


0-8 


§   Flow 

Deflection 


Angle 


Figure  *4.10 


NOTE   (1)   Flow  direction  has  changed 

(2)   V}  >  Vo  as  with  normal  shock 

V]_  must  be  supersonic  (since  V]_n  is  supersonic 
V"2  can  now  be  supersonic  —  if  V-^  is  large 
enough 

It  should  be  clear  that  the  Normal  Shock  Relations  (and 
Tables)  can  be  used  IF,  AND  ONLY  IF,  proper  care  is  taken. 


Dividing  both  sides  by  a-,  ,  and  noting  that  a-,  =  a-,n>  we 
have  : 


V 


In 
Lln 


V- 


sin  6 
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Recalling  the  definition  of  Mach  Number,  this  becomes 
M    =   M   sin  6 

Static  pressure  ratios,  static  temperature  ratios,  etc.  — 
can  now  be  taken  from  the  Normal  Shock  Tables,  since  these 
were  unaltered  by  the  superposition  of  V,  on  the  original 
shock  picture. 

c.   Now,  we  know  that  Mn   >  1 

'  In  — 

thus :    M-.  sin  e  >_  1 

Let  us  consider  the  range  of  possible  e  for  a  given  M: 
The  minimum  e  will  occur  when  IYL  sin  e   =1 

6  .    =   sin"   77- 
min  M, 


Recall  that  this  is  the  same  expression  that  was  developed 
for  the  Mach  Angle  y.   Hence,  the  Mach  Angle  is  the 
minimum  possible  shock  angle.   Note  that  this  is  a  limiting 
condition  and  really  no  shock  exists  for  this  case  since 
then  M]_n  =  1.0.   The  maximum  value  that  8  can  achieve  is 
obviously  90°.   This  is  another  limiting  condition  and 
represents  our  familiar  normal  shock. 

Notice  that  as  the  shock  angle  6  decreases  from  90°  to  the 
Mach  angle  y,  M^n  decreases  from  Mj  to  1.   Since  the 
strength  of  a  shock  is  dependent  upon  the  normal  Mach 
number  we  have  the  means  to  produce  a  shock  of  any 
strength  equal  to  or  less  than  the  normal  shock. 

c.  Study  Example  6.1  on  pages  105-106  of  the  text. 

d.  We  will  novj  try  to  relate  the  deflection  angle  (6)  to  the 
shock  angle  (e)  starting  with  the  continuity  equation. 

PlVln  =  p2V2n  Ci..5,«) 
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so 


2n 
In 


(•4.5,5) 


and  from  equation  (*4.5,1)  and  figure  6.5 


Vt   =   V1  cos 


Vp  cos  (e  -  6) 


(•4.5,6) 


so 


V, 


COS  6 


cos(e-5 ) 


sin  9 
sin(6-6 ) 


tan  6 
tan(e-6) 


(•4.5,7) 


From  our  Normal  Shock  Relations  in  Unit  4.2,  we  can  sh 

that  for  a  Perfect  Gas  (as  a  function  of  M,  ): 

In 


ow 


(y+1)  M 


In 


(Y-1)  \n      +    2 


(•4.5,8) 


but   M,   =  M,  sin 
In    1 


(*4.5,8a) 


Combining  equations  *4.5,7,  *4.5,8  and  *4.5,8a  yields 


(y+1)  Mi   sin' 
(y-1)M12  sin2i 


+  2 


tan  9 
tan(9-6) 


(•4.5,9) 


So   9   =   f(6,  M,  y) 

or  for  a  given  M,  and  y >    the  shock  angle  (e)  is  a  function 
of  the  deflection  angle  (6).   But  note  that  this  is  a 
transcendental  function.   So  for  supersonic  flow  along  a 
wall  we  cannot  solve  for  9  as  a  function  of  6 ,  M  and  y 
EXPLICITLY,  but  we  can  obtain  an  explicit  solution  for 
6  =  f(6,  M,  y)« 
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tan  6=2  cot 


M-,2  .  2 

1  sin  i 


M,  (y+cos  26)  +  2 


(•4.5,10) 


Hence  we  can  construct  a  chart  such  as  in  Appendix  C  to 
relate  these  three  in  accordance  with  Equation  *4.5,10. 


e.   Now  read  the  remainder  of  Section  6.2,  and  carefully 
study  the  examples. 


f.   Read  Section  6.3  in  the  text. 


g.   Read  Section  6.5.   We  summarize  here  some  of  the 
important  concepts : 

1.  Flow  always  turns  "toward"  an  oblique  shock  front. 

2.  For  given  values  of  6  and  M,  two  values  of  0  may 
exist . 

a.  A  large  pressure  ratio  results  in  a  strong 
shock  and  subsonic  M2 . 

b.  A  small  pressure  ratio  results  in  a  weak  shock 
and  supersonic  M2 . 

3.  A  maximum  value  of  6  exists  for  a  given  Mach  Number. 

4.  For  6  >5max  ,  a  detached  shock  results  (see  Figures 
6. 8  and  6.9  in  text ) . 
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UNIT  5  -  PRANDTL-MEYER  FLOW  -  OBJECTIVES 


The  student  shall  be  able  to: 

1.  Know  how  Entropy  Changes  and  Pressure  Ratios  vary  with 
deflection  angles. 

2.  Explain  how  finite  turns  (with  finite  pressure  ratios) 
can  be  accomplished  isentropically . 

3.  Show  Prandtl-Meyer  flow  (both  expansions  and  compressions) 
on  a  T-S  diagram. 

H.      Demonstrate  the  development  of  the  relation  between 
MACH  number  (M)  and  flow  turning  angle  (v). 

5.  Show  hov;  tables  can  be  developed  for  Prandtl-Meyer  Flow 
by  the  introduction  of  a  reference  state. 

6.  Explain  the  governing  boundary  conditions  and  show  the 
results  when  shock  waves  and  P-M  waves:  (a)  reflect  off 
physical  boundries,  and  (b)  reflect  off  "free" 
boundaries . 

7.  Draw  the  wave  forms  created  by  flow  over  rounded  and/or 
wedge-shaped  wings  as  the  angle  of  attack  changes.   Be 
able  to  solve  problems  of  this  type. 

8.  Describe  and  sketch  what  occurs  as  fluid  flows  past  a 
smooth  concave  corner  and  a  smooth  convex  corner, 
explaining  why  an  expansion  shock  CANNOT  occur. 

9.  Explain  with  the  aid  of  diagrams  the  flow  conditions  at 
the  outlet  of  a  supersonic  nozzle  created  by  under- 
expansion  and  overexpansion. 

10.  Solve  typical  problems  using  PRANDTL-MEYER  flow  tables. 

11.  Demonstrate  by  a  sketch,  an  understanding  of  the 
operating  conditions  of  a  fixed-geometry  diffuser. 
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UNIT  5  -  PRANDTL-MEYER  FLOW 


5.1  -  Introduction  to  Prandtl-Meyer  Flow 


You  saw  from  the  "Special  Problem"  that  for  values  of  M, 
(in  a  NORMAL  shock  process)  which  are  close  to  1.0,  the 
value  of  the  quantity  (M]_2  -  1)3  which  we  will  call  'm' , 
is  small.   It  was  shown  that  for  small  values  of  'm', 
the  pressure  change  across  the  shock  (AP)  is  proportional 
to  the  first  power  of  'm' .   We  also  saw  that  the  entropy 
change  (  As)  is  proportional  to  the  third  power  of  'm' . 


Now  consider  equation  *4.5,9  from  the  last  Unit,  but 
inverted. 


tan(e-6 ) 
tan   e 


(Y-l)M12sin2e  +  2 
(Y+l)M12sin2  6 


(*4.5,9a) 


We  can  rearrange  this  to  show: 


1 


2    2 
M,  sin  9 


(v+l)tan(e-6) 
2tan  6 


(y-D 


(»4.5,9b) 


By  applying  the  auspicious  trignometric  relations  and 
re   ranging  we  can  show  that  this  becomes: 


M12sin20  -  1 


(y+D  m  sin  9  sin  6 


M 


cos(e-6 ) 


(*5.1,D 


Notice  that  for  small  values  of  <5  where  we  can  approximate 
sin  6  as  6,  and  cos(6-<j>)  as  cos  6: 


M12sin26  -  1 


il+li  Ml2  tan 


(*5.1,2) 


You  saw  from  equation  *4.5>8  as  applied  in  the  "Special 
Problem" : 


P2  -  P, 


2y 


Y  +  1      1 


(m/  -  1) 


(•5.1,3) 


115 


Note  that  M1  in  equation  *5.1,3  represents  the  normal 
mach  number  (M1  )  and  since  from  unit  4.5(b)  we  saw 


Mln   =   Ml  sin 


Thus  Equation  (*5.1,3)  as  applied  to  oblique  shocks  is 


P2  -  P 


^j-  (M12sin2e  -  1) 


(•5.1,4) 


Thus  we  can  now  see  that  for  small  deflection  angles  (and 
hence  weak  shocks)  combining  Equations  (*5.1,1|)  and 
(*5.1,2)  yields: 


yM 


1 


(M12-l)1/2 


(*5.1,5) 


c.   We  can  draw  the  conclusion  from  this  that  for  very  weak 
oblique  shocks  the  pressure  change  across  the  shock  is 
proportional  to  the  deflection  angle. 


In  summary:   AP 


m 


where :   m  =  (M 


In 


1) 


But  also    AP  <*  m 
and        As  ^  m 


Thus 


As 


AP" 


Now  suppose  we  let  the  deflection  angle  become  infinitely 
small.   What  happens  to  As?   If  as  is  essentially  zero 
for  an  inf intesimally  small  turn,  we  can  accomplish  a 
finite  turn  isentropically  by  a  series  of  such  infinites- 
imally  small  turns,  and  the  finite  turn  will  have  a 
finite  pressure  change. 

Furthermore,  since  this  situation  has  no  losses,  it 
is  reversible .   Thus  we  can  have  compressions  or  expansions 
depending  on  the  boundary  conditions. 
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d.   Now  read  Sections  7.1,  7.2  and  7 • 3  in  the  text.   Note 

carefully  the  reasoning  which  leads  to  the  Prandtl-Meyer 
expansion  fan  pictured  in  Figure  7 . 8 . . . and  it  is  ISENTROPIC 
In  this  isentropic  (or  Prandtl-Meyer)  expansion  process: 

1)  Mach  number  Increases . 

2)  Mach  angle  Decreases . 
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5.2-  Flow  Equations 


a.   Read  section  1 .h    in  the  text,  through  page  125 


b.   Notice  that  the  tables  are  referenced  to  the  boundary 
condition  where  M0  =  1.0  and  v  =  C 
is  the  Mach  number  AFTER  the  turn. 


condition  where  M0  =  1.0  and  v  =  0.   Also  notice  that  M 


M.=  1.0 


Figure  *5.1 


Nov;  consider  Figure  7-10  in  the  text.   Suppose  you  were 
given  M]_  =  2.13  and  asked  to  find  M2 ,  given  the  angle 
V2~V]_  =  ^0°.   From  the  Prandtl-Meyer  flow  tables 
Appendix  D,  you  can  find  the  "virtual"  angle  v\y 


would  have  produced  a  Mach  number  of  M]_ 
then  add  this  value  (30°)  to  the 
consult  the  table  f or  v =  70°,  we 


value  of 
find  M2  : 


2.13 

VI  -v 

■   5.3: 


m 

which 
If  we 
and 


d.   Since  this  is  an  isentropic  process,  we  know  that  Ptl=Pt2 
and  Ttl=Tt2  >    so  we  ™ay  utilize  the  ISENTROPIC  tables 
(appendix  A)  to  solve  for  static  pressure  and  temperature 


Now  read  the  remainder  of  Section  7 .  ^    in  the  text  and 
follow  through  the  examples  carefully. 


Read  Section  7.5  in  the  text.   Note  that  for  the 
compression  process  the  assumption  of  Prandtl-Meyer 
flow  (hence  Isentropic)  is  only  valid  in  a  limited  region 
close  to  the  wall.   The  discussion  on  concave  corners 
(Section  7.3  in  the  text)  shows  that  for  many  situations 
involving  a  smooth  concave  turn,  the  flow  characteristics 
away  from  the  boundary  or  wall  can  be  found  by  considering 
the  turn  as  a  sharp  concave  corner.   The  transition  from 
point  to  point  close  to  the  wall  approaches  the  condition 
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of  a  series  of  Mach  waves  as  shown  in  Figure  7.5,  but  as 
we  consider  points  away  from  the  wall,  we  find  the 
infinitesmally  small  compression  waves  coalesce  into  a 
finite  oblique  shock.   Figure  7.6  in  the  text  is  not  an 
accurate  d  piction  of  the  coalescing  process.   If  the 
area  close  to  the  wall  is  enlarged,  as  in  Figure  *5.2. 


Figure  *5.2 


We  find  an  envelope  of  Mach  lines  forms  as  the  Mach  lines 
coalesce  forming  a  finite  shock  wave  at  the  Inclination 
Angle,  (e)  corresponding  to  the  overall  deflection  (6) 
and  the  initial  Mach  number.   Since  boundary  layer  and 
other  real  gas  effects  at  the  wall  are  often  neglected 
in  basic  engineering  considerations,  we  often  consider 
only  the  flow  away  from  the  wall.   However  if  the  turn 
is  sufficiently  smooth  and  gradual,  the  Prandtl-Meyer 
flow  region  may  be  significant,  and  the  flow  within  a 
reasonable  distance  from  the  boundary  may  be  treated  as 
Prandtl-Meyer  flow  with  reasonable  accuracy. 
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5-3  -  Wave  Reflections 


a.   Read  Sections  7.6  and  7.7  in  the  text 


b.   Review  the  summary  in  Section  6.5  of  the  text  which  puts 
great  emphasis  on  the  importance  of  boundary  conditions. 
Let  us  now  consider  the  two  principle  physical  conditions 
that  govern  wave  behavior.   We  shall  accomplish  this  by 
analyzing  the  operation  of  a  converging-diverging  axi  - 
symmetric  nozzle  depicted  in  Figure  *5.3  which  is  operatini 
between  its  SECOND  and  THIRD  critical  points.   Obviously 
at  a  physical  wall,  the  flow  must  be  parallel  to  the  wall. 
We  may  view  the  central  streamline  as  though  it  were  a 
"wall",  hence  (if  we  treat  the  flow  as  one  dimensional 
we  may  consider  the  flow  in  region   1   to  be  parallel  to 
the  centerline  as  depicted.   For  the  second  condition 
consider  the  boundary  of  the  "free  jet"  as  depicted  by 
the  dotted  line  in  region  (2)  .   There  must  be  pressure 
equilibrium  along  this  boundary,  so  ?2   ~   ^ambient-   We 
can  now  follow  from  region  to  region  and  by  matching  the 
required  pressures  or  flow  direction  angles,  we  may 
determine  a  great  deal  more  about  the  flow  conditions. 


free  jet  boundary 


Figure    *5«3 


p_  p_  p_  p 

2      lA  _  "6      'amb 

P>P       >P 

"3      'amb     "5 


l?>Pa 


amb 
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c.  If  we  operate  this  nozzle  with  a  pressure  ratio  between 
2nd  and  3rd  critical  it  is  obvious  that  we  need  a  compres- 
sion process  at  or  beyond  the  exit  in  order  for  the  flow 
to  end  up  at  the  ambient  pressure.   However,  a  normal 
shock  at  the  exit  will  produce  too  strong  a  compression. 
What  is  needed  is  a  shock  process  that  is  weaker  than  a 
normal  shock  and  the  oblique  shock  has  been  shown  to  be 
just  this.   Thus,  at  the  exit  we  observe  oblique  shock 

(A)  at  the  appropriate  angle  so  that  P~  =  P 

d.  We  recall  that  the  flov;  across  an  oblique  shock  is  always 
deflected  toward  the  shock  and  thus  the  flow  in  region  (2) 
is  no  longer  parallel  to  the  centerline.   Wave  front  (B) 
must  deflect  the  flov;  back  to  its  original  axial  direction 
This  can  easily  be  accomplished  by  another  oblique  shock, 
(an  alternate  way  of  viewing  this  is  that  the  oblique 
shocks  from  both  upper  and  lower  lips  of  the  nozzle  "pass 
through  each  other"  when  they  meet  at  the  centerline.   If 
one  adopts  this  philosophy  one  should  realize  that  the 
waves  are  altered  in  the  process  of  traveling  through  one 
another. ) 


e.  Nov;,  since  ?2   =  ^ambs  passage  of  the  flow  through  oblique 
shock  (B)  will  make  P3  >  Pamb  and  region  (3)  can  not  have 
a  free  surface  in  contact  with  the  surroundings.   Conse- 
quently a  wave  formation  must  emanate  from  the  point  where 
wave  (B)  meets  the  free  boundary  and  the  pressure  must 
decrease  across  this  wave.   An  "expansion  shock"  would 
fulfill  our  requirements  but  we  know  that  no  such  animal 
can  exist.   We  now  realize  that  wave  form  (C)  must  be  a 
Prandtl-Meyer  expansion  so  that  ? u    =  P  ,  . 

f.  However,  passage  of  the  flow  through  the  expansion  fan  (c) 
causes  it  to  turn  away  from  the  centerline.   Thus,  as  each 
wave  of  the  P-M  "expansion  fan"  meets  the  centerline  a 
wave  form  must  emanate  to  turn  the  flow  parallel  to  the 
axis  again.   If  wave  (D)  were  a  compression  in  which 
direction  would  the  flov;  turn?   We  see  that  to  meet  the 
boundary  condition  of  flow  direction  wave  (D)  must  be 
another  P-M  expansion.   Thus  the  pressure  in  region 
is  less  than  ambient. 

g.  Can  you  now  reason  that  to  get  from  (5)  to  (6)  and  meet 
the  boundary  condition  imposed  by  the  free  boundary,  (e) 
must  consist  of  P-M  compression  waves.   Similarly  (F) 
must  consist  of  P-M  compression  waves  In  order  to  turn 
the  flov;  in  region  (7)  to  match  the  direction  of  the 
"wall".   Now,  is  P7  equal  to,  greater  than,  or  less  than 
Pamb?   You  should  realize  that  conditions  in  region  Q) 

are  similar  to  those  in  region  (J)  and  so  the  cycle  repeats 
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The  situation  described  above  (operation  between  2nd  and 
3rd  critical)  is  referred  to  as  OVEREXPANSION.   It  is  also 
described  in  section  8.3  in  the  text. 


From  this  example  we  may  draw  some  general  conclusions 
about  reflections: 

1)  Reflections  from  a  Physical  (or  pseudo-physical 
boundary  where  the  boundary  condition  concerns 
the  flow  direction)  are  the  same  "family".   That 
is,  shocks  reflect  as  shocks,  compression  waves 
reflect  as  compression  waves,  and  expansion  waves 
reflect  as  expansion  waves. 

2)  Reflections  from  a  "free!"  boundary  (where  pressure 
equilization  occurs)  are  of  the  opposite  family, 
i.e.,  compression  waves  reflect  as  expansion  waves, 
and  vice  versa. 


The  angle  through  which  the  flow  turns  from  (3)  to  (J) 
is  determined  by  the  pressure  change  required,  whejreas 
the  angle  from  (5)  to  (5)  depends  on  the  flow  in 
paralleling  the  "wall". 


k.   Let  us  now  examine  an  UNDEREXPANDED  nozzle, 
operation  below  3rd  critical. 


This  is 


Figure  *5-^ 
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Note  that  the  flow  leaving  this  nozzle  has  a  pressure 
greater  than  ambient  and  the  flow  is  parallel  to  the 
axis.   Reflect  back  to  Figure  *5.3  showing  the  overexpanded 
nozzle  and  you  will  see  that  the  condition  is  exactly  the 
same  as  region  (J)  •   Thus  the  flow  patterns  are  the  same 
from  this  point  on.   This  situation  is  also  described  in 
Section  8.3  of  the  text. 


1.   We  should  note  that  at  the  jet  edge  in  a  real  gas  a 

SHEAR  layer  exists  as  a  result  of  the  velocity  gradient 


m.   Now  read  Section  7.8  in  the  text. 
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5.4  -  Supersonic  Diffuser 

a.  We  will  now  examine  the  design  of  a  fixed  geometry 
supersonic  diffuser.   Prom  our  studies  of  supersonic  nozzles 
we  would  assume  that  a  converging-diverging  section  would  do 
the  trick  -  and  indeed  it  will.   However,  there  are  some 
practical  operating  difficulties  that  must  be  considered. 

b.  Suppose  we  design  the  inlet  diffuser  for  an  airplane 
that  will  fly  at  Mach  1.86.   From  the  isentropic  tables  we 
see  that  the  area  ratio  corresponding  to  this  Mach  No.  is 
1.507.   We  shall  construct  the  diffuser  with  an  area  ratio 
(inlet  to  throat)  of  1.500.   Now  let  us  follow  the  operation 
of  this  diffuser  as  the  aircraft  takes  off  and  accelerated 
to  its  design  speed.   Note  that  as  the  flight  speed  reaches 
M   *  .43  the  diffuser  becomes  choked.   (Check  the  subsonic 
portion  of  the  isentropic  tables  for  the  above  area  ratio.) 
This  condition  is  shown  in  Figure  *5.5a.   Now  increase  the 
flight  speed  to  M0  =  .6.   "Spillage"  of  external  diffusion 
occurs  as  indicated  in  Figure  *5-5b.   As  M0  is  increased  to 
1.0  there  is  a  further  decrease  in  the  "capture  area" 

(area  of  the  flow  at  the  free  stream  mach  number  that 
actually  enters  the  diffuser). 

c.  As  we  increase  M0  to  supersonic  speeds  a  detached 
shock  wave  forms.  See  figures  *5«5d  and  *5-5e.   Note  that 
at  the  higher  flight  speeds  the  shock  moves  closer  to  the 
inlet  as  less  external  diffusion  is  required  to  reach 

M]_  =  .43  at  the  inlet.   Also  note  that  it  is  necessary  to 
fly  at  slightly  greater  than  MQ  =  4.19  in  order  to  have  the 
shock  attached  to  the  inlet  as  shown  in  figure  *5.5f 
(Check  the  shock  tables  to  substantiate  this.)   If  we  now 
increase  MQ  to  4.2  the  shock  moves  very  rapidly  past  the 
throat  and  forms  in  the  divergent  section  downstream  of 
the  throat  as  shown  in  figure  *5-5g.   This  is  referred  to  as 
"swallowing  the  shock"  and  the  diffuser  is  said  to  be 
"started".   Under  these  conditions  we  no  longer  have  Mach 
1.0  in  the  throat.   (Can  you  compute  the  Mach  number  that 
exists  in  the  throat?)   We  can  now  slowly  decrease  the  flight 
speed  to  the  design  condition  of  M  =  1.86  and  the  shock  will 
move  to  a  position  just  downstream  of  the  throat  and  occur 
at  a  Mach  number  of  just  slightly  greater  than  1.0.  (Thus 
we  have  a  very  weak  shock  and  negligible  losses.)   See 
figure  *  5«5h. 

d.  Two  comments  can  now  be  made  on  the  above  performance: 

1)   In  order  to  "start"  the  diffuser,  which  was 
designed  for  1.86,  it  is  necessary  to  "over- 
speed"  the  plane  to  a  Mach  number  of  4.2. 


124 


2)   If  the  plane  slows  just  slightly  below 
its  design  speed  —  or  perhaps  minor  air 
disturbance  might  cause  it  to  drop  below 
1.86  -  the  shock  will  pop  out  in  front  of 
the  inlet  and  the  diffuser  must  be  "started" 
all  over  again. 


in  (a)  thru  (f) 
M2=  1.0 


0 


Mr  .43 


9 


(a) 


M,=  43 


(b) 


Figure    *5«5 


M0=l. 


(C) 


M  =  1.2 


Ml 


(d) 


M0  =  2.0 


M,=43^- 


(e) 
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M0=4.19 
M,=  .43 


(f) 


Mo>4.20 


1.86  <  M<4.20 


(9) 


(h) 


e.   It  should  be  obvious  that  either  of  these  situations 
cannot  be  tolerated  and  for  this  reason  one  does  not  see 
fixed  geometry  converging-diverging  dif fusers .   We  can 
not  summarize  the  behavior  of  fixed  geometry  diffusers 
as  is  depicted  in  figure  *5-6. 


shock  'swallowed' 


Figure  *5-6 


condition  to 
swallow  shock 


design  area    ratio 


tached  shock 
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5.5  -  Airfoils 

a.  Read  Section  8.5  in  the  text. 

b.  Follow  carefully  through  the  examples  in  section  8.5 


127 


UNIT  6  -  FANNO  FLOW  -  OBJECTIVES 


The  student  shall  be  able  to: 

1.  List  assumptions  and  restrictions  applicable  to  FANNO 
FLOW  and  relate  these  to  real  physical  situations. 

2.  Sketch  a  Fanno  line  in  the  h-s  and  T-s  planes,  identifying 
the  regions  of  sonic,  supersonic  and  sub-sonic  flow. 

3.  Describe  the  variation  of  static  and  stagnation  pressure, 
static  and  stagnation  temperature,  density  and  velocity 
as  flow  varies  along  a  FANNO  line  (for  both  subsonic 

and  supersonic  flow). 

*J .   State  what  is  meant  by  "choked  flow". 

5.  Describe  the  effects  of  changing  duct  length  in  a  choked 
Fanno  flow  situation  in  supersonic  and  subsonic  flow. 

6.  Define  "friction  factor",  "equivalent  (or  hydraulic) 
diameter",  "absolute  and  relative  roughness",  "absolute 
and  kinematic  viscosity",  and  "Reynolds  number",  and 
know  how  to  determine  each. 

7.  Starting  with  basic  principles,  derive  expressions  for 
property  ratios  such  as  T2/T]_,  P2/Pl»  etc.  in  terms  of 
Mach  number  and  specific  heats,  in  a  Fanno  flow  situa- 
tion,  (for  Honors  credit) 

8.  Describe  (including  h-s  diagram)  how  the  Fanno  tables 
are  developed  by  use  of  a  reference  location. 

9.  Demonstrate  the  ability  to  solve  typical  Fanno  flow 
problems  by  use  of  the  appropriate  tables  and/or 
relations . 
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UNIT  6  -  FANNO  -FLOW  -  STUDY  GUIDE 


6.1  -  Introduction  to  Fanno  Flow 

a.  We  have  to  this  point  touched  only  briefly  on  the  subject 
of  friction  losses.   This  unit  will  analyze  the  effects  of 
flow  in  a  constant  area  duct  with  friction.   The  results  of 
this  analysis  have  many  direct  engineering  applications  in 
compressible  fluid  flow,  and  will  pave  the  way  for  increased 
understanding  of  the  effects  of  friction  on  fluid  flow 
under  any  circumstances.   To  simplify  the  analysis  we  will 
assume  that  the  flow  is  adiabatic  and  no  shaft  work  is 
added  to,  or  extracted  from  the  flow.   This  assumption  is 
generally  quite  valid  when  dealing  with  reasonably  short 
ducts  where  no  special  attempt  is  made  to  transfer  heat  or 
work.   Additionally,  we  shall  assume  steady,  one-dimensional 
flow  with  no  appreciable  change  in  potential.   This  type  of 
Flow  is  known  as  "Fanno  Flow". 

b.  We  shall  first  consider  the  effect  of  these  assumptions 
on  the  basic  equations.   Then  the  general  variation  of  fluid 
properties  with  Mach  number  will  be  considered.   Finally, 
working  equations  for  Fanno  flow  with  a  perfect  gas  will  be 
developed,  from  which  we  may  see  how  values  are  tabulated 
with  respect  to  a  reference  condition  as  has  been  done  with 
shocks  and  isentropic  flow.   Thus  we  will  simplify  the 
business  of  problem  solving.   For  your  own  benefit,  summarize 
the  assumptions  that  we  are  making  in  Fanno  Flow. 
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6.2-  Basic  Development  of  FANNO  Fl 


ow 


a.   Consider  first  the  continuity  equation  for  steady  one- 
dimensional  flow  as  developed  in  Unit  1: 

m  =   PAV  =  (£  (*6.2,1) 

Since  the  cross-sectional  area  is  constant,  this  reduces 
to: 

Pv  =  (j;  (*6.2,2) 

We  will  assign  a  new  symbol,  "G"  to  the  quantity  ' pV 
which  is  referred  to  as  the  "mass-velocity".   So  by 
definition : 


PV   e   G  (*6.2,3) 

What  are  the  units  of  "G"? 


b.   We  will  now  resurrect  the  energy  equation,  which  for 
steady  flow  between  any  two  points  in  the  flow  may  be 
written  in  terms  of  stagnation  enthalpy: 

htl   +   ^2   "  ht2   +   W1^2  (*2-1>7) 

but  since 


ql+2   =   Wl->2   =   ° 


then 

h,   «  h.    =  h,  =  (t  (*6.2,4) 

t]_        t2        t     T 

We  can  write  this  equation  in  terms  of  static  conditions 
as : 

ht   .  h  +  g-     +  g  -   <f  (.2.1.4) 
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Nov;  assuming  the  potential  remains  essentially  constant 
and  utilizing  equation  *6.2,3  we  can  rewrite  equation 
*2.1,H  as: 


P  2g, 


=   h, 


"  4 


(•6.2,5) 


Since  for  a  given  flow,  G  and  h,  are  constant,  we  see  that 
for  Panno  flow  h  and  p  are  uniquely  related.   We  can  now 
plot  this  relationship  for  various  values  of  G  (these  are 
called  Fanno  Lines)  in  the  h-v  plane  (recalling  that 
v  =  1/p): 


s3>  s2>s 


t    »  « 

\ 

v      v    » 

\     v 

.  \ 

v     v 

v          \        » 

\ 

\ 

\           \ 

\ 

\  - 

~^^^^~"==:^----~^     \ 

^ \S 

x/\\ 

V         V         \       -*J=5    (t 

V  \*    \    ^  ^ 

\v       \o    -^       \              ^ 

\          \                  XZ.    "-» 

\o  W       V 

1          \o 

! L \_ >. 

Figure  *6.1 


v(=  Vp  ) 


If  the  fluid  is  known,  one  can  also  plot,  as  in  Figure 
*6.1,  lines  of  constant  entropy  (dotted  lines)  which  come 
from  thermodynamic  considerations.   We  can  then  replot 
the  "Fanno  lines"  in  the  more  familiar  h-s  plane  as  shown 
in  Figure  *6.2. 
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h 


ht=$ 


Figure  *6.2 


c.   We  may  now  state  in  general  that  the  combination  of  the 
continuity  equation  with  the  energy  equation  and  the 
thermodynamic  equation  of  state  (which  established  the 
lines  of  constant  entropy)  form  the  locus  of  points  on  a 
given  Panno  Line . 

One  factor  stands  out : 


Since  this  is  ADIABATIC,  i.e.,  dSe  =  0,  the  only 
way  entropy  can  be  generated  is  by  LOSSES,  so  a 
process  can  only  move  toward  INCREASING  values  of  S. 


-p* 
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6.3  -  FANNO  FLOW  versus  MACH  NUMBER 


a.   We  shall  note  that  due  to  the  irreversible  nature  of 

FANNO  Flow,  a  FANNO  LINE  consists  of  two  distinct  regions 
separated  by  a  Limiting  Point. 


h 


Decreasing 
^  Pressure 


limiting  point 

Figure  *6.3 
Fanno  line 


We  note  that  in  the  upper  branch  a  process  (which  must 
always  result  in  INCREASING  ENTROPY)  results  in  DECREASED 
PRESSURE.   From  Figure  *6.1,  we  say  what  as  enthalpy 
decreases,  so  must  density  decrease,  hence  velocity  must 
INCREASE! 


b.   Now  consider  for  yourself  what  happens  on  the  LOWER 

branch.  Draw  an  arrow  on  the  Lower  Branch  of  Figure  *6.3 
to  indicate  which  way  a  process  must  proceed.  As  we  move 
along  this  branch  of  the  Fanno  line: 


Enthalpy 
Density 
Velocity 
Pressure 
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c.   Let  us  now  consider  the  Limiting  Point .   We  are  seeking 

to  explore  a  condition  where  ds  =  0.   Recalling  the  energy 
equation  (neglecting  Potential  changes). 


\    =   h   +   4-  -  £  C*2-1'^ 


Differentiate  this  to  obtain: 

0  =  dh  +  ^  (*6.3,1) 


Also  recall  that 

PV  =  (j:  (*6.2,2) 

Differentiating  this  we  obtain: 

pdV  +  Vdp   =   0  (*6.3,2a) 


or 


dV  =   _  Vdfi.  (^6. 3,2b) 

P 


Substitute  equation  *6.3,2b  into  Equation  *6.3,1  and  solve 
for  dh: 


dh   = 

g 


™p_  (*6.3,3) 

c 


Recall  the  thermodynamic  property  relation 

Tds   =   dh  -  vdp  (*6.3,^a) 


Tds   =   dh 

p 


^£  (*6.3,i»b) 
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Combining  Equations  *6.3,3  and  *6.3>^b  we  find  that 


(*6.3,5) 


This  expression  is  valid  at  ANY  point  on  a  Fanno  Line. 
We  now  apply  equation  *6.3,5  to  the  limiting  point  where 
ds  =  0,  then  we  find  at  the  limiting  point: 


IE 
P 


V2dp 
gcp 


or 


V 


.  d£ 
'c  dp 


At 

Limit 

Point 


8p_ 

gC  3p 


s=constant 


We  recognize  this  as  the  expression  for  a   (see  section 
2.3  in  the  text).   Hence  at  the  limiting  point: 


V  =  a 


or 


M  =  1.0 
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6.4  -  WORKING  RELATIONS  for  FANNO  FLOW  in  a  PERFECT  GAS 

a.   From  the  continuity  equation  in  steady,  one-dimensional 
flow  in  a  constant  area  duct: 

G  =  pV  =  <£        .  (*6.2,2) 

From  the  Perfect  Gas  Law: 


_P_ 

RT 


and  from  the  definition  of  Mach  number 


V  =   Ma 
and  for  a  Perfect  Gas 


a  =  aJ  YgcRT 
Equation  (*6.2,2)  can  be  re-written: 

=    (j:  (*6.4,i) 


PM 


V 


RT 
Ygc 


b.   Now  we  turn  to  the  Energy  Equation.   We  showed  that 
between  two  points  in  a  Fanno  Line: 


h.    =   h,  (*6.iJ,2) 

tl  t2 


Recall  that  for  a  perfect  gas: 


\  =  Vt  (*6-"'3) 
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So 


T, 


T    =  Tt  =  <£  (*6.H,4) 


We  also  saw  that  for  a  Perfect  Gas 


T.   =   T(l  +  *zi  M2)  (3.6)  in  text 


2 


so  between  tv;o  sections 


T.d+^M,2)   =   T9(1  +  ^M92)  (*6.455a) 


2    1  '      X2XJ-     2    2 


or  as  a  temperature  ratio 


T2 

1    +   *|i  f^2 

Tl 

1    +   x^i  M    2 

(*6.4,5b) 


c.   By  expressing  the  relation  in  Equation  *6.4,1  between  two 
points  in  the  flow  we  can  write 


P  M      PpM 

(•6.4,6a) 


yTh        V^ 


This  can  be  restated  as  a  ratio  of  the  static  pressures : 

1/2 

-1  =   _1   _2  (*6.4}6b) 

Px     M2  \  T  /  ' 
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or,  utilizing  equation  *6.4,5b  we  can  now  obtain  an 
expression  for  a  ratio  of  the  static  pressures  in  terms 
of  y   and  Mach  Number: 


1/2 


M. 


Y  -1    2-i 
1  +  - — -  M 

2    1 

Y  -1     2 

1  +  - — -  M 
-J-     2   n2  J 


(•6.4,7) 


Have  you  seen  some  of  these  equations  before?  For  what 
kind  of  flow?  What  assumptions  have  you  fed  into  these 
equations  so  far? 


d.   Now  recall  the  property  relation 


Tds   =   dh  - 


dP 


What  restrictions,  if  any,  are  attached  to  this  expression? 
If  we  now  recall  from  Section  1.11  in  the  text  that: 


dh   =   c  dT 


(*6.4,8) 


and 


C    =   R[-\] 
P  Y-l 


(*6.4,9) 


we  can  show  that 


Tds   =   [  -V  ]  RdT  -  — 
y-l  P 


(*6. 4,10a) 


and  we  can  rearrange  this  to  show 


R        Y-l 


dT 
T 


dP 
P 


(*6.4,10b) 
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Now  integrate  Equation  *6.4,10b  between  two  points  in  the 
flow  to  show  that 


S2"S1 
R 


T2 

7=i    ln  % 


In 


(•6.3,10c) 


Utilizing  the  expressions  we  have  just  developed  for  the 
static  pressure  ratio  and  the  static  temperature  ratio 
we  can  now  express  the  entropy  change  as  a  function  of  y 
and  Mach  Number: 


1/2 


«       q 

1   +   ^^  Mn  2 

2       1 
R 

=      -*-=-  in 

Y-1 

2         1 
2      u2 

ln 


1  +  ^^  M 
x     2    1 

1  +  ^Z=-   M 


(*6. 4,11a) 


You  can  now  combine  the  terms  of  this  expression  to  show 
that 


S0-S-,        M 
_± — ±  =   in  — 

R  Mx 


1  +  ^^  M 
2    1 

1  +  I|i  M22 


y    l 

y-1  "  2 


(*6. 4,11b) 


or 


S2-Sl 

R 


=   ln 


M, 


1  +  y~ -  M  ' 
-1     2    1 

1  +  ^  M2: 


y  +  1 

2T7=ry 


(6.4,11c) 


Refer  back  to  the  section  where  we  developed  the  stagnation 
pressure  ratio  as: 


AS 
R 


(•3.3,15) 


You  should  satisfy  yourself  that  this  relation  applies  to 
Fanno  Flow. 
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If  we  now  convert  Equation  *6.4,llc  to 


S2"S1 
R 


=   In 


M. 


1  +  %^=-  M 

1     2  2 

1  +  ^CZ-L  M 

2  1 


Y  +  l 
2(y-D 


(*6.4,12) 


Now  take  the  anti-log  of  both  sides  to  obtain 


As. 
R 


t2 


M. 


1  +  ^~-   M22 
1  +  ^~   M  2 


-,   Y  +  l 

2i7^ry 


(•6.4,13) 


We  now  have  the  means  to  obtain  all  the  properties  at  a 
downstream  point  (g)  if  we  know  all  the  properties  at 
some  upstream  point   1   and  the  Mach  Number  at  point  \2) 
But  we  need  to  be  able  to  predict  the  Mach  Number  at  (2) 
from  the  physical  set-up  of  the  flow  (i.e.,  Duct  length, 
etc .  ) . 


e.   In  Unit  1,  we  developed  an  expression  for  the  Momentum 

Equation  as  it  applies  to  steady,  one-dimensional  flow  in 
any  fluid 


dP    fdxV 
P    D  2g 


+  -£-  dz  + 
gc 


VdV 


=   0 


(*1.6,16) 


If  we  now  apply  the  perfect  gas  law  and  the  definition 
of  Mach  Number  to  this  equation,  we  can  show  (neglecting 
Potential  change): 


Y-   <RT)   + 


fdx  Myg^RT 


D 


2g, 


MVygpRT"  dV 
Sc 


0   (*6.4,l4) 


140 


This    can  be   rewritten   as: 


.2 


f     +      ™*   M_y      +      xdM2      +      XM2    dT      =      0  (t6.i»,15) 

e 


These  two  terms  result  from 
eliminating  dV  in  the  last  term 
of  Equation  *6.4,l4.   Can  you 
show  this? 


f.   We  must  first  evaluate  dP/P  as  a  function  of  Mach  Number. 
Prom  continuity: 


™   =  $  C6.4.6) 


v^ 


Therefore : 


In  P   +   In  M  -  ~   In  T   =   In  (j;  (*6.4,l6a) 


Differentiate  to  show: 


f  +  f  -  |f  -  0  (.6.4.16b) 


Now  we  must  obtain  an  expression  for  dT/T  as  a  function 
of  Mach  number.   From  the  energy  equation  we  obtained: 


Tt 


or 


In 


T(l  +  *~  M2)   =  <£ 


T   +   ln(l  +  ^M2)   =   ln(j;  (*6. 4,17a) 
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Again  differentiate   to   obtain 


dT  d(l    +  1=1  m2) 

V      +      v-1      2  =         ° 

1    +   ^-  1VT 


(*6.4,17b) 


Now  combine  equation  *6.4,17b  with  equation  *6.4,l6b 
and  substitute  into  the  momentum  equation  in  the  form 
of  equation  *6.4,15  to  show: 


fdx 
D 


d(l  +  Y~1  M2) 
(  ill  )      f   2   n  ; 

2Y      !  +  izi  M2 


2dM 

2 

YM 


(  I±i  )  dM' 


2y 


(*6.4,l8) 


M' 


We  can  now  integrate  equation  *6.4,l8  from  point  (T)  to 
point  (2)  in  the  flow  (see  Figure  *6.4)  to  obtain: 


f(x2-xl}    y+1  - 
D =  27  ln 


1  +  X-l  m22 " 

1  +  £==■  M 

-1  +   2    1  J 


1 


_M2     M 


1_1    y+1 

2     2Y 
1  J 


M, 


In 


M. 


T 


(•6.4,19) 


Figure  *6.4a 


Flow 


© 


x,  — ♦ 


©    © 


M=  1 


Figure  *6.4b 
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Knowing  the  duct  length  (x2  -  x]_),  flow  rate,  fluid  and 
duct  material  (hence  the  friction  factor),  configuration 
(to  obtain  De )  and  all  conditions  at  one  point,  we  may 
now  find  the  conditions  at  the  other  point. 

g.   The  solution  of  problems  (as  in  previous  analyses)  is 

simplified  by  the  introduction  of  a  reference  condition. 
The  obvious  choice  is  the  limiting  condition  where 
M  =  1.0.   If  vie  let  point  (1)  be  the  "*"  condition  (i.e., 
M  =  1.0;  reached  by  Fanno  Flow)  then: 


1 
P1  =  P* 

Point  (2)  may  be  any  other  point  in  the  flow.   Equation 
*6.^(,5b  may  now  be  reduced  to: 


1+1 

T    _       2 

T*"  y-1   2 

1       1  +  ^~  M 


(*6.4,20) 


Using  a  similar  approach,  obtain  the  expressions  for: 


pr  =     ?  (*6.4,2l) 


and 


P, 
— L  =       ?  (*6.iJ,22) 

V 


1^3 


You  should  also  be  able  to  corroborate  that 


1/2 


p 

V* 

V 

1 

M 

"  2(1   +   ^zl  M2) 

^  ' 

y    +    1 

(«6.i»,23) 


h. 


Now  check  your  expressions  for  P/P*  and  P,/P.  *  by 
solving  for  any  arbitrary  value  of  Mach  (subsonic  or 
supersonic)  and  using  a  value  of  y  =  1 . ^ .   Then  compare 
your  solutions  against  the  tabularized  values  in  Appendix 
E  in  the  text. 


M 


6.5-  Applications  of  Fanno  Fl 


ow 


a.   Applying  the  reference  condition  to  equation  *6.4,19 
(i.e.,  M1  =  1,  M2  =  M)  yields: 


f(x-x*) 
D. 


=      I 


+  1 


n  +  xzi  m2i 


In 


1 


+  1 


4llnM2-I 


-    1 


M 


(*6. 4 ,2^; ) 


Since  x*  is  always  greater  than  x  (see  Figure  *6.4b), 
it  makes  sense  to  change  all  signs  in  Equation  *6.4,24 
to  obtain  a  positive  expression  which  we  can  simplify  to 


f(x*-x) 
D 


=  I 


+  1 


2y 


In 


1+1  m2 

2 


1  +  .^~-  M' 


M 


-  1 


(*6.4,25) 


The  quantity  (x*-x)  now  represents  the  maximum  duct  length 
which  may  be  added  to  the  given  duct  and  still  maintain 
the  same  flow.   This  expression 


f  L 


max 


D. 


is  tabulated  in  Appendix  E  in  the  text 


b.   In  order  to  solve  problems  we  must  be  able  to  calculate 
friction  factor  (f)  and  equivalent  (or  "hydraulic") 
diameter  (De).   From  Unit  1  you  may  recall  the  definition 
of  these  two  quantities. 


PV< 


D 


For  a  circular  duct,  D   is  simply  the  diameter  (D). 


11J5 


The  friction  factor  may  seem  a  bit  intractible  since  the 
shear  force  (t-)  is  difficult  to  measure.   However,  you 
may  recall  from  your  previous  courses  in  Fluid  Mechanics 
that  friction  factor  is  a  function  of  Reynolds  number  (Nr) 
and  "relative  roughness"  which  we  can. define  as  e/De  where 
"e"  is  an  average  value  for  the  absolute  roughness  of  a 
given  material.   Table  *6.1  gives  some  typical  values  of 
"e"  (which  has  the  dimension  of  "feet")  for  various 
materials.   Thus  the  relative  roughness  (e/De)  is  a 
dimensionless  parameter. 


TABLE  *6.1 
TYPICAL  ROUGHNESS  VALUES 

Material  "e"  (in  'feet') 

glass,  brass,  copper,  lead  smooth  <  .0001 

steel,  wrought  iron  .00015 

asphalted  cast  iron  .0004 

galvanized  iron  .0005 

cast  iron  .00085 

wood  staves  .002 

concrete  .01 

riveted  steel  .03 


d.   Similarly,  Reynolds  number  'is  a  dimensionless  parameter 
which  relates  size  and  viscosity.   Reynolds  number  is 
defined  as: 


N„   =   £ 


VD 


R      »*c 

where  y  is  "absolute  viscosity" 

If  we  define  a  "kinematic  viscosity" 


ygc 
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Then 


NR   " 


VD 
v 


Table  3  on  page  387  in  the  text  gives  typical  values  for 
"vg   "    f°r  various  gases.   With  this  information,  and 
knowing  the  conditions  at  a  given  point  in  the  flow,  we 
can  calculate  the  density  and  velocity  (remember  G  =  pV 
=  (£ )  hence  Reynolds  number.   Then  using  experimentally 
derived  curves  such  as  depicted  in  Figure  *6.5,  the 
friction  factor  may  be  determined. 


-Laminar  flow  -  f  ~ 


64 


Figure  *6.5 


(log  SCALE) 


e.   Now  read  sections  9.1  and  9.2  in  the  text.   Study  care- 
fully example  9.1. 


f.   Now  consider  the  following  example  problem: 

Given  a  circular  duct  with  these  specifications 
Diameter  -  6  inches 
Length   -  100  feet  (Ax) 

Material  -  Galvanized  Iron 

-* L-. 

max 


Figure  *6 . 6 


1^7 


It  is  desired  to  have  this  duct  deliver  800  cubic 
feet/minute  of  air  at  70°F  and  18  psia.   Find  the  inlet 
conditions  required  at  point  (T)  . 

First,  let  us  summarize  a  step-by-step  problem  solution 
technique : 

1)  Sketch  the  physical  problem. 

2)  Establish  and  designate  locations  where  conditions 
are  known  and/or  desired. 

3)  Compute  the  Equivalent  Diameter. 

4)  Find  friction  factor. 

5)  From  equation  *6.4,19  (or  the  Tables  using  FL   /D, 
if  known)  find  Mach  number (s). 

6)  Use  relations  or  Tables  to  find  additional  properties 
desired. 


Now  let  us  apply  this  approach  to  our  example  problem. 
Steps  1)  and  2)  are  done  in  Figure  *6.6. 
Recall  for  a  circular  duct  that  the  Equivalent  Diameter 
is  the  pipe  diameter. 


D   =   .5   feet 
e 

From  the  Table  *6.1  for  Galvanized  Iron  the  value  for 
e/D  is 


D 


To  find  Reynolds  number  we  must  know  p  and  V.   Since  we 
know  the  conditions  at  (2)  ,  we  can  compute 

V2 


M2      =      5 


2 


a, 


p 


^YgcR  T2     =     49<y/5^0"    =     1128     ft/se. 

p 
4t       =      .0918      lbm/ft3 


2  R  ^ 


D2  TT  ~  2 


l2  "T"  15" 
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ApVp   =   Volumetric  flow  rate 


So 


V2   =  ft/sec 


V2 
M0   =   —   =   .602  X      .60 
2     a2 


NR 


D 


Now  you  should  be  able  to  determine  from  the  attached 
"Friction  Factor  for  Pipes"  chart: 


f   =   .0197  -   .02 


So 


fAx      (.02)(100)       h    nn 
"d-  "  UW~       ~     4*°° 


f  L 
-From  the  Appendix  E,  for  Mp  =  .6,     max   _    liqnRl 


thus 


f  Lx     f  L2    fAx 
"~1T~   =   ~~D~~  +   ~ 


Now  consult  Appendix  E  to  find   M, 
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Now  knowing  Mj_   and  M2  you  can  use  the  property  ratios 
to  obtain  the  properties  at  point  (l)  . 


T   T* 

T    =   -^  —  T9   =   556  °R    ;    P. 


Note  that  we  do  not  have  to  bother  with  subscripts  on 
the  *  quantities  (such  as  T*  and  P*)  since  for  any  given 
Fanno  Line  there  is  only  one  reference  location. 


g.   Now  read  section  9-3  in  the  text.   Pay  close  attention 
to  the  flow  situations  as  described  in  Figure  9.16. 
We  will  not  delve  into  the  subject  of  Iso-Thermal  Flow, 
but  if  time  permits,  you  are  encouraged  to  finish  reading 
Chapter  9  in  the  text. 
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Unit  7.  -  RAYLEIGH  FLOW  -  OBJECTIVES 


The  student  shall  be  able  to: 

1.  State  the  assumptions  and  restrictions  utilized  in 
the  analysis  of  constant  area  flow  with  simple  heat 
transfer  (Rayleigh  Flow). 

2.  Given  the  general  equations  of  energy,  continuity, 
momentum  and  state,  demonstrate  how  the  assumptions  in 
objective  (1)  modify  or  eliminate  each  term. 

3.  Sketch  a  Rayleigh  line  in  the  P-v  plane  together  with 
lines  of  constant  entropy  and  lines  of  constant  temper- 
ature.  Indicate  the  points  of  maximum  temperature 

and  maximum  entropy,  and  the  directions  of  increasing 
entropy  and  temperature. 

h .      Sketch  a  Rayleigh  line  on  a  T-s  diagram,  indicating 
the  regions  of  Sonic,  Subsonic  and  Supersonic  flow. 

5.  Describe  the  changes  in  fluid  properties  which  occur 

as  one  moves  along  a  Rayleigh  line  for  the  case  of  heating 
and  the  case  of  cooling  in  both  subsonic  and  super- 
sonic flow. 

6.  Correlate  the  reference  points  of  maximum  temperature 
and  maximum  entropy  on  a  T-s  diagram  of  a  Rayleigh 
line  v/ith  those  on  a  P-v  diagram. 

7.  Explain  what  is  meant  by  "thermal  choking". 

8.  Explain  by  a  T-s  diagram  how  the  flow  adjusts  to  the 
addition  of  heat  in  a  constant  area  duct  in  a  "thermally 
choked"  flow  situation  in  both  subsonic  and  supersonic 
flow. 

9.  Solve  typical  Rayleigh  flow  problems  by  the  use  of 
tables  and  equations. 
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Unit  7.  -  RAYLEIGH  FLOW  -  STUDY  GUIDE 


7 . 1   Introduction  to  Rayleigh  Flow 

a)  Constant  area  flow  with  significant  heat  transfer, 
addition  or  removal,  is  called  "Rayleigh  flow".   We 
will  approach  the  analysis  of  this  flow  condition 
in  a  manner  similar  to  that  used  in  previous  flows 
such  as  normal  shocks,  Fanno  flow,  etc.   We  will  first 
develop  the  basic  equations  as  applied  to  the  Rayleigh 
flow  situation.   Then  perfect  gas  relations  will  be 
introduced.   Finally,  we  will  introduce  a  reference 
state  which  permits  tables  to  be  developed,  which  will  in 
turn  simplify  problem  solution  techniques. 

b)  Read  Section  10.1  in  the  text. 
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7.2  -  Basic  Relations 


a.   We  can  summarize  the  assumptions  made  in  Rayleigh  flow 
as  follows : 

1)  Steady  flow 

2)  One-dimensional  flov; 

3)  Negligible  change  in  potential  energy 

4)  No  work  (shaft) 

5)  Constant  area 

6)  Negligible  friction 

This  last  assumption  warrants  some  comment  since  internal 
irreversible  effects  (dSj.)  are  obviously  present.   What 
is  really  being  said  is  that  since  we  are  dealing  with 
significant  heat  transfer  the  change  in  entropy  from  this 
cause  (dSe )  will  be  relatively  large.   Thus  the  assumption 
that  is  actually  being  made  is  that: 


dS^  >>  dS. 
e      i 

and  we  can  thus  neglect  dS.  and  say  that 

dSe  «   dS  (*7.-2>l) 

b.   From  the  continuity  equation,  we  will  recall 
pAV   =   rn  =  <£ 
but   A  =  (t  hence  we  can  say 

pV  =  cj: 

Do  you  recall  this  expression?   Do  we  have  a  name  for 
this  constant? 


c.   The  energy  equation  for  steady,  one-dimensional  flow  can 
be  written  as: 


h,   +  q     =   h    +  w  (*7.2,2) 

tl     1    d  Z2  Sl+2 
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For  our  flow        =  0  ,  and 

Sl->2 


t1    ^l-*2       tp 


(*7.2,3) 


A  Note  of  Caution: 

This  is  the  first  major  flow  category  for 
which  'h^'  has  not  been  constant.   Consider 
a  perfect  gas.   Will  'T  '  be  constant? 


d.   Recall  our  basic  momentum  equation  in  differential  form 


dP   +    fdx_Vl  +  VdV  +   J£_ 
P  De   2gc         gc         sc 


dz   =   0 


(*7.2,4) 


Under  our  assumptions  this  becomes: 


dP   VdV  _   n 
—  +  —   u 

P    g„ 


(•7.2,5a) 


or 


dP  +  p 


VdV 


=   0 


(*7.2,5b) 


From  continuity  we  know  that  pV  =  G  =  constant 
the  momentum  equation  can  be  written  as 


Thus 


dP  +  S2V  =   0 


(•7.2,6) 


This  can  easily  be  integrated  to  give 


P  + 


GV 

S„ 


(*7.2,7) 


151! 


Or,  by  substituting  for  the  velocity,  this  can  be 
written  as 


P  + 


G 


gcp 


Constant 


Cf7.2,8) 


or 


P  +  £-  v 


=   Constant 


(•7.2,9) 


For  a  given  flow  rate  this  indicates  a  linear  relationship 
between  P  and  v.   Try  sketching  equation  (*7.2,9)  in  the 
P-v  plane.   This  line  is  called  a  "Rayleigh  line." 


Hint: 

y  =  mx  +  b 


Figure  *7 • 1 


As  an  aside,  rewrite  equation  ( *7 .2,7)  as 

PA  +  Olav)_v     =   Constant 

Sc 


or 


PA   + 


m  V 
So 


Constant 


Compare  this  expression  with  equation  (4.2,5)  in  the  unit 
on  shocks.   This  constant  is  sometimes  called  the  "impulse 
function"  and  is  given  the  symbol  "F"  or  "I"  by  various 
authors . 

If  we  add  lines  of  constant  temperature  to  our  plot 
(this  is  easy  if  we  are  dealing  with  a  perfect  gas) 
we  will  have  a  diagram  as  shown  in  Figure  *7.2. 
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Figure  *7-2 


v 


Consider  a  heating  process  which  moves  from  A  to  B  along 
the  Rayleigh  line.   If  we  add  more  heat,  we  move  farther 
along  and  the  temperature  increases  more.   When  point  C 
is  reached,  the  fluid  has  reached  a  maximum  temperature. 
Can  we  move  farther  along  this  Rayleigh  line?   What  is 
the  limiting  factor? 

Recall  that  as  we  add  heat  we  are  increasing  the  entropy 
of  the  fluid  (dSe  =  dq/T).   Thus  we  had  better  investigate 
lines  of  constant  entropy.   When  these  lines  are  plotted 
on  Figure  *7.2  they  are  very  similar  to  the  lines  of 
constant  temperature.   To  get  an  accurate  picture  of  what 
this  plot  will  look  like  we  must  investigate  the  relative 
slopes  of  'T  =  Constant'  and  'S  =  Constant'  lines  in  the 
P-v  diagram.   It  will  be  much  easier  to  carry  out  this 
investigation  if  we  consider  perfect  gases. 

For  a  ' T  =  Constant'  line 


Pv 


RT  =  constant 


Pdv  +vdP 


and 


dP 
dv 


F 
v 


(*7.2,10) 
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For  an  'S  =  Constant'  line 


Pv   =   constant 


vYdP  +  PyvY  1dv   =   0 


and 


dP 
dv 


=   -  Py 


V 


Y-l 


V 


P 
"  Y  V 


(•7.2,11) 


Comparing  equations  (*7.2,10)  and  (*7.2,11)  and  noting 
that  y  is  always  greater  than  1.0,  we  see  that  the 
'  S  =  constant'  line  has  a  greater  negative  slope  and 
thus  these  lines  will  plot  as  shown  in  Figure  *7-3. 


Figure  *7-3 


g.   We  now  see  that  we  can  add  more  heat  and  go  past  the 
point  of  maximum  temperature  until  we  reach  the  point 
of  maximum  entropy.   Let  us  pause  for  a  moment  and  see 
if  this  appears  reasonable. 

Normally  we  think  of  the  addition  of  heat  as  causing  the 
fluid  density  to  decrease.   This  requires  the  velocity 
to  increase  since  pV  =  constant.   The  increase  in  velocity 
requires  a  certain  increase  in  kinetic  energy.   Thus, 
some  of  the  heat  which  is  added  goes  into  increasing  the 
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kinetic  energy  of  the  fluid  and  some  increases  the 
enthalpy  (temperature)  of  the  fluid.   As  we  add  more  heat 
the  required  velocity  increase  (and  corresponding  kinetic 
energy  increase)  becomes  larger  and  larger.   Eventually 
we  reach  a  point  where  all  of  the  heat  energy  added  is 
needed  for  the  kinetic  energy  increase.   At  this  point  we 
reach  a  maximum  enthalpy  (temperature).   Further  addition 
of  heat  causes  the  kinetic  energy  to  increase  by  an  amount 
greater  than  the  heat  energy  being  added.   Thus  from  this 
point  on  the  enthalpy  (temperature)  decreases  to  provide 
the  proper  energy  balance. 


h.   Can  we  go  past  the  point  of  maximum  entropy?   Not  by  heat 
addition.   An  alternate  way  of  stating  this  is  that  no 
more  heat  can  be  added  to  the  system  if  we  wish  to 
maintain  the  same  flow  rate  (i.e.,  stay  on  the  same 
Rayleigh  line ) .   Thus,  we  see  that  a  limit  point  has  been 
reached  and  we  say  that  the  flow  is  "choked."   We  next 
proceed  to  investigate  this  limit  point. 

At  the  point  of  maximum  entropy  the  slope  of  the  's  = 
constant'  line  is  the  same  as  that  of  the  Rayleigh  line. 
Recall  the  equation  of  the  Rayleigh  line : 

P  +  =   Constant 

Differentiate  to  obtain 


(«7.2,12) 


dP   _   -^  =   -p2V2 
dv      gc        gc 


Recall  the  slope  of  the  's  =  constant'  line: 


|f  =   -  y|  =   -  YP2RT  (*7.2,11) 


We  now  equate  these: 

2  2 
-  £_£-  =   _  yp2rt  (*7.2,13a) 

gc 

and  solve  for 

V2   =   gcyRT  (*7.2,13b) 
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You  immediately  recognize  this  as  sonic  velocity  and 
thus  the  'limiting  point'  divides  the  Rayleigh  line  into 
a  subsonic  branch  and  a  supersonic  branch  (similar  to 
the  situation  in  Fanno  flow) .   Subsonic  flow  is  on  the 
left  and  supersonic  flow  is  on  the  right  of  the  limiting 
point  in  the  P-v  diagram. 


Another  interesting  fact  can  be  noted  to  exist  at  the 
limiting  point.   Recall  the  entropy  relation  from 
Section  1.10: 


TdS   =   dh  - 


dP 


(1.22) 


The  differential  form  of  equation  (10.2)  may  be  written 


dP   =   - 


pVdV 
Sc 


(•7-2.14) 


The  property  relation  in  equation  (1.22)  can  thus  be 
expressed  as : 


TdS   =   dh 


VdV 


(•7.2,15) 


At  the  'thermal  choke  point'  where  M  =  1,  dS  =  0,  so: 


0   =   dh  + 


VdV 


(•7.2,16) 


Integrate  this  to  obtain: 
r2 


h  + 


\r 


2g, 


Constant 


(•7.2,17) 


We  have  already  shown  that  this  is  the  expression  for 
'stagnation  enthalpy'  (neglecting  potential). 


h, 


=  h  + 


2g, 


4 


(•7.2,18) 
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Since  at  this  'limiting  point',  h,  =  constant,  dh,  must 
equal  zero.   Thus  h,  is  a  maximum  at  this  point. 

k.   The  above  situations  become  more  clear  if  the  Rayleigh 

line  is  plotted  in  the  h-s  plane  along  with  its  associated 
stagnation  curves.   This  has  been  done  in  Figure  *7.^« 


\./ 


Figure  *7.H 


Rayleigh  Line 


Now  correlate  this  plot  with  the  previous  information 
and  satisfy  yourself  that  the  static  and  stagnation  lines 
are  properly  drawn.   Compare  this  figure  to  Figure  *7«3. 
Can  you  add  lines  of  constant  pressure  to  Figure  *7-^? 
Note  that  for  a  perfect  gas  this  is  equivalent  to  a  T-s 
diagram.   In  which  direction  does  a  cooling  process  move 
along  the  subsonic  branch?   How  do  heating  and  cooling 
processes  move  on  the  supersonic  branch? 
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7- 3  -  Working  Relations  for  Raylegh  Flow  -  Perfect  Gas 

a.  We  now  proceed  to  develop  property  relations  of  a  function 
of  Mach  number  and  specific  heat  ratio.   These  will  be  the 
working  equations  used  for  problem  solving.   We  introduce 
the  usual  perfect  gas  assumptions  as  we  have  done  many 
times  previously. 

b.  From  the  momentum  equation  (*7«2,7)  we  can  write: 

v2 
P  +  ^r-  =  <t  (*7.3,D 

fec 


but 


V2   =   M2a2   =   M2YgcRT  (*7-3,2) 


and 

P  -  ^  (*7.3,3) 

Show  that 

P(l  +  yM2)  =  <f  (*7.3,4) 

Thus  the  pressure  ratio  between  the  two  points  is 

P0     1  +  yM-,2 

/  =  ^  (*7.3,5) 

rl      1  +  YM2 

How  would  you  obtain  a  ratio  for 


pt2 

99 


Is  this  related  to  the  entropy  change  as  we  have 
previously  always  been  noting? 
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c.   We  saw  that  the  continuity  equation  could  be  expressed  as 


PV 


Constant 


Show  that  this  can  be  written  as: 


PM 


■  <t 


(*7.3,6) 


thus  between  two  points 


T, 


'  M 

2  ^2 

>  M 

1  Ul 


(•7.3,7) 


Combining  this  with  equation  ( *7 • 3  , 5  )  show  that 


T. 


1  +  yM^ 

1  +  YM2' 


2   M  2 

2 

2 
IVL 


(•7.3,8) 


d.   In  Unit  3  we  saw  that 


Tt   =   T  [1  +  '1-~   M2] 


(3.6)  in 

text 


Now  you  can  show  that 


T 


t2 


T. 


l  +  ym 


i 


1  +  YM2' 


M, 


M- 


1  +  1—$L   m  2 

x     2    2 

2   Ul 


(»7.3,9) 


e.   In  summary,  considering  any  two  points  in  the  flow  system, 
if  we  know  the  properties  at  point  (l)  and  one  property 
at  point  (2)  we  can  compute  all  other  conditions  at 
point  (2)  . 
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Can  we  predict  what  will  occur  downstream? 

Yes,  if  we  have  information  concerning  the  heat  transfer. 

From  the  energy  equation  we  know  that 


\    +  ql+2   =   ht2 


(•7.2,3) 


Thus  in  general    Q-i^p   =   hj- 


2 


and  for  perfect  gases: 


^2   "   Cp(Tt2  "  V 


(•7.3,10) 


If  we  know  the  rate  of  heat  transfer  we  can  easily  solve 
for  Tt2  and  then  use  equation  (*7.3»9)  to  obtain  M2 ;  and 
all  the  other  properties  are  easily  computed. 


The  above  developments  assume  cp  =  constant. 
In  some  cases  where  heat  transfer  rates  are  extremely 
high  with  large  temperature  changes  resulting,  Cp  may 
vary  enough  to  warrant  using  an  average  value  for  cp. 
If,  in  addition,  significant  variations  in  y  occur,  it 
will  be  necessary  to  return  to  the  basic  equations  and 
derive  new  working  relations,  treating  y    as  a  variable. 
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7.4 


Reference  State  and  Tables 


a.   We  again  introduce  a  reference  where  the  Mach  number  is 
unity  —  as  reached  by  the  process  under  consideration  — 
i.e.,  by  Rayleigh  flow.   We  designate  this  the  '*' 
condition.   Note,  in  Figure  *7.5>  how  this  point  is 
introduced  in  both  the  T-s  diagram  and  the  physical 
diagram,  for  a  subsonic,  heating  example. 


M<1 


(b   q  d) 


M  =  1 


Figure  *7.5 


Complete  the  following  diagrams  (Figure  *7«6)  for  a 
supersonic,  cooling  example. 


M>1 


6        / 


& 


Figure  *7.6 


b.   In  the  previous  equations,  let  point  (T)  be  the  reference 
(or  ' *') condition  and  let  (J)  be  any  other  point  in  the 
flow.   Then  . . . 


M]_  =  1.0 
P1  =  P* 

ptl=  Pt*'  Tt] 


T1  =  t* 


M2  =  M 
P2  =  P 

Ptt-  P 


t> 


T2  =  T 
Tt2=  Tt 


16>J 


So 

p      pir  2  l  '<q'ij 

1  1  +  yM 

Derive  the  expression  for  T/T*  : 


T2     T 
*7 


T*       ?  (*7.4,2) 


In  like  manner,  expressions  for  the  other  ratios  can  be 
obtained  and  all  of  these  are  found  tabulated  in  Appendix 
P  in  the  text. 


c.   Pick  a  value  for  'M'  and  let  y    be  1.4.   Compute  the 

corresponding  value  of  T/T*  and  check  your  value  in  the 
tables . 


d.   Now  read  Section  10.2  and  study  the  examples  carefully. 
We  will  cover  the  remainder  of  Section  10  in  the  next 
unit  which  deals  with  the  correlation  of  Rayleigh  and 
Fanno  flow. 


e.   Note  that  in  Figure  10.5  in  the  text,  the  mass  flow  rate 
adjusts  to  a  new  Rayleigh  line  if  heat  is  added  to  a  duct 
in  which  the  flow  is  "choked." 


f.   It  will  prove  profitable  to  prepare  a  summary  of  the 
property  changes  that  are  taking  place  in  the  subsonic 
and  supersonic  regimes  for  the  heating  and  cooling 
processes.   Try  to  fill  in  the  spaces  in  Table  *7.1 
without  reference  to  the  tables  in  Appendix  F.   Use  basic 
equations,  special  equations  and/or  diagrams  that  have 
been  developed  for  Rayleigh  flow.   For  each  property, 
indicate  whether  that  property  is  increasing,  decreasing 
or  remaining  constant  for  each  regime. 
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TABLE    *7.1 


Property 

Heating 

Cooling 

M  <  1 

M  >  1 

M  <  1 

M>  1 

Stagnation 
Enthalpy 
(and  temperature 
if  a  gas) 

Entropy 

Static 
Pressure 

Velocity 

Mach  Number 

Static 
Enthalpy 
(and  temperature 
if  a  gas) 

Stagnation 
Pressure 
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Unit  8  -  CORRELATION  OP  FLOWS  -  OBJECTIVES 


The  student  shall  be  able  to: 

1.  Compare  similarities  and  differences  among  Fanno 
Flow,  Rayleigh  Flow  and  normal  shocks. 

2.  Sketch  on  the  same  h-s  diagram  a  Rayleigh  Line,  a 
Fanno  Line  and  a  normal  shock  process  (all  for  the 
same  flow  per  unit  area). 

3.  Solve  typical  problems  involving  various  sequential 
combinations  of  Fanno  Flow,  Rayleigh  Flow  and/or 
normal  shocks. 
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Unit  8  -  CORRELATION  OF  FLOWS  -  STUDY  GUIDE 

8. 1   Introduction. 

a)  In  the  previous  developments,  we  have  utilized  the 
basic  equations  relating  to  fluid  flow  and  have  de- 
veloped relations  which  have  considered  the  effects 
of  only  one  of  these  'factors'  at  a  time.   V/e  can 
gain  a  great  deal  more  flexibility  in  problem  solving 
by  considering  various  sequential  combinations  of 
Fanno  Flow,  Rayleigh  Flow  and  Normal  Shocks. 

b)  Another  type  of  problem  involves  simultaneous  addition 
(or  removal)  of  heat  with  friction  effects  of  the 
same  order  of  magnitude,  and  may  even  also  include 
changes  in  area.   Exact  treatments  for  problems  of  this 
type  usually  involve  tedious  numerical  integration 

as  explicit  solutions  can  not  be  obtained.   Your  text 
touches  on  some  of  these  problems  in  sections  9.^, 
9.5  and  10.5.   You  are  encouraged  to  read  these  sections 
after  you  have  completed  this  unit . 

c)  In  this  unit  we  wish  to  bring  out  some  of  the  similar- 
ities and  differences  among  Fanno  Flow,  Rayleigh  Flow 
and  Normal  Shocks.   Thus  we  shall  assume  steady,  one- 
dimensional  flow  in  a  constant  area  duct  with  negligi- 
ble changes  in  potential  and  no  shaft  work. 
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8.2  -  Fanno  Flow  with  a  Normal  Shock 


a)   Recall  from  Unit  7  the  definition  of  the  "Impulse" 
or  "Thrust"  Function: 


PA  + 


m  V 


(•7.2,7) 


We  saw  that  in  a  shock  the  two  states  (before  and 
after  the  shock)  displayed  the  same  flow  per  unit  area, 
the  same  stagnation  enthalpy  and  the  same  impulse 
function.   To  corroborate  this,  re-examine  equations 
*4.2,2;  *4.3}3  and  *4.2,6. 


b)   A  Fanno  Line  represents  states  with  the  same  flow  per 

unit  area  (m/A),  and  the  same  'stagnation  enthalpy'  (h^ ) . 
As  the  flow  passes  along  a  Fanno  Line,  it  is  experiencing 

the  effects  of  ?   Because  of  this  the  'Impulse 

Function'  does  not  remain  constant.   Suppose  we  now  plot 
the  'Impulse  Function'  versus  Mach  number  for  a  Fanno 
Line.   Such  a  plot  is  pictured  in  Figure  *"8.1. 


Mrl.O 


Figure  *8.1 


M 


c)   Note  that  for  every  point  on  the  supersonic  branch, 
there  is  a  corresponding  point  on  the  subsonic  branch 
where  the  Impulse  Function  has  the  same  value.   It  should 
be  obvious  that  points   A   and   B   satisfy  all  criteria 
for  not  only  the  Fanno  Line  but  also  for  the  Normal  Shock. 


d)   Thus  we  can  imagine  a  shock  between  points   A   and   B 
in  Figure  *8.1.   Would  this  shock  process  occur  from 
A  to   B   or  from  B  to   A  ? 
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e)   Consider  a  constant  area  duct  such  as  in  Figure  *8.2a 
which  is  experiencing  Fanno  Flow.   We  can  now  visualize 
the  possibility  of  a  shock  in  the  duct  as  pictured 
between  points  (2)  and  (3)  .   This  situation  meets  the 
conditions  of  flow  along  the  supersonic  branch  of  a 
given  Fanno  Line  from  point  (T)  to  point  \^2)    and  flow 
along  the  subsonic  branch  of  the  same  Fanno  Line  from 
(3)  to  (4)  .   We  can  also  plot  this  situation  in  the 
h-s  plane  as  shown  in  Figure  *8.2b. 


FloW' 


=£>   M>  1 


(D 


® 


M<  1 


(i) 


i) 


Figure  *8.2a 


h 


Figure  *8.2b 
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8.3-  Rayleigh  Flow  with  a  Normal  Shock 


a)   A  Rayleigh  Line  represents  states  which  exhibit  the 
same  flow  per  unit  area  and  the  same  impulse  function 
but  do  not  have  a  constant  stagnation  enthalpy,  due  to 
heat  transfer.   We  can  picture. a  Rayleigh  Line  and  its 
corresponding  Stagnation  Enthalpy  line  in  the  h-s  plane 
as  shown  in  Figure  *8.3« 


h 


Total 

(or  Stagnation) 

conditions 


Rayleigh  Line 

Figure  *8.3 


Notice  that  for  every  point  on  the  supersonic  branch 
of  the  Rayleigh  Line  there  is  a  corresponding  point  on 
the  subsonic  branch  which  has  the  same  stagnation 
enthalpy.   (Also  note  that  the  converse  is  not  always 
true).   Would  q2->-*  be  greater  than,  less  than  or  equal 

to  ^o^. . ?  Why  (°r  w^y  n°t )  ? 

b)  As  with  the  Fanno  Flow  and  a  Normal  Shock,  we  can  now 
surmise  the  existence  of  a  duct  experiencing  Rayleigh 
Flow  and  a  Normal  Shock  forming  as  pictured  in  Figure 
*8.4a.  The  corresponding  process  is  shown  in  the  h-s 
plane  in  Figure  *8.4b.  Note  that  the  flow  from  point  (T) 
to  point  (2)  and  the  flow  from  points  (J)  to  (4)  lie  on 
the  same  Rayleigh  Line. 
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Rayleigh 
Flow 


Flow 


=£> 


M>  1 

/ 


<b 


0 


Rayleigh 
Flow 


M<  1 


+ 


<b 


q- 


(4) 


Figure    *8.4a 


h 


Figure    (8.4b 
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8.  H      Rayleigh  Flow,  Fanno  Flow  and  Normal  Shocks 


a) 


Consider  nov;  the  combination  of  Rayleigh  Flov;,  Fanno 
Flow  and  a  shock  process  all  for  the  same  mass  flow 
per  unit  area,  as  pictured  in  the  h-s  plane  in  figure 
*8.5. 


Fanno  Line 


h 


Rayleigh  Line 


Figure  *8.5 


b)  We  have  seen  that  points  (T)  and  (2)  have  the  same 
mass  flow  rate,  the  same  impulse  function,  and  the  same 
stagnation  enthalpy;  and  thus  can  lie  on  a  Fanno  Line, 

a  Rayleigh  Line,  and  at  opposite  sides  of  a  Normal  Shock 

c)  From  the  units  on  Rayleigh  and  Fanno  flows  we  saw 

that  when  the  "choked  flov;"  condition  was  reached  along 
the  subsonic  flow  branch,  the  mass  flow  rate  re- 
adjusted to  provide  for  flov;  along  a  new  Rayleigh 
Line  (or  Fanno  Line)  if  more  heat  (or  duct  length, 
hence  friction)  was  added  to  the  flow.   In  the  super- 
sonic regime  the  readjustment  can  (and  usually  does) 
take  the  form  of  a  normal  shock.   In  this  case  we  do 
not  necessarily  have  to  move  to  a  new  Fanno  or  Rayleigh 
Line  since  once  subsonic  flow  has  been  obtained  we  can 
add  much  more  heat,  or  introduce  much  more  friction, 
as  the  case  may  be.   Many  problems  can  be  considered 
by  matching  these  three  types  of  flov;  processes 
sequentially,  knowing  the  flow  conditions.   Let  us  now 
consider  qualitatively  two  such  possible  combinations. 

d)  Consider  the  combination  in  figure  *8.6a.   Assume  that 
only  when  heat  is  being  transferred  may  friction 
effects  be  considered  negligible. 
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Flow 


(b 


Fanno 
Flow 


=> 


Rayleigh 
Flow 


(a) 


<\ 


Fanno 
Flow 


® 


Figure  *8.6a 


® 


If  this  entire  flow  system  is  in  the  subsonic  regime 
the  h-s  diagram  would  be  that  shown  in  Figure  *8.6b. 


h 


Figure  *8.6b 


We  may  describe  the  process  as  follows: 

1)  From  (T)  to  (2)  the  flow  is  along  a  Fanno  Line 
hence  ' h-^'  and  'm/A'  are  constant. 

2)  From  (2)  to  (3)  'm/A'  remains  unchanged  as  does 
'  'F'  (i.e.,  'PA  +  mV/gc'),  but  due  to  the  heat 

transfer,  'h^'  is  increasing. 

3)  From  (3)  to  (4)  the  flow  is  along  a  different 
Fanno  Line  even  though  it  has  the  same  'm/A'. 
The  reason  is  that  the  stagnation  enthalpies  for 
the  two  Fanno  Lines  are  different. 


e. 


Now  suppose  we  consider  the  situation  in  Figure  *8.6a 
but  with  supersonic  flow  coming  in  at  point  (T)  .   Plot 
this  situation  in  the  h-s  plane  in  Figure  *8.6c  assuming 
that  the  flov;  is  not  choked. 
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Figure 
*8.6c 


f)   V/e  can  easily  imagine  the  amount  of  heat  added  from 

(2)  to  (3)  exceeding  that  required  to  "choke"  the  flow, 
in  which  case  a  shock  could  possibly  form.   Many 
interesting  possibilities  could  result  depending  upon 
the  boundary  conditions  imposed  on  the  system.   Are  the 
original  stagnation  conditions  maintained?   If  the  flow 
is  subsonic  at  the  exit,  what  pressure  must  exist? 


g)   It  might  prove  useful  to  summarize  the  effect  on  fluid 
properties  of  flow  in  ducts  in  the  supersonic  and 
subsonic  regimes  for  the  various  types  of  flow. 


Table  *8.1 


Property 

Rayleigh  Plow 

Fanno  Flow 

Normal 
Shock 

M>1 

Heating 

Coo: 

.ing 

M<1 

M>1 

M<1 

M>1 

M<1 

M>1 

Pressure 

D 

I 

I 

D 

D 

I 

I 

Mach 

I 

D 

D 

I 

I 

D 

to  M<1 

Velocity 

I 

D 

D 

I 

I 

D 

D 

Temperature 

I/D 

I 

I/D 

D 

D 

I 

I 

Entropy 

I 

I 

D 

D 

I 

I 

I 

Stag.  Press. 

D 

D 

I 

I 

D 

D 

D 

Stag.  Temp. 

I 

I 

D 

D 

Density 

D 

I 

I 

D 

D 

I 

I 

Impulse  Fct. 

D 

D 

Check  this  table  with  the  information  developed  in  the 
previous  units  and  satisfy  yourself  that  the  information 
is  correct. 

h)   Mow  read  Section  10.3  in  the  text.   If  time  permits, 
glance  at  Sections  9  •  *J ,  9-5  and  10.5. 
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APPENDIX  B  - 

STUDENT  EVALUATION  QUESTIONAIRE + 

AE  30^3,  "Gas  Dynamics"  -  Quarter  IV,  1971-72 

Please  answer  these  questions  as  completely  and  candidly 
as  possible.   The  results  will  in  NO  fashion  effect  your 
grade,  but  will  be  utilized  in  assessing  the  value  and 
desirability  of  this  method  of  instruction. 

Where  you  are  given  a  scale  on  which  to  indicate  your  opinion, 
a  descriptive  work  (or  words)  is  provided  to  define  each 
extreme.   Please  consider  the  scale  linear  with  the  mean 
midway  between  2  and  3-   You  may  thus  circle  1 ,  2 ,  3  or  4 
to  indicate  your  own  attitude  relative  to  the  extremes. 
Any  amplifying  comments  you  may  wish  to  add  to  any  question 
will  be  greatly  appreciated.   (Use  the  back  of  the  question- 
aire  or  attach  additional  pages,  as  desired) 

1.   How  many  previous  courses  (of  any  duration,  here  at  NPS 
or  elsewhere)  have  you  taken  which  were... 


SELF  STUDY  |Ave .  1 .  91  ?   SELF  PACED   lAve .  2.91  ? 

PROGRAMMED  175%  had  had  many]  ? 

Indicate  your  opinion  of  the  utility  of  the  following 
educational  innovations  (NOTE:  These  features  may  or 
may  NOT  have  been  included  in  this  course.) 


a. 

Defined 

Obj  ectives 

(no 

value) 

1 

2 

3 

1] 

(great 

value) 

b. 

Self-Study 

Concept 

(no 

value) 

1 

2 

3 

4 

(great 

value ) 

c . 

Self-Paced 

Instruction 

(no 

value ) 

1 

2 

3 

4 

(great 

value ) 

d. 

Study  Guide 

(no 

value ) 

1 

2 

3 

4 

(great 

value ) 

e . 

Check  Tests 

(no 

value) 

1 

2 

3 

i| 

(great 

value ) 

f. 

Group 

Discussions 

(no 

value ) 

1 

2 

3 

A 

(great 

value) 

g- 

Individual 

Tutoring 

(no 

value ) 

1 

2 

3 

k 

(great 

value ) 

h. 

Audio- 

Tutorial 

(no 

value) 

1 

2 

3 

k 

(great 

value) 

(Tape-Slides) 

i. 

Final  Exam 

(no 

value) 

1 

2 

3 

H 

(great 

value ) 

Boxed  values  indicate  the  average  results  of  Questionnaires 
returned  to  the  writer. 
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STUDENT  EVALUATION  QUESTIONAIRE  (page  2.) 


The  following  questions  pertain  to  the  materials/methods 
used  in  this  course.   Please  consider  each  question  as  it 
reflects  on  the  course  as  a  whole,  then  circle  or  underline 
your  opinion.   Then  below  each  question,  identify  any  unit 
(or  section)  to  which  your  general  answer  does  not  apply  and 
explain  why.   You  will  probably  find  it  useful  to  have  your 
course  materials  at  hand  while  considering  these  questions. 


3. 


4. 


5. 


6. 


7. 


How  would  you  rate  the  clarity  of  the  Course  OBJECTIVES? 

(obscure)    1   2   3   ^   (clear) 
Exceptions : 


Average  rating  -  3.58 


Did  you  find  the  Study  Guides  understandable? 

(obscure)    1   2   3   ^   (clear) 
Exceptions : 


Average  rating  -  3-51 


What  is  your  opinion  of  the  organization  of  the  Study 
Guides? 


(confused)   1 
Exceptions : 


2   3^   (logical) 


Average  ratins 


3.^2 


Were  sufficient  graphs/diagrams/etc.  provided  in  the 
Study  Guides? 

(far  too  few)   1   2   3   4   (far  too  many) 
Exceptions : 


Average  rating  -  2 . k2 


Were  there  enough  example  problems  provided  in  the  Study 
Guides? 

(far  too  few)   1   2   3   4   (far  too  many) 
Examples : 


Average  rating  -  2.17 
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STUDENT  EVALUATION  QUESTIONAIRE  (page  3.) 


8.   Were  there  enough  homework  problems  provided  to  cover 
the  material? 


(far  too  few)    12   3   4   (far  too  many) 
Exceptions : 


Average  rating  -  3-04 


9.   Was  a  correct  portion  of  the  quarter  allocated  for  the 

completion  of  each  unit  in  the  course: 

(too  little  time)   1   2   3   4   (too  much  time) 
Exceptions : 


Average  rating  -  2.58 


10.   Outside  of  the  scheduled  problem  sessions,  on  the  average, 
how  many  times  did  you  consult  the  instructor  for 
advice/assistance,  per  unit?   

Exceptions : 


Average  -  1  time 


11.   Did  you  seek  advice/assistance  from  your  classmates? 


(seldom) 
Exceptions 


12   3   4   (regularly) 


Average  rating  -  2.5 


12.   How  do  you  feel  the  course  materials  aided  you  in  pre- 
paring for  the  check  tests? 


(useless ) 
Exceptions : 


12   3   4   (invaluable) 


Average  rating  -  3-58 


13.   How  beneficial  were  the  check  tests  in  preparing  for 
the  major  quizzes/final  exam? 


(useless ) 
Exceptions : 


12      3      4      (invaluable) 


Average   rating  -    3. 08 
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STUDENT  EVALUATION  QUESTIOMAIRE  (page  4 . ) 


14 


15 


16 


17 


18 


19 


Estimate  the  percentage  of  your  total  study  time 
expended  on  this  course.   


Average  rating  -  25% 


Do  you  feel  any  unit  was  excessively  long  or  short? 
(If  so  please  identify  and  comment.) 


Unit  one  -  2 
Unit  four-  2 


feel  the  amount  you  learned  using  this 
self-study  method  compares  with  what  you 
would  have  learned  in  a  conventional  lecture  course? 


How  do  you 
controlled 


(less  in  this  course:  about  the  same;  more  in  this  course) 


1 


0 


H 


For  several  reasons  this  course  was  run  on  a 
control led -pace  basis. 


a. 


Would  you  have  liked  the  course  better  if  it  were 
completely  self-paced?  


Yes 
No 


b. 


If  it  were  run  on  a  completely  self-paced  basis, 
when  do  you  think  you  would  have  finished  the  course? 


on  time  -  7 
early  -  2 
late     -  2 


If  you  took  a  self-paced  course  and  you  were  among  the 
first  students  to  successfully  finish  a  unit,  would 
you  be  willing  to  serve  as  a  tutor  to  those  having 
more  difficulty?   (please  comment ) 


Yes  -  11 

No   -   1 


Assuming  check  tests  would  be  retaken  until  mastery  of 
material  is  essentially  demonstrated,  and  that  failure 
of  a  check  test  would  not  be  held  against  you  -  would 
it  be  better  to  have  more  comprehensive  check  tests  and 
eliminate  the  major  quizzes?   (please  comment) 


Yes 
No 


-  6 

-  6 
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STUDENT  EVALUATION  QUESTIONNAIRE  (page  5.) 
20.   Would  you  prefer  check  tests  to  be  .. 


closed  book? 

7 

open  book? 

0 

part  open, 

part  closed 

k 

oral? 

1 

21.   To  what  extent  has  participation  in  this  course  improved 

your  ability  to  "decipher"  a  typical  textbook  or  technical 
paper? 


little 
much 


22 


23 


Given  the  choice  of  some  type  of  self-study  method  as 
opposed  to  instruction  using  the  conventional  lecture 
method,  which  would  you  prefer? 

(You  may  qualify  your  answer  in  any  way  you  choose, 
but  please  give  your  reasons.  This  is  your  chance 
to  comment  on  anything  you  feel  this  questionnaire 
overlooked . ) 


self  study 
lecture 


9 
2 


Would  you  like  to  take  all  your  courses  by  some  type  of 
a  self-study  method? 


no 

yes 

undecided 


THANK  YOU  FOR  YOUR  TIME  AND  EFFORT! 
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APPENDIX  C  - 

SELECTED  COMMENTS  EXTRACTED  FROM  STUDENT  QUESTIONNAIRES 
AND  CRITIQUE  SESSION. 


Question  19 

"...  yes...  eliminate  the  pressure  of  having  to  do 
well  at  a  specified  time..." 

"Yes.  Retakes  under  present  check  test,  hour  quiz  policy 
are  senseless  waste  of  time." 

"No  because  it  is  very  tiresome  getting  up  for  a  test 
every  week" 

"No.   The  check  tests  are  aids,  but  the  quizzes  are  for 
grades,  (forgive  my  poetic  license)" 


Question  21 


Comment  by 

author: 

In 

general  the  students  seemed  to 

miss  the 

intent 

of 

this 

question . 

Several 

comment 

ed  that 

they 

had 

no 

basis 

on  which 

to  answer 

Most 

answers 

seemed 

hesitant 

,  or  non- 

commital . 

Question  22 


"Self  study  involving  lesson  guides,  (audio-visual  perhaps) 
Prob  sets 
Followed  by  active  participation  seminar  to  reinforce 
important  points" 

"Prefer  a  lecture.   I  have  a  good  memory  so  I  can  keep 
up  with  previous  ideas  while  learning  others.   Most 
lecturers  speak  at  an  easily  grasped  pace  for  me." 


c 

omment 

by 

author : 

9  of 

the 

12  que 

stionnaires 

indicated 

a 

prefe: 

rence 

for 

this 

method 

or 

some 

mo 

dification 

of 

self- 

-paced 

inst 

ruct 

ion. 

Comment  from  Critique  Session 

"...  This  was  really  not  what  I  would  call  self-taught 
but  was  more  like  an  ideal  lecture  course,  where 
everyone  has  read  his  lesson  assignment  before  class 
and  the  material  is  re-inforced  by  the  professor  in  the 
lecture ..." 
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